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Leptogenesis is a class of scenarios where the baryon asymmetry of the Universe is produced from a 
lepton asymmetry generated in the decays of a heavy sterile neutrino. We explain the motivation for 
leptogenesis. We review the basic mechanism, and describe subclasses of models. We then focus on recent 
developments in the understanding of leptogenesis: finite temperature effects, spectator processes, and 
' in particular the significance of flavour physics. 
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1 The Baryon Asymmetry of the Universe 



1.1 Observations 

Observations indicate that the number of baryons (protons and neutrons) in the Universe is unequal to 
the number of antibaryons (antiprotons and antineutrons). To the best of our understanding, all the 
structures that we see in the Universe - stars, galaxies, and clusters - consist of matter (baryons and 
electrons) and there is no antimatter (antibaryons and positrons) in appreciable quantities. Since various 
considerations suggest that the Universe has started from a state with equal numbers of baryons and 
antibaryons, the observed baryon asymmetry must have been generated dynamically, a scenario that is 
known by the name of baryogenesis. 

One may wonder why we think that the baryon asymmetry has been dynamically generated, rather 
than being an initial condition. There are at least two reasons for that. First, if a baryon asymmetry had 
been an initial condition, it would have been a highly fine-tuned one. For every 6,000,000 antiquarks, there 
should have been 6,000,001 quarks. Such a fine-tuned condition seems very implausible. Second, and 
perhaps more important, we have excellent reasons, based on observed features of the cosmic microwave 
background radiation, to think that inflation took place during the history of the Universe. Any primordial 
baryon asymmetry would have been exponentially diluted away by the required amount of inflation. 

The baryon asymmetry of the Universe poses a puzzle in particle physics. The Standard Model (SM) 
of particle interactions contains all the ingredients that are necessary to dynamically generate such an 
asymmetry in an initially baryon-symmetric Universe. Yet, it fails to explain an asymmetry as large as 
the one observed (see e.g. [1]). New physics is called for. The new physics must [2], first, distinguish 
matter from antimatter in a more pronounced way than do the weak interactions of the SM. Second, it 
should depart from thermal equilibrium during the history of the Universe. 

The baryon asymmetry of the Universe can be defined in two equivalent ways: 

= (6.21 ±0.16) x 10" 10 , (1.1) 

o 

= (8.75 ±0.23) x 10" 11 (1.2) 

o 

where tig, rig, n 7 and s are the number densities of, respectively, baryons, antibaryons, photons and 
entropy, a subscript implies "at present time" , and the numerical value is from combined microwave 
background and large scale structure data (WMAP 5 year data, Baryon Acoustic Oscillations and Type 
la Supernovae) [3]. It is convenient to calculate Yab, the baryon asymmetry relative to the entropy 
density s, because s = g*(27r 2 /45)T 3 is conserved during the expansion of the Universe (g* is the number 
of degrees of freedom in the plasma, and T is the temperature; see the discussion and definitions below 
cqn (|13.6p ). The two definitions (II. ip and (|1.2p are related through Yab = (n-yo / 'so)t) — ry/7.04. A third, 
related way to express the asymmetry is in terms of the baryonic fraction of the critical energy density, 

= ps I petit- (1.3) 

The relation to r\ is given by 

f] = 2.74 x 10" 8 n B h 2 , (1.4) 

where h = H Q /100 km s^ 1 Mpc -1 = 0.701 ± 0.013 [4] is the present Hubble parameter. 

The value of baryon asymmetry of the Universe is inferred from observations in two independent ways. 
The first way is via big bang nucleosynthesis [5-7] . This chapter in cosmology predicts the abundances 
of the light elements, D, 3 He, 4 He, and 7 Li. These predictions depend on essentially a single parameter 
which is 77. The abundances of D and 3 He are very sensitive to 77. The reason is that they are crucial in 
the synthesis of 4 He via the two body reactions D(p, 7) 3 He and 3 He(D,p) 4 He. The rate of these reactions 
is proportional to the number densities of the incoming nuclei which, in turn, depend on r\: n(D) cx rj 
and n( 3 He) cx i] 2 . Thus, the larger 77, the later these 4 He-producing processes will stop (that is, become 
slower than the expansion rate of the Universe), and consequently the smaller the freeze-out abundances 
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of D and of 3 He will be. The abundance of 4 He is less sensitive to 77. Larger values of 77 mean that the 
relative number of photons, and in particular photons with energy higher than the binding energies of the 
relevant nuclei, is smaller, and thus the abundances of D, 3 He and 3 H build up earlier. Consequently, 4 He 
synthesis starts earlier, with a larger neutron-to-proton ratio, which results in a higher 4 He abundance. 
The dependence of the Li-abundance on 77 is more complicated, as two production processes with opposite 
77-dependencies play a role. 

The primordial abundances of the four light elements can be inferred from various observations. The 
fact that there is a range of 77 which is consistent with all four abundances gives a strong support to the 
standard hot big bang cosmology. This range is given (at 95% CL) by [5] 

4.7 x 1(T 10 < 77 < 6.5 x 1(T 10 , 0.017 < VL B h 2 < 0.024. (1.5) 

The second way to determine Q B is from measurements of the cosmic microwave background (CMB) 
anisotropics (for a pedagogical review, see [8,9]). The CMB spectrum corresponds to an excellent ap- 
proximation to a blackbody radiation with a nearly constant temperature T. The basic observable is the 
temperature fluctuation 0(n) = AT/T (n denotes the direction in the sky). The analysis is simplest in 
Fourier space, where we denote the wavenumber by k. 

The crucial time for the CMB is that of recombination, when the temperature dropped low enough 
that protons and electrons could form neutral hydrogen. This happened at redshift z roc ~ 1000. Before 
this time, the cosmological plasma can be described, to a good approximation, as a photon-baryon fluid. 
The main features of the CMB follow from the basic equations of fluid mechanics applied to perfect 
photon-baryon fluid, neglecting dynamical effects of gravity and the baryons: 

e + c 2 s k 2 e = 0, c s = *JJfp = (1.6) 

where c s is the sound speed in the dynamically baryon-free fluid {p and p are the photon energy density 
and pressure). These features in the anisotropy spectrum are: the existence of peaks and troughs, the 
spacing between adjacent peaks, and the location of the first peak. The modifications due to gravity and 
baryons can be understood from adding their effects to eqn (|1.6p . 

e + clk 2 e = F, Ca = 1 (1.7) 

a/3(1 + 3pb/4/j 7 ) 

where F is the forcing term due to gravity, and p B is the baryon energy density. The physical effect of the 
baryons is to provide extra gravity which enhances the compression into potential wells. The consequence 
is enhancement of the compressional phases which translates into enhancement of the odd peaks in the 
spectrum. Thus, a measurement of the odd/even peak disparity constrains the baryon energy density. A 
fit to the most recent observations (WMAP5 data only, assuming a ACDM model with a scale-free power 
spectrum for the primordial density fluctuations) gives (at 2 a) [4] 

0.02149 < fl B h 2 < 0.02397. (1.8) 

The impressive consistency between the nucleosynthesis (jl.5p and CMB (|1.8p constraints on the baryon 
density of the Universe is another triumph of the hot big-bang cosmology. A consistent theory of baryo- 
genesis should explain n B ~ 10~ 10 s. 



1.2 Ingredients and Mechanisms 

The three ingredients required to dynamically generate a baryon asymmetry were given by Sakharov in 
Ref. [2]: 

1. Baryon number violation: This condition is required in order to evolve from an initial state with 
Yab = to a state with Yab 7^ 0. 
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2. C and CP violation: If cither C or CP were conserved, then processes involving baryons would 
proceed at precisely the same rate as the C- or CP-conjugate processes involving antibaryons, with 
the overall effects that no baryon asymmetry is generated. 

3. Out of equilibrium dynamics: In chemical equilibrium, there are no asymmetries in quantum num- 
bers that are not conserved (such as B, by the first condition). 

These ingredients are all present in the Standard Model. However, no SM mechanism generating a 
large enough baryon asymmetry has been found. 

1. Baryon number is violated in the Standard Model, and the resulting baryon number violating 
processes are fast in the early Universe [10]. The violation is due to the triangle anomaly, and leads 
to processes that involve nine left-handed quarks (three of each generation) and three left-handed 
leptons (one from each generation). A selection rule is obeyed, 4 

AB = AL = ±3. (1.9) 

At zero temperature, the amplitude of the baryon number violating processes is proportional to 
e _87r 1 9 [11], which is too small to have any observable effect. At high temperatures, however, 
these transitions become unsuppressed [10]. 

2. The weak interactions of the SM violate C maximally and violate CP via the Kobayashi-Maskawa 
mechanism [12]. This CP violation can be parameterized by the Jarlskog invariant [13] which, when 
appropriately normalized, is of order 10~ 20 . Since there are practically no kinematic enhancement 
factors in the thermal bath [14-16], it is impossible to generate Yab ~ 10~ 10 with such a small 
amount of CP violation. Consequently, baryogenesis implies that there must exist new sources of 
CP violation, beyond the Kobayashi-Maskawa phase of the Standard Model. 

3. Within the Standard Model, departure from thermal equilibrium occurs at the electroweak phase 
transition [17,18]. Here, the non-equilibrium condition is provided by the interactions of particles 
with the bubble wall, as it sweeps through the plasma. The experimental lower bound on the Higgs 
mass implies, however, that this transition is not strongly first order, as required for successful 
baryogenesis. Thus, a different kind of departure from thermal equilibrium is required from new 
physics or, alternatively, a modification to the electroweak phase transition. 

This shows that baryogenesis requires new physics that extends the SM in at least two ways: It must 
introduce new sources of CP violation and it must either provide a departure from thermal equilibrium 
in addition to the electroweak phase transition (EWPT) or modify the EWPT itself. Some possible new 
physics mechanisms for baryogenesis are the following: 

GUT baryogenesis [19-28] generates the baryon asymmetry in the out-of-equilibrium decays of 
heavy bosons in Grand Unified Theories (GUTs). The Boltzmann Equations (BE) for the bosons and 
the baryon asymmetry are studied, for instance, in [29-31]. BE are also required for leptogenesis, and in 
this review, we will follow closely the analysis of Kolb and Wolfram [31]. The GUT baryogenesis scenario 
has difficulties with the non-observation of proton decay, which puts a lower bound on the mass of the 
decaying boson, and therefore on the reheat temperature after inflation 5 . Simple inflation models do not 
give such a high reheat temperature, which in addition, might regenerate unwanted relics. Furthermore, 
in the simplest GUTs, B + L is violated but B — L is not. Consequently, the B + L violating SM 
sphalerons, which are in equilibrium at T 10 12 GeV, would destroy this asymmetry 6 . 

Leptogenesis was invented by Fukugita and Yanagida in Ref . [35] . New particles - singlet neutrinos 
- are introduced via the seesaw mechanism [36-40]. Their Yukawa couplings provide the necessary new 

4 This selection rule implies that the sphaleron processes do not mediate proton decay. 

5 Bosons with m > Tr^^t could nonetheless be produced during "preheating" [32,33], which is a stage between the end 
of inflation and the filling of the Universe with a thermal bath. 

6 A solution to this problem [34], in the seesaw model, could be to have fact L violation (due to the righthanded neutrinos) 
at T > 10 GeV. This would destroy the L component of the asymmetry, and the remaining B component would survive 
the sphalerons. 
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source of CP violation. The rate of these Yukawa interactions can be slow enough (that is slower than 

H, the expansion rate of the Universe, at the time that the asymmetry is generated) that departure from 
thermal equilibrium occurs. Lepton number violation comes from the Majorana masses of these new 
particles, and the Standard Model sphalcron processes still play a crucial role in partially converting the 
lepton asymmetry into a baryon asymmetry [41]. This review focuses on the simplest and theoretically 
best motivated realization of leptogenesis: thermal leptogenesis with hierarchical singlet neutrinos. 

Electroweak baryogenesis [17,42,43] is the name for a class of models where the departure from 
thermal equilibrium is provided by the electroweak phase transition. In principle, the SM belongs to this 
class, but the phase transition is not strongly first order [44] and the CP violation is too small [14, 15]. 
Thus, viable models of electroweak baryogenesis need a modification of the scalar potential such that 
the nature of the EWPT changes, and new sources of CP violation. One example [45] is the 2HDM 
(two Higgs doublet model), where the Higgs potential has more parameters and, unlike the SM potential, 
violates CP. Another interesting example is the MSSM (minimal supersymmetric SM), where a light 
stop modifies the Higgs potential in the required way [46, 47] and where there are new, flavour-diagonal, 
CP violating phases. Electroweak baryogenesis and, in particular, MSSM baryogenesis, might soon be 
subject to experimental tests at the CERN LHC. 

The Affleck-Dine mechanism [48,49]. The asymmetry arises in a classical scalar field, which later 
decays to particles. In a SUSY model, this field could be some combination of squark, Higgs and slepton 
fields. The field starts with a large expectation value, and rolls towards the origin in its scalar potential. 
At the initial large distances from the origin, there can be contributions to the potential from baryon or 
lepton number violating interactions (mediated, for instance, by heavy particles). These impart a net 
asymmetry to the rolling field. This generic mechanism could produce an asymmetry in any combination 
of B and L. 

Other, more exotic scenarios, are described in Ref. [1]. 

I. 3 Thermal leptogenesis: advantages and alternatives 

There are many reasons to think that the Standard Model is only a low energy effective theory, and 
that there is new physics at a higher energy scale. Among these reasons, one can list the experimental 
evidence for neutrino masses, the fine-tuning problem of the Higgs mass, the dark matter puzzle, and 
gauge coupling unification. 

It would be a particularly interesting situation if a new physics model, motivated by one of the reasons 
mentioned above, would also provide a viable mechanism for baryogenesis. Leptogenesis [35] is such a 
scenario, because it is almost unavoidable when one invokes the seesaw mechanism [36-40] to account 
for the neutrino masses and to explain their unique lightness. Indeed, the seesaw mechanism requires 
that lepton number is violated, provides in general new CP violating phases in the neutrino Yukawa 
interactions, and for a large part of the parameter space predicts that new, heavy singlet neutrinos decay 
out of equilibrium. Thus, all three Sakharov conditions are naturally fulfilled in this scenario. The 
question of whether leptogenesis is the source of the observed baryon asymmetry is then, within the 
seesaw framework, a quantitative one. 

In the context of the seesaw extension of the SM, there are several ways to produce a baryon asymme- 
try. They all have in common the introduction of singlet neutrinos Ni with masses Mi that arc usually 
7 heavier than the electroweak breaking scale, Mi ^> v u . They may differ, however, in the cosmologi- 
cal scenario, and in the values of the seesaw parameters. (The number of seesaw parameters is much 
larger than the number of measured, light neutrino parameters.) A popular possibility, which this review 
focuses on, is "thermal leptogenesis" with hierarchical masses, Mi <C Mj > i [35]. The Ni particles are 
produced by scattering in the thermal bath, so that their number density can be calculated from the 
seesaw parameters and the reheat temperature of the Universe. 

Thermal leptogenesis has been studied in detail by many people, and there have been numerous clear 
and pedagogical reviews. Early analyses, focusing on hierarchical singlet neutrinos, include [51-54] (see 
the thesis [55] for details, in particular of supersymmetric leptogenesis with superficlds). The importance 

7 An exception is the ^MSM [50] . 
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of including the wave function renormalisation of the decaying singlet, in calculating the CP asymmetry, 
was recognised in [56]. Various reviews [57,58] were written at this stage, a pedagogical presentation 
that introduces BE, and discusses models and supersymmetry can be found in [58]. More detailed 
calculations, but which do not include flavour effects, are presented in [59,60]: thermal effects were 
included in [60], and [59] gives many useful analytic approximations. The importance of flavour effects 
was emphasized in [61] for resonant leptogenesis with degenerate N t , and in [62-64] for hierarchical 
singlets. An earlier "flavoured" analysis is [65], and flavoured BE are presented in [66]. The aim of this 
review is to pedagogically introduce flavour effects, which can change the final baryon asymmetry by 
factors of few or more. Some previous reviews of flavour effects can be found in the TASI lectures of [67], 
and in the conference proceedings [68-70]. They are also mentioned in Refs. [71,72]. 

A potential drawback of thermal leptogenesis, for hierarchical masses M^'s, is the lower bound on 
Mi [73,74] (discussed in Section IS~2|) . which gives a lower bound on the reheat temperature. This bound 
might be in conflict with an upper bound on the reheat temperature that applies in supersymmetric 
models with a "gravitino problem" [75-81]. The lower bound on Mi can be avoided with quasi-degenerate 
Mi's [56,61,82-85], where the CP violation can be enhanced in iV,- — Nj mixing. This scenario is discussed 
in Section [101 

Other leptogenesis scenarios, some of which work at lower reheat temperatures, include the following: 

• "Soft leptogenesis" [86-88] , which can work even in a one-generation SUSY seesaw. Here the source 
of both lepton number violation and CP violation is a set of soft supersymmetry breaking terms. 

• The Affleck Dine mechanism [48,49,89] where the scalar condensate carries lepton number. For a 
detailed discussion, see, for instance, [89]. 

• The Ni's could be produced non-thermally, for instance in inflaton decay [90,91], or in preheating 
[92,93]. 

• The singlet sneutrino N could be the inflaton, which then produces a lepton asymmetry in its 
decays [94]. 

• Models with a (e.g. flavour) symmetry that breaks below the leptogenesis scale, and involve addi- 
tional heavy states that carry lepton number [95]. 

1.4 Reading this review 

This review contains chapters at different levels of detail, which might not all be equally interesting to 
all readers. Analytic formulae, which can be used to estimate the baryon asymmetry, can be found in 
Appendix 1161 Appendix 1121 is a dictionary of the notation used in this review. A reader interested in a 
qualitative understanding and wishing to avoid Boltzmann Equations (BE), can read Section H perhaps 
Section [3[ Sections [4] 15.11 and 15. 2i and browse Sections [9] and [101 

For a more quantitative understanding of thermal leptogenesis, some acquaintance with the BE is 
required. In addition to the Sections listed above, readers may wish to review preliminary definitions in 
Appendix [T21 read the introduction to simple BE of Sections 16.11 and 16.21 and browse Sections [7] and [5J 
and Appendix [T4l 

Aficionados of leptogenesis may also be interested in Section 16 - 31 in which more complete BE are 
derived, and Appendix [TBI that describes the toy model which motivates our claims about flavour effects. 

The layout of the review is the following: In Section [2] we briefly review the seesaw mechanism and 
describe various useful parameterizations thereof. Section [3] introduces the non-perturbative B + L vio- 
lating interactions which are a crucial ingredient for baryogenesis through leptogenesis. Section 3] gives 
an overview of leptogenesis, using rough estimates to motivate the qualitative behaviour. Flavour depen- 
dent quantities (CP asymmetries and efficiencies) are introduced already in this section: flavour blind 
equations in section [4] would be inapplicable later in the review. The calculation of the CP asymmetry 
and bounds thereon are presented in Section [5[ Subsections 15.31 and 15.41 discuss the implications of CPT 
and unitarity, relevant to a reader who wishes to derive the BE. Simple BE that describe interactions 
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mediated by the neutrino Yukawa coupling are derived in Section [6.21 using definitions from Section [6. II 
(Formulae for the number densities and definitions of the rates, calculated at zero temperature, can be 
found in Appendix [T3"l ) In Section |6~31 the most relevant scattering interactions, and a large set of 0(1) 
effects, are included. The next three sections refine the analysis by implementing various ingredients 
that are omitted in the basic calculation. Section [7] is an overview of finite temperature effects. Sec- 
tion [5] analyzes the impact of spectator processes. (These are often described by chemical equilibrium 
conditions, which are reviewed in Appendix [Til ) Section [5] is devoted to a detailed discussion of flavour 
effects. In section 1101 we present some variations on the simplest framework: non-hierarchical singlets 
(subsection 110. 2[) , soft leptogenesis (subsection 110. 3[) , Dirac leptogenesis (subsection 110. 4p , leptogenesis 
with scalar triplets (subsection 1 1 . 5|) . and with fermion triplets (subsection 110. 6|) . Our conclusions and 
some prospects for future developments are summarized in Section 111! We collect all the symbols that 
are used in this review into a table in Appendix 1121 where we also give the equation number where they 
are defined. The basic ingredients and consequences of kinetic equilibrium are reviewed in Appendix 1131 
while those of chemical equilibrium are reviewed in Appendix [TJJ The BE used in the main text are 
not covariant under transformations of the lepton flavour basis. Appendix [15] discusses, via a toy model, 
covariant equations for a "density matrix" , and explains how a flavour basis for the BE is singled out by 
the dynamics. Approximate solutions to simple BE are given in Appendix 1161 
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2 Seesaw Models 



Measurements of fluxes of solar, atmospheric, reactor and accelerator neutrinos give clear evidence that 
(at least two of) the observed weakly interacting neutrinos have small masses. Specifically, two neutrino 
mass-squared differences have been established, 

Ato^ = (7.9 ± 0.3) x 10~ 5 eV 2 = m 2 ol , |Am^ 2 | = (2.6 ± 0.2) x 10" 3 eV 2 = m 2 tm . (2.1) 

If these are "Majorana" masses, then they violate lepton number and correspond to a mass matrix of the 
following form: 

An„ = -jV c a [m] a i3Vf3 +h.c. (2.2) 

The mass terms [m] a| a break the SU(2)l gauge symmetry as a triplet. Consequently, they cannot 
come from renormalizable Yukawa couplings with the Standard Model Higgs doublet. (This stands 
in contrast to the charged fermion masses which break SU(2)l as doublets and are generated by the 
renormalizable Yukawa couplings.) Instead, they are likely to arise from the dimension five operator 
(£ a 4>)(£i3(j)). Seesaw models are full high energy models which induce this effective low energy operator at 
tree level. Seesaw models are attractive because they naturally reproduce the small masses of the doublet 
neutrinos. Explicitly, if the exchanged particle has a mass M which, by assumption, is much heavier than 
the electroweak breaking scale v u , then the light neutrino mass scale is v 2 /M, which is much smaller 
than v u , the charged fermion mass scale. 

There are three types of seesaw models, which differ by the properties of the exchanged heavy particles: 

• Type I: SU(S) x SU(2) x [/(l)-singlet fermions; 

• Type II: 5t/(2)-triplet scalars; 

• Type III: 5J7(2)-triplet fermions. 

We now describe these three types of seesaw models in more detail, with particular emphasis on the type I 
seesaw, which is well motivated by various extensions of the SM. We also comment on the supersymmetric 
seesaw framework. 

2.1 Singlet fermions (Type I) 

In the type I seesaw [36,37,39,40], two or three singlet fermions iVj (sometimes referred to as "right- 
handed neutrinos") are added the Standard Model. They are assumed to have large Majorana masses. 
The leptonic part of the Lagrangian can be written in the mass basis of the charged leptons and of the 
singlet fermions as follows: 8 

C = +[h]^(£^ c *) e ^-[A]; fc (I Q 0*)iV fc -^M J -Ar? + h.c. ) (2.3) 

where cj> is the Higgs field, with vev v u ~ 174 GeV, and the parentheses indicate anti-symmetric SU(2) 
contraction: 

(4<n = (*r Q ,er Q ) (J "o 1 ) ( t*° ) ' (2 ' 4) 

Twenty-one parameters are required to fully determine the Lagrangian of eqn (|2.3p . To see this, notice 
that it is always possible to choose a basis for the N^s where the mass matrix M is diagonal M — Dm, 
with three positive and real eigenvalues. Similarly, one can choose a basis for the i a and eu.fi such that 
the charged lepton Yukawa matrix h is diagonal and real (in particular, hh^ = D^) giving other three 
parameters. Then, in this basis, the neutrino Yukawa matrix A is a generic complex matrix, from which 

8 The Yukawa indices are ordered left-right, and the definition of A is aimed to reproduce the Lagrangian that corresponds 
to the superpotential of eqn 112.211 1. 
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three phases can be removed by phase redefinitions on the £ a , leaving 9 moduli and 6 phases as physical 
parameters. Therefore there are in total 21 real parameters for the lepton sector. See [96] for a more 
elegant counting, in particular of the phases. 

If the effective mass matrix to of the light neutrinos is normalized as in eqn (12. 2p then, in the charged 
lepton mass basis, it is given by 

[m} a p = [\} ak M~ 1 [X} 0k vl. (2.5) 

One has to take care of using a consistent set of definitions of the mass and the Higgs vev [97]. In 
particular, note the factor of 1/2 in eqn (|2.2p and our use of v u = 174 GeV, both of which are important 
in the relation (|2.5p . These choices are also important for the proof of eqn (14. 8p . so v u — 174 GeV should 
be used in the definition of rh and of eqn (|4.7p . 

The leptonic mixing matrix U is extracted by diagonalizing [m] : 

[m] = U*D m U\ (2.6) 

where D rn = diagjmi, m2, TO3}. The matrix U is 3 x 3 and unitary and therefore depends, in general, on 
six phases and three mixing angles. Three of the phases can be removed by redefining the phases of the 
charged lepton doublet fields. Then, the matrix U can be conveniently parametrized as 

U= L>-diag(l,e M ,e l/3 ) , (2.7) 

where a and /3 are termed "Majorana" phases (if the neutrinos were Dirac particles, these two phases 
could be removed by redefining the phases of the neutrino fields), and U can be parametrized in a fashion 
similar to the CKM matrix: 



U = 



C13C12 C13S12 s 13 e~ lS 

-C23S12 - S 2 3Si 3 Ci 2 e lS C23C12 - S23Sl3Sl2e 1 ' 5 S23C13 
S23S12 - C23Sl3Cl2e I<5 -S23C12 - C 23 Sl3'Sl2e 1 ' 5 C23C13 



(2- 



where = cos#y and Sij = sin 9ij. There are various other possible phase conventions for U (see e.g 
the first appendix of [98] for a list and translation rules). 

In the leptonic sector of the SM augmented with the Majorana neutrino mass matrix of eqn {23}, 
there are twelve physical parameters: The 3 charged lepton masses m e , m^, m T , the 3 neutrinos masses 
mi,m2,TO3, and the 3 angles and 3 phases of the mixing matrix U. Seven of these parameters are 
measured (m e , m^, m r , Am^, | Am§ 2 |, S12, ^23)- There is an upper bound on the mixing angle S13. The 
mass of the lightest neutrino and three phases of U are unknown. 

In addition, there are 9 unknown parameters in the high scale theory. They are, however, relevant to 
leptogenesis. 



2.1.1 Parametrizing the seesaw 

The usual "top-down parametrization" of the theory, which applies at energy scales A <; Mj, is given 
in eqn (|2.3p . To relate various parametrizations of the seesaw, it is useful to diagonalize A, which can be 
done with a bi-unitary transformation: 

A = vlD x V R (2.9) 

Thus in the top-down approach, the lepton sector can be described by the nine eigenvalues of Dm,D\ 
and Dh, and the six angles and six phases of Vl and Vr. In this parametrization, the inputs are masses 
and coupling constants of the propagating particles at energies A, so it makes "physical" sense. 

Alternatively, the (type I) seesaw Lagrangian of eqn (|2.3p can be described with inputs from the left- 
handed sector [99]. This is referred to as a "bottom- up parametrization", because the left-handed 
(5J7(2)i-doublet) particles have masses ^ the weak scale. D^, U and D m , can be taken as a subset of 
the inputs. To identify the remainder, consider the i basis where to is diagonal, so as to emphasize the 
parallel between this parametrization and the previous one (this is similar to the N basis being chosen 
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to diagonalize M). If one knows AA^ = W^ l D\Wl in the D m basis, then the N masses and mixing angles 
can be calculated: 

M" 1 - D^W L D m WlD^ = V a D]£vg (2.10) 

In this parametrization, there are three possible basis choices for the i vector space: the charged lepton 
mass eigenstate basis (-Dh), the neutrino mass eigenstate basis (-D m ), and the basis where the A is diagonal. 
The first two choices are physical, that is, U rotates between these two bases. -Dh, D m and U contain 
the 12 possibly measurable parameters of the SM seesaw. The remaining 9 parameters can be taken to 
be D\ and Vl (or W = VlU*). If supersymmetry (SUSY) is realized in nature, these parameters may 
contribute to the Renormalization Group (RG) running of the slepton mass matrix [100]. 

The Casas-Ibarra parametrization [101] is very convenient for calculations. It uses the three 
diagonal matrices Dm, D rn and -D h ht , the unitary matrix U defined in eqn (|2.6[) and a complex orthogonal 
matrix R. In the mass basis for the charged leptons and for the singlet fermions, which is used in the 
Lagrangian of eqn Q2.3[) . R is given by 

R = D m 1 / 2 U T \D M 1/2 v u . (2.11) 

R can be written as R — R diag{±l, ±1, ±1} where the ±1 are related to the CP parities of the iVj, and 
R is an orthogonal complex matrix 



R 



C13C12 C13S12 813 

-C23S12 ~ S23S13C12 C23C12 - S23S13S12 S23C13 
S23S12 - C23S13C12 -S23C12 - C23S13S12 C23C13 



(2.12) 



where cy = cosz.y, Sy = sinz-y, with complex angles. 

An alternative parametrization [102], that can be useful for quasi-degenerate N or v masses, is given 

by 

R = Ox exp{iA}, (2.13) 

where O is a real orthogonal matrix, and A a real anti-symmetric matrix. 

In summary, the lepton sector of the seesaw extension of the SM can be parametrized with Dh, the 
real eigenvalues of two more matrices, and the transformations among the bases where the matrices are 
diagonal. The matrices-to-be-diagonalized can be chosen in various ways: 

1. " top-down" - input the N sector: Dm, D x ->\, and Vr and Vl- 

2. " bottom-up" - input the vl sector: D m , Dx\t, and Vl and U. 

3. "intermediate" - the Casas-Ibarra parametrization: Dm, D m , and U and R. 



2.1.2 The two right-handed neutrino model 

The minimal seesaw models that are viable have only two A^'s [103]. Such models, known as "two right- 
handed neutrino" (2RHN) models, have a strong predictive power. The 2RHN model can be thought 
of as the limit of the three generation model, where iV 3 decouples from the theory because it is either 
very heavy or has very small couplings, that is, (\X 3a \ 2 M^ 1 )/(\X ia \ 2 M^ 1 ) — > for i = 1,2. In the 2RHN 
models, M is a 2 x 2 matrix, and A is a 3 x 2 matrix. It was originally introduced and studied with a 
particular texture [103, 104], and later studied in general (see e.g. [105] for a review). Here we follow the 
parametrization of the general model used in [106]. 

The Lagrangian of the lepton sector is of the form of eqn (|2.3|) . but with the sum over the singlet 
indices restricted to two generations (a, — e,fi,r; j,k = 1,2). The model has fourteen independent 
parameters, which can be classified as follows: 5 masses (m e , m^, m T , Mx, M%), and 9 real parameters in A 
(three phases can be removed from the six complex elements of A by phase redefinitions on the doublets). 

In this model, the lightest neutrino is massless, mx = 0, and so its associated phase vanishes, a = 0. 
This situation implies that there are ten parameters that can be determined by low energy physics, 
instead of the usual twelve. 
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[Vr] (2-14) 



In the mass eigenstate bases of charged leptons and singlet fermions, the 3x2 matrix A can be written, 
analogously to eqn (|2.9[) . as follows: 

[00 

[Vl^ A 2 
. A 3 

where [Vr] is a 2 x 2 unitary matrix with one angle and one phase, [Vl] is a 3 x 3 matrix which can be 
written as U* , with 4 of the six off-diagonal elements of W vanishing. 

As mentioned above, the main prediction of this model is that one of the doublet neutrinos is massless. 
This excludes the possibility of three quasi-degenerate light neutrinos, and leaves only two other options: 

1. Normal hierarchy, with m\ — 0, m 2 = rn so \, m 3 = m atm . 

2. Inverted hierarchy, with 777,3 = 0, m± t 2 ~ m atm , 7712 — mi = TOg ol /(2rn a t m ). 

2.2 Triplet scalars (Type II) 

One can generate neutrino masses by the tree level exchange of 5i7(2)-triplet scalars [40,107-110]. These 
S'{7(2)-triplets should be color-singlets and carry hypercharge Y = +1 (in the normalization where the 
lepton doublets have Y = — 1/2). In the minimal model, there is a single such scalar, which we denote 
by T. The relevant new terms in the Lagrangian are 

C T = -M T \T\ 2 + ~ {[X L ] aP £ a £ f3 T + M T A <^P + h.c.) . (2.15) 

Here, Mt is a real mass parameter, Az, is a symmetric 3x3 matrix of dimcnsionless, complex Yukawa 
couplings, and A^ is a dimcnsionless complex coupling. 

The exchange of scalar triplets generates an effective dimension-5 operator ££<f>(j) which, in turn, leads 
to neutrino masses. The triplet contribution to the neutrino masses, mn, is 

[mnU = (2-16) 

This model for neutrino masses has eleven parameters beyond those of the Standard Model: 8 real 
and 3 imaginary ones. Of these, 6 + 3 can in principle be determined from the light neutrino parameters, 
while 2 (Mt and lA^I) are related to the full high energy theory. 

The model involves lepton number violation because the co-existence of Al and A^ does not allow a 
consistent way of assigning a lepton charge to T, and new sources of CP violation via the phases in Al 
and X^. 

The supersymmetric triplet model [111] must include, for anomaly cancellation, two triplet superfields, 
T and T, of opposite hypercharge. Only one of these couples to the leptons. The relevant superpotential 
terms are 

W T = M T Tf + ({\ L ] al3 L a TL/3 + \ Hd H d TH d + X„ u H u TH u ) , (2.17) 

leading to 

[mnW = [AlU^I (2.18) 

2.3 Triplet fermions (Type III) 

One can generate neutrino masses by the tree level exchange of S'C/(2)-triplet fermions T" [112-114] (i 
denotes a heavy mass eigenstate while a is an SU(2) index). These Si7(2)-triplets should be color-singlets 
and carry hypercharge 0. The relevant Lagrangian terms have a form that is similar to the singlet-fermion 
case (12.31) . but the contractions of the SU(2) indices is different, so we here show it explicitly: 

£ Ta = [\ T ] ak r a p JP a m - \MiT?T? + h.c. (2.19) 
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Here, M, are real mass parameters, while At is a 3 x 3 matrix of dimensionless, complex Yukawa couplings. 

The exchange of fermion triplets generates an effective dimension-5 operator ££<jxj) which, in turn, 
leads to neutrino masses. The triplet contribution to the neutrino masses, mm, is 

[mm\ aP = [XrUkM^Xpkvl. (2.20) 

As in the standard seesaw model, this model for neutrino masses has eighteen parameters beyond 
those of the Standard Model: 12 real and 6 imaginary ones. 

The model involves lepton number violation because the co-existence of At and Mk does not allow a 
consistent way of assigning a lepton charge to Tf., and new sources of CP violation via the phases in At- 



2.4 Supersymmetry and singlet fermions 

One of the motivations of supersymmetry (SUSY) is to cancel quadratically divergent contributions to 
the Higgs mass. In the seesaw extension of the SM, a large mass scale Mj is present. This results in 
an additional set of corrections of 0(A 2 M 2 ). 9 Cancelling the contributions from the large seesaw scale 
motivates the supersymmetric version of the model. 

The superpotential for the leptonic sector of the type I seesaw is 

W = \N?MiN? + {L a H u )[\] ai N? - (L a H d )[h] a E c a , (2.21) 

where L a and are respectively the SU{2) doublets and singlets lepton superfields, and H u and H4 are 
the Higgs superfields. (To resemble the SM notation, the scalar components of the Higgs superfields will 
be denoted as 4>u and <j>d-) The SU(2) contractions in parentheses are anti-symmetric, as in eqn (|2.4ll . 
Both Higgs bosons, <fi u and tfid, have vacuum expectation values (vevs): (4>i) = Vi. Their ratio is defined 
as 

tan /3 ee— . (2.22) 

Vd 

When including "flavour effects" in supersymmetric leptogenesis, the value of tan/3 is relevant, because 
the Yukawa coupling h e Ai r oc m e Ai r / cos 0. 

To agree with experimental constraints (while keeping supersymmetry as a solution to the m^-fine- 
tuning problem), it is important to add soft SUSY-breaking terms: 

\[fhl] a0 L* a L p + ... + ^N?[BM]ijNj + [AX] m (L a (b u )N^ + [Ah] 0a (L a <p d )E° + h.c, (2.23) 

where the ... represent soft masses-squared for all the scalars. In the thermal leptogenesis scenario that 
this review concentrates on, the soft SUSY breaking terms give 0(m| ugY /M 2 ) corrections, and can be 
neglected. In other mechanisms, such as soft leptogenesis [86-88] and Affleck-Dine leptogenesis [48,49,89], 
the soft parameters play a central role. 

The interesting feature of the SUSY seesaw, is that the neutrino Yukawa couplings may contribute 
to the RG running of the soft slepton mass matrix, and induce flavour-changing mass-squared terms. 
Consider, for example, the Type I seesaw, with universal soft masses at some high scale A: [m 2 ] Q /3 = 
m o^a/?5 [^.A] a i = Ao[A] a j. Then at the electroweak scale, in the flavour basis, the RG contributions to 
the off-diagonal elements of the mass-squared matrix can be estimated at leading log as follows: 

[mlU = -^f^ E A ™ l0 S ^t X h (2-24) 

i 

In general, the soft mass matrix is unknown. One can argue, however, that the flavour-changing mass- 
squared matrix elements (off-diagonal in the flavour basis) are at least of order eqn (|2.24p , because we do 

9 The requirement of no excessive fine tuning in the cancellation of these contributions has been used to set the bound 
Mi <> 10 7 GeV [115]. 
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not expect fine-tuned cancellations among different contributions. The flavour off-diagonal terms induce 
processes like fj, — > ej, t — > /Z7 and r — ► 67 [100], at rates of order [116] 

r^vr) _« 3 IKWI 2 . (1 + tan2/3)j (2 . 25) 



T{t a - ^w) Gf 



F m susy 

where m susy is a typical slepton mass. For reasonable values of the unknown parameters, one obtains 
predictions [101,116-118] that are in the range of current or near-future experiments [119-121]. 

The neutrino Yukawa couplings make real and imaginary contributions to the soft parameters. These 
phases give contributions to lepton electric dipole moments [122-124] which are orders of magnitude 
below current bounds, but possibly accessible to future experiments. 

The seesaw contribution to RG running of the slepton masses is relevant to leptogenesis, because the 
slepton mass matrix could be an additional low energy footprint of the seesaw model. As discussed in 
Section [2~T1 it is possible to reconstruct the Yukawa matrix A and the masses Mj, from the light neutrino 
mass matrix m u of eqn ()2.5|) . and the Yukawa combination AA^ that enters the RG equations. So it is 
interesting to study correlations between low energy observables and a large enough baryon asymmetry 
at low enough T re h e at- Early (unflavoured) studies can be found in refs. [125,126], and in many other 
works: [127] (degenerate^), [118] (hierarchical Ni), [102] (degenerate light neutrinos), and [128] (type II). 
Recent flavoured analyses can be found in refs. [129,130] (hierarchical Ni) and [131,132] (degenerate N). 
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3 Anomalous B + L Violation 



The aim of this section is to give a qualitative introduction 10 to the non-perturbative baryon number 
violating interactions that play a crucial role in leptogenesis. A similar discussion can be found in 
Ref. [67], while more details and references can be found, for instance, in Section 2 of Ref. [17], and in 
Refs. [138-141]. 

From a theoretical perspective, the baryon number B and the three lepton flavour numbers L a arc 
conserved in the renormalisable Lagrangian of the Standard Model. Furthermore, experimentally, the 
proton has not been observed to decay: r p 10 33 years [142, 143] . (For a review of proton decay, see [144] .) 
However, due to the chiral anomaly, there are non-perturbative gauge field configurations [11,145,146] 
which can act as sources for B + L e + + L T . (Note that B — L e — — L T is conserved.) In the early 
Universe, at temperatures above the electroweak phase transition (EWPT), such configurations occur 
frequently [10,22,147], and lead to rapid B + L violation. These configuations are commonly referred to 
as "sphalerons" [148-150]. 

3.1 The chiral anomaly 

For a pedagogical introduction to the chiral anomaly, see for instance Ref. [138]. 

Consider the Lagrangian for a massless Dirac fermion ip with U(l) gauge interactions: 



C = tlrfpn - iA^ - ^F^F^. (3.1) 



It is invariant under the local symmetry: 

^x) _> e ^ x U(x) , A^x) - A^x) + d^x), (3.2) 
It is also invariant under a global "chiral" symmetry: 

i/)(x) -> e^^x). (3.3) 

The associated current, 

ft = ^757^, (3.4) 

is conserved at tree level, but not in the quantum theory. This can be related to the regularization of 
loops — renormalization introduces a scale, and the scale breaks the chiral symmetry, as would a fermion 
mass (see, for instance, chapter 13 of Ref. [138]). Indeed, at one loop, one finds 

« = Ttr^ F ^ = ^F Pa F" V ■ (3-5) 

The right-hand side can be written as a total divergence involving gauge fields, and is related to their 
topology: it counts the "winding number" , or Chcrn-Simons number, of the field configuration. (An 
instructive 1+1 dimensional model, where the topology is easy to visualize, can be found in Ref. [141].) 
In four dimensions, the space-time integral of the right-hand side of eqn (|3.5|) vanishes for an Abelian 
gauge field, but can be non-zero for non-Abelian fields. 

In the context of leptogenesis, we are looking for an anomaly in the B + L current. Within the 
four-dimensional SM, it arises due to the SU(2) gauge interactions, which are chiral and non-Abelian. 
We neglect other interactions in the following (see Ref. [17] for a discussion of the effects of Yukawa and 
SU(3)c x U(1)y interactions). The fermions that are relevant to our discussion are the three generations 
of quark and lepton doublets: {^\} = {Ql^i where a, (3 are generation indices, a, b are colour indices, 
and A, B are SU(2) indices. The Lagrangian terms for the SU{2) gauge interactions read 

c = £^zy (a„ - i 9 - ^Ovi- (3.6) 



10 This Section is based on a lecture given by V. Rubakov at the Lake Louise Winter Institute, 2008. 
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It has twelve global £7(1) symmetries (one for each field): 



tfiOr) - e^fytol . (3.7) 
The chiral currents associated to these transformations, 

#=^1*7^1. (3-8) 
are conserved at tree level, but are "anomalous" in the quantum theory: 

d " j j = eh F ^ vA (3 ' 9) 

Let us define Q l (t) — J j Q l d?x, AQ Z — Q' l (+oo) — Q l (—oo), and let us suppose for the moment that 
there exist field configurations such that 



is a non-zero integer. This implies that fermions will be created, even though there is no perturbative 
interaction in the Lagrangian that generates them. One way [151] to understand "where they come from" 
is to think in the Dirac sea picture, and place the chiral fermions {V'i} m & n external gauge field for 
which the right hand side of eqn (|3.10|) is non-zero. In the ground state at t — > — oo, all the negative 
energy states are filled, and all the positive energy states are empty. As the fermions are massless, there 
is no mass gap at E = 0. At any given t, one can solve for the eigenvalues of the fermion Hamiltonian. 
See, for instance, Ref. [139] for a discussion. One finds that the levels move as a function of t: negative 
energy states from the sea acquire positive energy, and empty positive energy states could become empty 
sea states. In the case of the chiral SU(2) of the SM, one finds that, for each species of doublet, what 
was a filled left-handed state in the sea at t — * — oo, becomes a particle at t — * +oo. See figure |3~T1 
This "level-crossing" occurs for each type of fermion, so the gauge field configuration centered at t = 
in figure [3~T1 is a source for nine quarks and three leptons. 

Left-handed fermions 

h 




Figure 3.1: Evolution with time of the energy eigenstates of chiral fermions in a gauge field background 
with FF ^ 0. 



3.2 B + L violating rates 

At zero temperature, gauge field configurations that give non-zero J d 4 xFF correspond to tunneling 
configurations, and are called instantons [152] (for reviews, see e.gRefs. [140,141]). They change fermion 
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number by an integer N, so the instanton action is large: 



/ d* X F*F^ >\-Lj d*xF*F^ 



64ir 2 N 



The first inequality follows from the Schwartz inequality (see [141]). Consequently, the associated rate is 
highly suppressed, 

Y ^ e — (instanton Action) ^ ^— Atv / ct\y 

and the mediated B + L violation is unobservably small. Moreover, the instantons do not threaten 
the stability of the proton [11], because an instanton acts as a source for three leptons (one from each 
generation), and nine quarks (all colours and generations), so it induces AB = AL = 3 processes. Notice 
that the three quantum numbers B/3 — L a are not anomalous, so they are conserved in the SM. 

If the ground state of the gauge fields is pictured as a periodic potential, with minima labeled by 
integers, then the instantons correspond to vacuum fluctuations that tunnel between minima. With this 
analogy, one can imagine that at finite temperature, a thermal fluctuation of the field could climb over the 
barrier. The sphaleron [148-150] is such a configuration, in the presence of the Higgs vacuum expectation 
value. The B + L violating rate mediated by sphalerons is Boltzmann suppressed: 

r sph cx e - E ^ T , 

where E sp h = 2Bmw / ctw is the height of the barrier at T = 0, and 1.5 ^ B ^ 2.75 is a parameter that 
depends on the Higgs mass. 

For leptogenesis, we are interested in the B + L violating rate at temperatures far above the EWPT. 
The large B + L violating gauge field configurations occur frequently at T » mw [153-157]. The rate 
can be estimated as (see [158] for a recent discussion) 

~ 250 a^T . (3.11) 

This implies that, for temperatures below 10 12 GeV and above the EWPT, B + L violating rates are in 
equilibrium [158,159]. 
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4 A Toy Model of Thermal Leptogenesis 



In this section, the baryon asymmetry produced by thermal leptogenesis is estimated. The goal is to 
provide a basic understanding and useful formulae, while avoiding, at this stage, the details. The CP 
asymmetry in the singlet fermion decay is introduced as a parameter, and discussed in more detail in 
Section O The dynamics is estimated by dimensional analysis; Boltzmann equations appear in Section [6] 
In this section, we make the following simplifying assumptions: 

1. The lepton asymmetry is produced in a single flavour a; 

2. The A-masses are hierarchical: Mi ~ 10 9 GeV <§; M 2 ,M 3 . (The kinematics is simpler in an 
effective theory of a propagating Ni and effective dimension- five operators induced by N 2 and N 3 ). 

3. Thermal production of Ni and negligible production of A3. 

When any of these three assumptions does not hold, there are interesting modifications of the simplest 
scenario that we describe in this section. The effects of flavour arc discussed in section[9j the consequences 
of non-hierarchical Mj's (and the possible effects of ^2,3) are summarised in section 110.21 and other 
leptogenesis mechanisms were mentioned in section [1.31 

The basic idea is the following. Scattering processes produce a population of N±s at temperatures 
T ~ Mi. Then this N± population decays away because, when the temperature drops below Mi, the 
equilibrium number density is exponentially suppressed oc e~ Ml ^ T . If the N± interactions are CP violating, 
asymmetries in all the lepton flavours can be produced. If the relevant interactions are out-of-equilibrium, 
the asymmetries may survive. They can then be reprocessed into a baryon asymmetry by the SM B + L 
violating processes that have been discussed in Section [3] 

A unique feature of thermal leptogenesis, which distinguishes it from other "out-of-equilibrium decay" 
scenarios of baryogenesis, is that the same coupling constant controls the production and later disap- 
pearance of the population of iV's. As demonstrated in the seminal works [53, 133], a sufficient number 
density of N's can be produced via their Yukawa coupling A. The CP asymmetry in the processes that 
produce the N population is closely related to the CP asymmetry in the N decays, which wipe out the 
N population. In particular, in our toy model of hierarchical A^'s, the CP asymmetry in the scattering 
interactions by which the Ni population is produced is equal in magnitude but opposite in sign to the 
CP asymmetry in Ni decays. At first sight, this suggests that the final lepton asymmetry is zero. 11 A 
non-zero asymmetry survives, however, because the initial anti-asymmctry made with the N population 
is depleted by scattering, decays, and inverse decays. This depletion is called washout, and is critical 
to thermal leptogenesis. The importance of flavour effects in leptogenesis is a consequence of the crucial 
role played by washout: the initial state of washout interactions contains a lepton, so it is important to 
know which leptons are distinguishable. 

Our aim here is to estimate the asymmetry-to-entropy ratio by considering the Sakharov conditions. 
Each condition gives a suppression factor. The baryon asymmetry can be approximated as 



The first factor is the equilibrium Ni number density divided by the entropy density at T > Mi, of 
0(4 x 10~ 3 ) when the number of relativistic degrees of freedom is taken ~ 106, as in the SM. An 
equilibrium number density of Ai's is the maximum that can arise via thermal Ni production. It is 
produced when A Q i is large. A smaller Ni density is parameterized in r/ a . As concerns the other factors 
in (|4.1j) . we note the following: 

11 Notice that the potential cancellation is between the CP asymmetry in processes with N and £ a in the final state, 
such as X — > N£ a scattering, and the asymmetry in processes with N in the initial state and l a in the final state, such 
as TV — > tt>£ a - Only processes with i a in the final state can generate an asymmetry. In particular, there is no cancellation 
between the asymmetry produced in decays and inverse decays. It is intuitive, and straightforward to verify (see Section 
[BJ, that interactions with the lepton in the initial state can wash-out the asymmetry, but not produce it. 




(4.1) 
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1. e aa is the CP asymmetry in N± decay. For every l/e aa N± decays, there is one more £ a than there 
are s. 

2. rj a is the efhciency factor. Inverse decays, other "washout" processes, and inefficiency in Ni pro- 
duction, reduce the asymmetry by < r\ a < 1. In particular, rj a — is the limit of Ni interactions 
in perfect equilibrium, so no asymmetry is created. 

3. C describes further reduction of the asymmetry due to fast processes which redistribute the asym- 
metry that was produced in lepton doublets among other particle species. These include gauge, 
third generation Yukawa, and B + L violating non-perturbative processes. As we discuss in sections 
[8] and 19. 3[ C is a matrix in flavour space, but for simplicity we approximate it here as a single 
number. 

Formulae for e aa and r\ a can be located in this review, by consulting the table in section 1121 Our aim is 
now to estimate e Qa , rj a and C . 

4.1 CP violation (e aa ) 

To produce a net asymmetry in lepton flavour a, the N\ must have L Q -violating interactions (see section 
14. 3|) . and different decay rates to final states with particles or anti-particles. The asymmetry in lepton 
flavour a produced in the decay of Ni is defined by 

£aa ~ r(JVi — > 4>g)+r(N 1 — ► <t>i) [ ' 

where I denotes the anti-particle of t. N\ is a Majorana fermion, so N% = N% 12 . The asymmetry e aa is 
normalized to the total decay rate, so that the Boltzmann Equations are linear in flavour space. When we 
include additional lepton generations, we find that the CP asymmetry is a diagonal element of a matrix, 
so we give it a double flavour index already. 

By definition, |e QQ | < 1. Usually, it is much smaller than that. It is a function of the parameters 
of the Lagrangian (|2.3|) . This dependence is evaluated in section [5] The requirement that it is large 
enough to account for the observed baryon asymmetry (roughly speaking, |e aa | > 10 -7 ) gives interesting 
constraints on these parameters. 

4.2 Out-of-equilibrium dynamics (i] a ) 

The non-equilibrium which is necessary for thermal leptogenesis is provided by the expansion of the 
Universe: interaction rates which are of order, or slower, than the Hubble expansion rate H are not 
fast enough to equilibrate particle distributions. Interactions can be classified as much faster than H, 
much slower, or of the same order. For the purposes of making analytic estimates (and writing Boltzmann 
codes), it is convenient to have a single scale problem. The timescale of leptogenesis is -ff -1 , so we neglect 
interactions that are much slower than H . Interactions whose rates are faster than H are resummed into 
thermal masses, and impose chemical and kinetic equilibrium conditions on the distributions of particles 
whose interactions are fast. 

For the initial conditions, we assume that after inflation the Universe reheats to a thermal bath at 
temperature T re h oa t which is composed of particles with gauge interactions. A thermal number density of 
N% (n^ ~ n.y) is produced if T re h e at ~ Mi/ 5 [59,60], and if the production timescale for iVi's, l/r pro d, 
is shorter than the age of the Universe ~ 1/H. The N\ can be produced by inverse decays cf)£ a — * Ni 
and, most effectively, by 2 — > 2 scatterings involving the top quark or electroweak gauge bosons. We 
here neglect the gauge interactions (although > h t at 10 10 GeV) because it is simpler and formally 
consistent, and because the O(a) corrections do not give important new effects. N\ can be produced by 

12 In supersymmetry, Ni represents the chiral superfield, so one may wish to distinguish Ni from JVi. In that case, the 
e aa arise in corrections to a "D-term", and can be denned by replacing r(7Vi — > <j>£ a ) — * r(JVi — » <j>£ a ) m eqn l4.2l 
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s or t channel exchange of a Higgs : qhtn. — > 4> ~ * ?-aN or £ a tn — > — > ^lA^. So the production rate can 
be estimated (by dimensional analysis in zero temperature field theory) as 

r prod ~ £ . (4.3) 

a 

If r pro d > i? then, since h t ~ 1, the TVi total decay is also "in equilibrium": 

r ZJ >fr(r = Mi), (4.4) 

where 

r g = ]T r aa = J2 r(Nj. -> <R,, <R») = [At A ^ lMl , (4.5) 



and 

T 2 

ff(T = Mi) = I.663* 72 

m pl 



(4.6) 

T=Mi 

Here $* is the number of relativistic degrees of freedom in the thermal bath (see eqn 113. 8p . Within the 
SM, 5 * = 106.75, 

It is useful to introduce two dimensionful parameters [133] m and ro*, which are of the order of the light 
neutrino masses, and which represent, respectively, the decay rate To and expansion rate H(T — Mi): 

m «« = 2^ Ml = 8it m2 Td > 

a 1 

to* = 8tt-^% i?| T=Ml ~ 1.1 x 1(T 3 eV. (4.7) 

It can be shown [134] that 

fh > TOmin, (4.8) 

where m m ; n is the smallest light neutrino mass (this is relevant for degenerate light neutrino masses), 
and that "usually" m <; m so i [135] (see also [136]). The > H condition reads, in the language of fh 
and to* , simply as 

to > m* . (4-9) 

If indeed fh > m so i, then this condition is satisfied. This range of parameters is referred to as "strong 
washout" . The converse case, to < to* , is referred to as "weak washout" . 

In the strong washout scenario, at T ~ M 1; a thermal number density of N\ is obtained (njvj ~ n»), 
and the total lepton asymmetry Yj, ~ (any asymmetry made with the ATj is washed out). As the 
temperature drops and the iVi's start to decay, the inverse decays — * Ni, which can wash out the 
asymmetry, may initially be fast compared to H. Suppose that this is indeed the case for flavour a. Then 
the asymmetry in lepton flavour a will survive once the partial inverse decays from flavour a are "out of 
equilibrium" : 

T ID {4>t a -> Nx) \T aa e- Ml ' T <H = 1.66^7— (4.10) 

2 TOpi 

where the partial decay rate r QQ is defined in eqn (|4.5[) . and Tjd ~ e~ Ml / T T At temperature T a 
where eqn (|4.10p is satisfied, the remaining Ni density is Boltzmann suppressed, oc e~ Ml ^ Ta . Below T a , 
the Ni 's decay " out of equilibrium" , and contribute to the lepton flavour asymmetry. So the efficiency 
factor rj a for flavour a can be estimated as 

n Nl (T a ) -Mi/Tc _ m * r - ^ \ rAii\ 

*!<*— trr ^ 1\ T \ - ~ (?71 > TO*, TO QCe > TO*), (4.11) 

njvAi > Mi) TO a „ 
where m*/fh aa = H(T = Mi)/T aa . This approximation applies for to qq > to* ~ 10 _3 eV. 
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Now consider an intermediate case, where to > to* (strong washout), but rh aa < to*. In this case, 
the N\ number density reaches its equilibrium value, because it has large couplings to other flavours, but 
the coupling A Q i to the flavour we are interested in is small. The (anti-)asymmetry in flavour a is of 
order — e QQ n 7 . As the population of Ni decays at T ^ Mi, a lepton asymmetry ~ e QQ n 7 is produced. 
Consequently, at lowest order in rh aa , the lepton asymmetry vanishes. 13 A small part of the asymmetry 
(in flavour a) made during N\ production is, however, washed out before Ni decay. This fraction can be 
estimated as ~ — (TO aQ ,/m*)e QQ n 7 , yielding an efficiency factor 1 

rj a ~ — — (to > to*, rh aa < to*). (4-12) 

TO* 

In the weak washout scenario, not only rh aa < to*, as above, but also the total decay rate is 
small, to < to*. In this case, the Ni number density does not reach the equilibrium number density 
~ 7i 7 . Production is most efficient at T <~ Mi, when the age of the Universe satisfies 2tjj = 1/H, so 
njvj ~ r pro( jr[/n 7 ~ (rn/m*)n 7 . The cancellation, at lowest order, of the lepton anti-asymmetry and 
asymmetry is as in the intermediate case above, so that the efficiency factor can be estimated as 

i] a ~ — ^ — (rh < m*, m QQ < m*). (4.13) 



4.3 Lepton and B + L violation (C) 

The interactions of Ni violate L because lepton number cannot be consistently assigned to iVi in the 
presence of A and M. If L(N\) = 1, then A Q i respects L but Mi violates it by two units. If L(Ni) = 0, 
then Mi respects L but A Q i violates it by one unit. The Ni decay, which depends on both Mi and A Q i, 
does not conserve L. The heavy mass eigenstate is its own anti-particle, so it can decay to both £(f> and 
£(f>* . If there is an asymmetry in the rates, a net lepton asymmetry will be produced. 

The baryon number violation is provided by B + L changing SM non-perturbative processes [17] (see 
Section [3]) . Their approximate rate is given in eq. (|3.11|) and is faster than the Hubble expansion H in 
the thermal plasma for T 10 12 GeV. The asymmetry in lepton flavour a, produced in the N\ decay, 
contributes to the density of B/3 — L a , which in conserved by the SM interactions. In equilibrium, this 
excess of B — L implies (for the SM) a baryon excess [41]: 

a 

where Ya q is the asymmetry in B/3 — L a , divided by the entropy density. The value of C = 12/37 applies 
in the SM (see ean ll4.17j) . In the MSSM, it is 10/31 (see eqn ll4.18p . 

So far, the focus has been on the neutrino Yukawa interactions, which produce an asymmetry in the 
lepton doublets £ a . The SM interactions, which redistribute the asymmetries among other particles, are 
included in Section [H As discussed in section [973l these interactions usually give O(l) effects, which are 
parameterized with the A-matrix [65] that is derived in section 1141 One effect that can be explained 
already at this stage is the following. When the charged lepton Yukawa coupling h a is in chemical 
equilibrium, that is, when interaction rates such as T(gauge boson-fe^ «-» t + (j>) are fast compared to H, 
the lepton asymmetry in flavour a is shared between eRa and £ a . But only the part that remains in £ a 
is washed out by the neutrino Yukawa interactions, so there is a mild reduction in washout. 

13 This cancellation is discussed in more detail in ref. [59], when they consider production by inverse decays. This 
cancellation is absent in their later discussion when production by scattering is included, because CP violation in scattering 
was neglected. 

14 This estimate assumes the momentum distribution f(p) is thermal (see Appendix [13} . This is a usual assumption in 
leptogenesis, where Boltzmann Equations for the total number densities are solved. Differences that arise when the BE are 
solved mode by mode were studied in [137]. 
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4.4 Putting it all together 

An estimate for the baryon asymmetry can be obtained from eqn (|4.ip . with the prefactor C from eqn 
631: 

Y AB ~ 10~ 3 e aa r] a , 

where the CP asymmetry e aa is taken from eqn (|5.9[) or (|5.13p , and the efficiency factor is taken from 
eqn (|4.11[) . (|4.12[) or (|4.13[) . To obtain more accurate estimates, as can be found in Appendix [T6l the 
dynamics should be calculated via the Boltzmann Equations, introduced in Section [5] 

We can anticipate the flavour issues, which are discussed in Section [9j by supposing that there are 
CP asymmetries in all flavours. Then, in the strong washout case for all flavours, we obtain 

Y AB ~ 

i] a ~ 10 3 ^ , aa (flavoured, strong washout) (4-15) 

Tn act 

a a 

where the flavours summed over are presumably the charged lepton mass eigenstates. Alternatively, one 
might choose l a — 1^, the direction in flavour space into which Ni decays (see eqn I5.11| ). Then en is 
the total CP asymmetry e and T(N — > <j>£i) is the total decay rate Td- One obtains 

Yab ~ 10 _3 m* 4- (single flavour, strong washout), (4.16) 
m 

which is simpler but different. In Appendix ll5.4l (see also Section[9]) we discuss why and when the charged 
lepton mass basis is the relevant one. 
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5 CP Violation 



In section 15.11 we review how to calculate a CP asymmetry in N± decays, including the vertex and 
wavefunction [160] contributions. In Section 15.21 wc calculate e aa for hierarchical N, and give the 
formulae for less hierarchical N, as calculated in [56]. It was shown in [73,74], that for hierarchical Ni, 
there is an upper bound on the CP asymmetry proportional to M\. Section [5.21 contains a derivation of 
this bound, which gives a strong constraint on models, and a discussion of various loopholes. Sections 
I5.3H5.4I discuss general constraints from S-matrix unitarity and CPT invariance, which have implications 
for the generation of a cosmological asymmetry from decays. In particular, we explain a subtlety in the 
analysis: the contribution of an on-shell N\ to scattering rates should be subtracted, as it is already 
included in the decays and inverse decays. 

5.1 CP violation in N\ decays 

The CP asymmetry in lepton flavour a, produced in the decay of N\, is defined in eqn (14. 2j) : 

r(JVi->^ a )-r(iVi-^^ a ) 



,-<t> 




Figure 5.1: The diagrams contributing to the CP asymmetry e aa . The flavour of the internal lepton 
£p is summed. The internal tp and Higgs 4> are on-shell. The X represents a Majorana mass insertion. 
Line direction is "left-handedness" , assigning to scalars the handedness of their SUSY partners. The loop 
diagrams on the first line are lepton flavour and lepton number violating. The last diagram is lepton 
flavour changing but "lepton number conserving" , in the sense that it makes no contribution to the total 
CP asymmetry e. It is suppressed by an additional factor Mi/M 2 ,3 [see eqn (15.131) ]. 

The CP asymmetry e aa arises from the interference of tree-level (subscript 0) and one-loop (subscript 
1) amplitudes. This is discussed further in section 13751 As noted in [160], it is important to include all 
the one loop diagrams, including the wavefunction corrections. The tree and loop matrix elements can 
each be separated into a coupling constant part c and an amplitude part A: 

M = M + Mi = cqAo + ciA\ . (5.2) 

For instance, in the tree level decay of figure (|5.1|) . 

c = A* 1 Aq(N — > (j>H a ) = u ta P R u N . (5.3) 

The matrix element for the CP conjugate process is 

M = c^Ao + c\Ax, (5.4) 
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where 15 \Ai\ 2 — \Ai\ 2 . Thus the CP asymmetry can be written 

/ \cqAq + ciAi\ 2 5 fflt,* - J \c* Aq + c\Ai\ 2 5 dU e , 



2 E/3 / Ico^ol 8 dU, 
Imfac?} 2 Jlm{AoAl}5 dU e ^ 



Eakol 2 J\A \ 2 SdU e 



(5.5) 



where 



and Pi, Pf are, respectively, the incoming four-momentum (in this case P)v) and the outgoing four- 
momentum (in this case +pi). The loop amplitude has an imaginary part when there are branch cuts 
corresponding to intermediate on-shell particles (see Cutkosky Rules in [161], or eqn (|5.19|) ). which can 
arise in the loops of figure [57X1 when the <j> and 1$ are on-shell: 



2Im{A^i} = A Q (N -> 04) [^* (N i0<t>') 6' dn^AHl'^' <P£ a ). (5.7) 

J 

Here 0' and fjj are the (assumed massless) intermediate on-shell particles, and dRy p is the integration 
over their phase space. 

5.2 e aa and the lower bound on Mi 

In the limit M2, M3 ^ Mi, the effects of N2, N3 can be represented by an effective dimension-5 operator. 
In the diagram of 15. 1[ this corresponds to shrinking the heavy propagator in the loop to a point. For 
calculating e aa , the Feynman rule for the dimension-5 operator can be taken oc [m]/v 2 . (There is a 
contribution to [m] from N\ exchange, which is not present in the dimension-5 operator that is obtained 
by integrating out N2 and N3. But the iVi-mediated part of [m] makes no contribution to the imaginary 
part for e aa .) Then, we obtain for the relevant coupling constants 

c = A* 1 ci = 3 A/31 [m*]f3 a /v 2 . 



The factor of three comes from careful book-keeping of weak SU(2)l indices; The dimension-5 operator 
is 

^K0o - el^)(v% - e^ + ) + h.c. (5.8) 

This leads to a Feynman rule 2(d"5^ + 5^S")[m} a p / (2v 2 ) for the vertex i/£(/> o i/£0o or e i,0 +e L0 + ; but 
to a Feynman rule — [m\ pa /v 2 for f^0 o eJ0 + or ^£0 o e^0 + . Summing over all possible lepton/Higgs 
combinations in the loop gives the factor of three. It can also be seen in the theory with propagating 
iV^: The charged and the neutral components of the intermediate cj)' and i'^ contribute in the Ni wave- 
function correction, giving a factor of 2, but only the charged or the neutral <p' and l'^ appear in the 
vertex correction. 

To obtain the amplitude ratio, in the case of hierarchical TV's, we take Aq(1'q4>' — > 4>£ a ) — v^Plu^, 
and after spin sums we obtain: 16 

" ' Im {[A] al [m*A] Ql } (5.9) 



Wirvl [At A] 11 



15 In the CP conjugate amplitude, A, the ua spinors are replaced by vi spinors. Since, however, U(Ui = = Vgvg, the 
magnitude 1 2 is the same. 

16 The various 2's for the initial and final state averages and sums are discussed around eqn H13.18I I. They cancel in the 
ratio and can be ignored here. 
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where there is no sum on a in this equation. 
The upper bound 



3MiTO n 



1 '''max / u Ni | nJVi 



—a/BZ 1 +B- 1 - (5.10) 

where B££ ee r(iVi -► xy)/^ and 

m m ax is the largest light neutrino mass, can be derived by defining 

the unit vector 

i Nl = - Xal (5.ii) 

and using |m • < |m max ^Vi|- 

The upper bound on |e aa | can be used to obtain a lower bound on M\ and the reheat temperature of 
the Universe for thermal leptogenesis with hierarchical iVj. The first calculations [73,74] used the total 
asymmetry 17 e — ^ Q e aa ; the differences are discussed around eqn (15.17|) . The estimate for the baryon 
asymmetry, eqn (|4.1|) . combined with eqn (14.14|) . will match the observed asymmetry of eqn (|1.2p for 
J2 a e aaila ~ 10~ 7 . The efficiency factor, < r/ a < 1, is usually 0.1 (see section[6]), which implies 

e aa ;>ht 6 . 

Taking m max = m atm in the upper bound of eqn (|5.10p . we find 

Mi £ 10 9 GeV. (5.12) 

A more precise bound can be obtained numerically. The CP asymmetries e aa can be smaller when r\ a 
is larger, that is, when there is less washout, which occurs when Tu ~ H(T = Mi). For this range of 
parameters, analytic approximations are not so reliable (see the plots of section 15.4. 3[) . 

The bound of eqn (|5.12[) is restrictive, so it is worth repeating the list of assumptions that lead to it 
and to enumerate some mechanisms that evade it: 

1. The bound applies for non-degenerate heavy neutrinos. The CP asymmetry can be much larger for 
quasi-degenerate N i: with Mi — M 2 ~ r_o (see section [T0.2.2p . 

2. The bound actually applies only for strongly hierarchical heavy neutrinos. To prove the upper bound 
on e, for arbitrary A compatible with the observed [m], requires M 2j z > lOOMi [162,163]. For a 
milder hierarchy, it is possible to tune A such that a dimension-seven operator gives a contribution to 
e that can be as high as Mj 5 j '(MfMa) and thus possibly exceed the bound [164], but no contribution 
to [m]. If such tuning is neglected, the bound is "usually" good for M 2 ^ > lOMi. 18 

3. The bound can be evaded by adding particles and interactions [165]. Some of the possibilities are 
the following: 

• In multi-Higgs models, the vev(s) of the scalars that appear in the light neutrino mass matrix 
may be unknown. The CP asymmetries e aa increase as these vevs become smaller. 
For instance, writing v u = 174 sin/3 GeV in a model with two Higgs doublets, the bound 
becomes Mi <; sin 2 /3 x 10 9 GeV, which can be significantly weaker if sin/3 <C 1 [166]. It is 
interesting to note, in this regard, that this is not the case in the supersymmetric Standard 
Model, 19 in spite of its being a two-Higgs model. The reason is that in this framework tan /3 > 
1. If there were extra Higgs doublets, or a non-analytic term LH^N with large tan/3, then 
thermal leptogenesis could be successful for lower Mi and T re heat values. 



17 Notice that [74] used e for the MSSM, which is twice as big, see eqn 1 15. 151 . The bound given in [74] therefore had 87r 
rather than Wit in the denominator. 

18 This can be seen from [163], where the procedure of scanning of the parameter space is described. There is no similar 
information in [162], where contrary claims are made. 

19 It is desirable, within the supersymmetric framework, to avoid the bound 1 15. 121 1 because it may be in conflict with 
upper bounds on T re h ea t from the gravitino problem. 
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In "inverse seesaw" models [167], which contain additional singlet fermions and scalars, the 
light neutrino masses are proportional to the unknown vev(s) of the additional scalar(s). For 
sufficiently small vev(s), hierarchical leptogenesis can work down to the TeV scale [168-170]. 

• A minimal extension that works down to the TeV, is to add a singlet [171]. 

• Including a fourth lepton generation allows thermal leptogenesis at T ~ TeV [172-174]. 

• Consider the case that the number of heavy singlet neutrinos is larger than three. Their 
contribution to the effective operator (£(f>)(£(f>) has no effect on the bounds on the flavoured 
asymmetries e aa . In contrast, the contribution of many singlets in weak washout can enhance 
the final baryon asymmetry, allowing thermal leptogenesis at T rc hcat that is a factor ~ 30 
lower than in the case of three singlet neutrinos [175]. In the case where leptogenesis is 



unflavoured, extra singlets weaken the upper bound on 
3Mim 3 /(167re 2 ) [175,176]. 



|e| from 3Mi(m3 — m\)/ (\Qirv 2 ) to 



4. For degenerate light neutrinos, m max > m atm , so the individual flavour asymmetries can be larger. 
However, the efficiency factor is smaller, so for m max ^ eV (a conservative interpretation of the 
cosmological bound [177-179]) the lower bound on Mi is similar to eqn (|5.12[) [62,180,181]. 

One can go beyond the effective theory and incorporate the states as dynamical degrees of 



freedom. For a not-too degenerate Ni spectrum, Mi 
in Section [TO. 2. 2j) . one obtains 



Mj ^> Td (the case of Mi — Mj ~ To is discussed 



1 



1 



(8tt) [AtA] n 
1 



]Tlm {(A* 1 )(A^A)i J Aa i j}g (a 



(8tt) [At A] 



(5.13) 



where 

and, within the SM [56], 

g(x) = • 



Xj = Mf/Ml 



I , , / 1 + x 

+ 1 - (1 + x) In 



1-x 



2^i 6x 3 / 2 



(5.14) 



In the MSSM, N% decays to a slepton + Higgsino, as well as to lepton + Higgs. The sum of the 
asymmetries to leptons and to sleptons is about twice larger than the SM asymmetry [56]: 



+ In [1 + 1/x] 



2X 3 / 2 



(5.15) 



The first line of eqn (|5. 13[) corresponds to the diagrams on the first row of figure 15.11 while the second 
line [56, 63, 66, 95] corresponds to the diagram of the second row. This contribution violates the single 
lepton flavours but conserves the total lepton number, and thus it vanishes when summed over a: 



1 



(8tt) [At A] ii 



^Tlml^tA)!,] 2 }^) . 



(5.16) 



As discussed in Sections 15.21 and 19.4.31 the upper bound on the flavoured CP asymmetries e aa can be 
used to obtain a lower bound on the reheat temperature. Here we discuss the upper bound on the total 
CP asymmetry e = J2 a f »o 01 e q n (SHU) [73, 74]: 



e < 



3 (m n 



16tt 



x /3(TO,TO max ,m min ) 



(5.17) 
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where m max (m m in) is the largest (smallest) light neutrino mass, and /? ~ 1 can be found in [162]. 

The interesting feature of the bound (|5.17[) . is that it decreases for degenerate light neutrinos. This 
was used to obtain an upper bound on the light neutrino mass scale from unflavoured leptogenesis 
[135, 182] (discussed in Section f9. 4. ip . and explains the interest in the form of the function j3. However, 
the maximum CP asymmetry in a given flavour is unsuppressed for degenerate light neutrinos [62], so 
flavoured leptogenesis can be tuned to work for degenerate light neutrinos, as discussed in Section f9. 4. II 



5.3 Implications of CPT and unitarity for CP violation in decays 

S-matrix unitarity and CPT invariance give useful constraints on CP violation (see e.g. [24,31,183]). This 
section is a brief review of some relevant results for CP violation in decays. CP transforms a particle £ a 
into its antiparticle which we represent as t a . 

Useful relations, between matrix elements and their CP conjugates, can be obtained from the unitarity 
of the S-matrix S = 1+iT: 

1 = SS f = (l+iT)(l-iT f ) (5.18) 

which implies that iT a b = *T£ a — [TT^] a fc. Assuming that the transition matrix T can be perturbatively 
expanded in some coupling constant A, it follows from 

|T afc | 2 - |T ba | 2 = -2 lm{[TT^] ab T* ba } + \[TT% b \ 2 (5.19) 

that CP violation in a tree process, such as Nj decay, can first arise in the loop corrections. Notice 
that the unstable Ni is being treated as an asymptotic state (the unitary S-matrix is defined between 
asymptotic states); this approximation requires some care, as discussed around eqn (|5.24[) . 
CPT, which should be a symmetry of Quantum Field Theories, implies 

|M(a-» b)\ 2 = \M(b^ a)\ 2 (5.20) 

where iA4(a — > o)(27r) 4 <5 4 (X" pi — X)/ 1 ?/) i s the iT& matrix element from an initial state of particles 
{ai(pi), ...a n (p n )} to a final state of particles {&i(gi), ■ ■■b m (q m )}. In particular, for a Majorana fermion 
N±, which is its own antiparticle, 

\M{N -> £ a <f))\ 2 = \M{fi a -» N)\* (5-21) 

Following many textbooks (for instance, section 3.6 of [183]), one can show from unitarity and CPT 
that the total decay rate of a particle X and its antiparticle X are the same. The unitarity condition 
E m <X|S|6)(6|St|X) = 1 implies 

\M(X - 6)| 2 = ]T \M{b - X)\ 2 , (5.22) 

where the sum is over all accessible states b. Combined with the CPT condition of eqn (|5.20[) (with 
a = X), one obtains, as anticipated, 

£ \M(X - b)\ 2 = \M(X b)\ 2 . (5.23) 

{b} {b} 

(Notice that {6} = {&}.) It is nonetheless possible to have a CP asymmetry in a partial decay rate. In 
the case of Ni, which decays to 4>£ a and <f>i a , the asymmetry of eqn (|5.ip can be non-zero. 

N\ can be approximated as an asymptotic state, for unitarity purposes, if its lifetime is long compared 
to the S-matrix timescale. This timescale can be identified as 1/^/Skm (where Skin is a Lorentz invariant 
measure of the center of mass energy of the process, for instance the Mandelstam variable s for 2^2 
scattering), because, in calculating (for instance) a decay rate, one squares the S-matrix element, using 

f n. m\ / n m \ 

5>-5> $>-!> )Vt, (5.24) 
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where V and r are the volume of the box and the time interval in which the interaction takes place. For 
finite r, there must be an uncertainty in the energy conservation <5-function of order 1/r, so the 5(E) 
makes sense for ^/Skin 3> 1/r. Consequently, unitarity is satisfied for matrix elements with Ni in the 
initial or final state when <C y/s^n ~ M± (the narrow width approximation). 

One must take care, however, to subtract from scattering rates into true asymptotic states the contri- 
bution of on-shcll iVi 's (to avoid double counting) . This is usually done in the narrow width approximation 
(see e.g. ref. [60] for a clear discussion). Then one can check the result by verifying the CPT and unitarity 
constraints. Here we do the converse: use CPT and unitarity to guess what should be subtracted (see 
section 3.8 of [183] for a more complete analysis). 

To see how this works, consider the process <fi£ a — > anything, at 0(A 4 ) in the rate. The possible final 
states are N±, with one-loop corrections, and (j)tp or (pip, at tree level. Then one can make the following 
three observations: 

1. Unitarity and CPT imply (see eqn (|5.23[) with the initial state X — <p£ a ) that there can be no CP 
asymmetry in this total rate: 

\M((f>£ a -> anything)) 2 = \M(tj>l a -> anything)) 2 . (5.25) 

2. For leptogenesis to work, we need e aa ^ which, by CPT, implies that there is a CP asymmetry 
in the partial rates to N\\ 

\M($l a -> Ni)\ 2 - \M{cf>£ a -► A^i)| 2 + 0. (5.26) 

3. From the unitarity constraint (|5.19p . there should be no CP asymmetry to cancel f|5 . 26[) in the 
tree-level scatterings. 

The apparent contradiction arises because on-shell s-channel exchange of N% is included in the scat- 
tering, so we have counted it twice. This on-shell part, also referred to as "Real Intermediate State", 
should therefore be subtracted from the scattering: 

\M(^£ a anything)| 2 = \M{(j>i a N)\ 2 +^ (\M{cj>£ a -> 0£p)\ 2 - \M° s (<f>£ a <^)| 2 



+ \M(cf>£ a -> cj>£ (j )\ 2 - \M os (<f>£ a ^)| 2 ), (5.27) 

where A4° s stands for the on-shell contribution to the amplitude. It is simple to check that the asymmetry 
(|5.25[) vanishes as required if the subtracted matrix element squared, denoted by M! , 

\M'{ct>t a ^4>ip)\ 2 = \M{ct>£ a - 4>ip)\ 2 - |A4 os (cR* -» ^)| 2 , (5.28) 

is taken as follows: 

\M'(<f>£ a ^$tp)\ 2 = \M^£ a ^4>£ (j )\ 2 -\M{cjy£ a - N)\ 2 B^. (5.29) 

where B^j is the branching ratio for N — > <j>£/3- This is the subtracted matrix-element-squared one 

obtains in the narrow width approximation. It will be useful for writing Boltzmann Equations for the 
lepton asymmetry. 

5.4 CP violation in scattering 

Scattering processes are relevant for the production of the Ni population, because decay and inverse 
decay rates are suppressed by a time dilation factor oc M\/T. The N\ = N% particles can be produced 
by s-channel </>-exchange in qt c — ► N£ a and qi c — ► N£ a , and by i-channel </>-exchange in q£ a — > Nt c , 
t c £ a -> Nq, £J C -> Nq and q£ a -► Nt c . 
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In this section, we explicitly calculate the CP asymmetry in scattering processes, for the case of 
hierarchical Nj, and show that it is the same as in decays and inverse decays [64]. This result was 
found in refs. [61,85, 184] for the case of resonant leptogenesis. To introduce CP violation in scattering 
into the Boltzmann Equations, one must correctly include all processes of order h 2 A 4 with the on-shell 
intermediate state A^s subtracted out [185]. This is done in section |6~31 following the analysis of [185]. 

For simplicity, we work at zero temperature, in the limit of hierarchical singlet fermions. This means 
we follow the framework of subsection !5.21 that is, we calculate in the effective field theory with particle 
content of the SM +N\ , where the effects of the heavier N2 and N3 appear in the dimension- five operator 
of eqn (f5~8|) . 

We define the CP asymmetries in AL — 1 scattering (mediated by s- and i-channel Higgs boson 
exchange) as 

~s _ <j(t c q^N£ a )-a(q¥^N£ a ) 
a + a 

H _ a{qN -> ¥£ a ) - 5{qN -> t c I a ) 

a(g£ a -> FN) - d{q£ a -> t c N) 
= ^T~= ' ( 5 - 31 ) 

a + a 

where the cross-sections in the denominator are summed over flavour. The initial state density factors 
cancel in the ratio, so the cross-sections <r, a can be replaced by the matrix elements squared |.M| 2 , 
integrated over final state phase space J dH. Separating the tree and loop matrix elements into a coupling 
constant part c and an amplitude part A, as in eqn (|5.2|) . the CP asymmetry can be written as in eqn (|5.5|) . 
The loop amplitude has an imaginary part when there are branch cuts corresponding to intermediate 
on-shell particles, which can arise here in a bubble on the N line at the N4>£ a vertex, e.g. for s-channel 
Higgs exchange: 



lm{A (t c q -» N£ a )A\{t c q -» N£ a )} = A (t c q -» N£ a ) J A* (t c q -> U^^dR 1 Alit 1 ^ 1 -» N) . (5.32) 

Here (j>' and £'a are the (assumed massless) intermediate on-shell particles, and dH' is the integration over 
their phase space. 

In the scattering process, cq = h t \* la and c\ = 3h t Xif3[m*]f3 a /v 2 , where ht is the top Yukawa coupling. 
The complex coupling constant combination in the scattering processes is clearly the same as in e aa for 
decays discussed in section [5Tl To obtain the amplitude ratio [the second ratio in eqn (|5.5p ]. we take, for 
instance, Aq(N — > <f)£ a ) — uiPlun- After performing straightforward spin sums, we find that it is the 
same for scattering and for N decay, so 

^aa — ^aa = £ aa- (5.33) 

CPT and unitarity are realized in the scattering process (qt c — » N£ a ) in a similar way to inverse 
decays. They should hold order by order in perturbation theory, so we work at order A 2 A 2 h 2 , and define 

\M(qt c ^ N£ a )\ 2 = |.M s | 2 (l + e QQ ), (5.34) 

where | | 2 oc A 2 /i 2 , and |A / J s | 2 e QQ oc A 2 A 2 ft- 2 . At order A 2 A 2 /i 2 , we should also include various 2^3 
tree diagrams without N% in the final state. Following the inverse decay discussion, one can write 

\M(qt c -» A^)| 2 = \M(qt c -» N£ a )\ 2 + ^ [\M(qt c -» £p<\>£ a )\ 2 - \M os (qt c -» £^£ a )\ 2 \ 



+ J2 [\M{qt c -> £^£ a )\ 2 - \M os (qt c -> l^£ a )\ 2 ] 



= \M S \ 2 (1 + e aa ) + \M{qt c -> £cf>£ a )\ 2 - \M S \ 2 (1 + e aa )^^ 
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+\M( q t c -> m a )\ 2 - \m s \ 2 {\ + e aa y-^-i 

= [\M(qt c -> lpcj>£ a )\ 2 + \M(qt c Z^)! 2 ] , (5.35) 

where, in the narrow width approximation, 

\M os (qt c -> ^0C)| 2 = -> M a )| J x (5.36) 

In eqn (|5.35[) . the CP asymmetry e aa has disappeared in the final result, so if we repeat the calculation 
for the CP conjugate initial state, qt c , we should obtain the same result, verifying that a CP asymmetry in 
qt c — * N£ a is consistent with CPT and unitarity. Furthermore, the final result of eqn (|5.35p is reassuring, 
because the unstable state N has disappeared. There is no CP violation in the total rate for qt c — > 
asymptotic (stable) final states, but CP violation in the partial rate to the unstable N is possible. This 
can be relevant to the final value of the baryon symmetry when some of the lepton flavours are weakly 
washed out. 
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6 Boltzmann Equations 



The lepton (and baryon) asymmetry produced via leptogcnesis, can be computed by solving the Boltz- 
mann equations (BE). These describe the out-of-equilibrium dynamics of the processes involving the 
heavy singlet fermions. The aim of this chapter is to derive the basic Boltzmann equations, restricted to 
the non-supersymmetric framework, and to the case where the processes that generate a lepton asym- 
metry involve just the lightest singlet fermion N±. Modifications from supersymmetry (see e.g. [186]) 
arc described in Scction riO.il Possible contributions from heavier singlet neutrinos 7V2,3 [65,187,188] are 
considered in Section Q 0.2. II 

To derive the BE one has to consider a large set of processes, as well as abundances and density 
asymmetries of many types of particles, and the use of notations in extended form can result in rather 
cumbersome equations. Thus, we start in Scction fe.ll bv introducing some compact notation. To make the 
navigation through the details in the following subsections easier, we also present the general structure 
of the final equations. 

Simple Boltzmann equations, taking into account decays, inverse decays, and 2 «-> 2 scatterings 
mediated by N% exchange, are derived in Section [6~2l and are given in eqns (|6.16l) and (|6.28|) . In Section 
16.31 we include scattering processes with Ni on an external leg, and we also discuss 1 <-» 3 and 2^3 
processes, since the CP asymmetries of their off-shell parts must be taken into account for consistency. 
The corresponding BE for the evolution of the Ni density is given in eqn (|6.51| . and for the relevant 
flavour asymmetry in eqns (|6.53p . (|6.54|) and (|6.56[) . 

We emphasize that to reach quantitatively accurate results, one has to take into account (i) some 
relevant interactions that do not involve the Majorana neutrinos, and (ii) flavour effects. These tasks are 
taken in the following sections. 



6.1 Notation 

In this subsection we introduce our notation. A brief introduction to particle number densities and rates 
in the early Universe can be found on Appendix [13] We denote the thermally averaged rate for an initial 
state A to go into the final state £?, summed over initial and final spin and gauge degrees of freedom, as 
(see eqn 113. 14f) 

7 ^ = j(A -> B) . (6.1) 
The difference between the rates of CP-conjugate processes is written as 

A 7 | = 7 |- 7 |. (6.2) 

We denote by n a the number density for the particle a, by its equilibrium density, and by s the 
entropy density (see Appendix [T3|) . We define: 

Y a = — , Va = ^q, Ay a = y a -y a . (6.3) 

Thus, we write all particle densities (Y a ) normalized to the entropy density. To simplify the expressions, 
we rescale the densities Y a by the equilibrium density of the corresponding particle (Y° q — n° q /s). We 
denote the asymmetries of the rescaled densities by Ay a . 

The difference between a process and its time reversed, weighted by the densities of the initial state 
particles, is defined as 

n m 

[A~B] = (Ry^hi (J] Vhhl (6-4) 
»=i j=i 

where the state A contains the particles a±, . . . , a n while the state B contains the particles b%, . . . ,b m . 
We consider only processes in which at most one intermediate state heavy neutrino N% can be on-shell. 
In these cases, a primed notation j'g (and [A <-> B]') refers to rates with the resonant intermediate state 
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(RIS) subtracted. In other words, for the process A — > B, we distinguish the on-shell piece "f os g from 
the off-shell piece 7'^: 

Yi^Tg-Vl (6-5) 
In the simple case where only 2^2 scatterings are considered, the on-shell part is just 

7° s ^ 7 £ 1 < 1 , (6-6) 

where B^ 1 is the branching ratio for Ni decays into the final state B. To include processes of higher 
order in the couplings, eqn (|6 .€>[) needs to be generalized [185]. 

We introduce from the start a set of BE that allow a proper treatment of flavour effects. To do 
that, we write down the BE for the evolution of the density of the heavy singlet fermions N and of the 
asymmetry for a generic lepton flavour a. Below T ~ 10 12 GeV (T ~ 10 9 GeV), reactions mediated by the 
r (/i) Yukawa couplings become faster than the Universe expansion rate, possibly resolving the flavour 
composition of these lepton doublets (see Appendix |15.4|) . As discussed in Appendix [15j if the charged 
lepton Yukawa interactions are fast compared to both H and Tiu , the equations of motion for the lepton 
asymmetry reduce to the BE in the flavour basis. Assuming that the flavour basis does not change during 
leptogenesis, we can work with simple projections onto flavour of all the relevant quantitites. However, 
we still adopt a double-index notation for most of the flavour-dependent quantities, as a reminder that 
they correspond to diagonal elements of matrices in flavour space. For example, we denote the density 
of leptons of flavour a by 

Y£ a = Y la +Y ea , (6.7) 

where Yi a is the density of the two gauge degrees of freedom in £ a . The inclusion of the density Y Ga for 
the right-handed charged lepton e Q is required when (some of) the L-conserving charged lepton Yukawa 
interactions become fast compared to the Universe expansion rate, since in this case they transfer part 
of the asymmetry to the right handed degrees of freedom (see section |5J| . 

Wc define the asymmetry in the densities of leptons and antileptons of flavour a as 

Y ALa = Y£ a - Yl a . (6.8) 

When the rates of charged lepton Yukawa interactions are negligible, the lepton asymmetry is stored only 
in the lepton doublets, and one simply has Y\L a — Yi a ~ Yg a . 

As we explain below, the following asymmetries are particularly useful in the context of leptogenesis: 

Y Aa = ^f--Y ALa) (6.9) 
where Yab is the baryon asymmetry to entropy ratio. For time derivative, we use 

y^f, (6.10) 

z az 

where 

z = Mi/T, H 1 = H(T = M 1 ). (6.11) 
We split the contributions to the evolution equation for Y£ a into three parts: 



Y*a = f Y £ a ) + (Y£ a ) + (y l qq ) . (6.12) 

V / I \ /ll V / sphal 



sphal 

1. (Y£ a )j includes contributions of 0(A 2 ) and of C(A 4 ). It is evaluated in section ROl 

2. {Y^ a )n includes contributions of C(A 2 ft. 2 , A 4 ft. 2 ) and 0(X 2 g 2 , A 4 g 2 ), with g a generic gauge coupling 
constant. It is evaluated in section [6^31 



3. (y£* a )sphai represents the change in the lepton densities due to electroweak sphalerons. 
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As concerns the sphaleron effects, although their precise rates are hard to estimate, it is known that 
below T ~ 10 12 GeV they are a source of rapid baryon number violation, but they leave B — L unchanged. 
More precisely, sphalerons generate the same change in the baryon and lepton number of each generation, 

(iAX a )sphal = gO^AB^phal, (6.13) 

leaving unchanged the charge densities of eqn (|6.9|) . Hence, it is convenient to write an equation directly 
for these quantities. By subtracting from eqn (|6.12[) the analogous equation for the density of antileptons 
Yj* a and by subtracting again the result from the equation that describes the evolution of the baryon 
asymmetry, (Yab)/3 = (YAB)sphai/3, one obtains the evolution equations 

*A«, = - (*AL a ) - (*AL Q ) 7/ (6.14) 

that do not depend on the sphaleron rates. 

The heavy fermions N\ are treated as quasi-stable particles on the time scale of the Universe expansion. 
This is justified by the fact that leptogenesis requires that the N% lifetime is of the order of the expansion 
time -fff 1 . The final baryon asymmetry depends on the density of the neutrinos N\ as a function of time, 
so a Boltzmann equation for Yn 1 is needed. It is convenient to split also this equation into two parts: 



Y, 



(Y Nl ) i+ (Y Nl ) n . (6.15) 



Ni 

\ / i \ / jj 

The term (Yn^i includes the contributions of terms up to 0(A 2 ) and is evaluated in section IB~2l The 
term (Yzvi)// includes the contributions up to 0(X 2 h 2 ) or 0(X 2 g 2 ) and is evaluated in section HT31 

6.2 The 0(A 2 ) and C(A 4 ) terms 

This section aims to obtain the basic BE which depend only on the neutrino Yukawa coupling A, including 
only the terms (Yal q )/ of eqn (|6.14|) and (Yp^)/ of eqn (|6.15|) . The SM gauge interactions are assumed 
to be fast, which ensures kinetic equilibrium for the particle distributions. All other SM interactions, in- 
cluding sphalerons, are neglected. The processes of Ni decay and inverse decay and two-to-two scattering 
mediated by the N^s are included. However, the latter is important mainly to subtract real intermediate 
states, but its effects are small in the temperature range T < 10 12 GeV and in a first approximation can 
be neglected (see the appendix of [62] for a brief discussion, and [59] for a detailed one). 
The Boltzman equation for Yn 1 , including only decays and inverse decays, is given by 

(^) = E { fat ~ N ^ + ^ N ^} 

= - 53 [y^ + t£) ~ v+vhiN? ~ vfVifiirtff - ~ (y*i - 1) 7w^2, (e.ie) 


where in the last expression we have approximated the cj> and lp number densities with their equilibrium 

densities, neglecting small corrections cx — it^^ivj^ ~ Hep) ^ na ^ are second order in the small 

quantity epp. We use 



for the thermally averaged two body N\ decay rate. 

To obtain the BE for the lepton asymmetry we should work to order A 4 , because e aa — C(A 4 ). At 
this order, the doublet leptons participate in the following three types of processes: 

(i) 1 *-* 2 processes: the (CP violating) decays Ni — > € Q 0, and inverse decays l a 4> — > N\\ 
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(ii) 2<->2 scatterings mediated by s-channel N% exchange: £ a 4> <-» £p(f> with (3 ^ a, and £ a <j) <- 

(Hi) 2^2 scatterings mediated by t- and u-channel iVi exchange: (f>4> <-> ^ Q ^/3 and 00 «-> £ot£ j3 • (Neglect- 
ing provisionally thermal effects, no RfS can appear in this case, and hence there are no on-shell 
contributions to be subtracted.) 

Accordingly, at this order the evolution equation for the density of the lepton flavour a reads 

(V), - (Yr)^ 2 + (yr) N 2 i 2 + ( Y r)Z 2 < ^ 

where 

= [Ni^£ a 4>], (6.19) 

= B^«4^+£N«W]', (6.20) 

( Y r) N 2 2 = ^{[<M"eJp] + 0. + s af i)[M*->t a tp]}. (6.2i) 

2^2 

We first focus on (i) of eqn (|6.19p and (ii) oi eqn (|6.20p that give rise to the source term for the asymmetry. 
The processes (Hi) of eqn (|6.21[) contribute only to the washout (at C(A 4 )) and will be added at the end 
of the section. As discussed in section [5.31 some care is required in combining (i) and (ii): In the 2 <-> 2 
scattering, the contribution from on-shell s-channel iVi-exchange is already accounted for by the decays 
and inverse decays. This contribution should be removed, to avoid double-counting, by using a subtracted 
\A4\ 2 (see eqn I5.29|) . as indicated by the primed notation in eqn. (|6.20[) . That is, we include only the 
off-shell part of the scatterings rate density 

7 i x = 7- - — Tw B- - 6.22 

that has a CP asymmetry of the same order as the CP asymmetries of decays and inverse decays. The 
scatterings £ a (f> <-> £p (f> (with (3 ^= a) are treated in a similar way. 

The equations for the comoving number densities of the lepton doublet, and of the anti- lepton doublet, 
are therefore 

{yA = [Nx~l a <f>]+Yyi$"la$ + Yy t rf*+ l * , l > ]', (6-23) 

1 13 

1 f3 f3 

Note that the sum in the last two terms of the two equations has been extended to include the contribu- 
tions from (3 = a which cancel in the difference between the process and its time reversed. Eqns (|6.23[) 
and (|6.24p can be written more explicitely as 



= ^x^+E(^v£{+i W 7t3-^EW{+'rtJ). (6-25) 


- ^<V + E(^7t! + ^7tj)-^E(7t^Vg). (6.26) 





Using eqns (16. 5[) and (|6.6p and CPT, we rewrite eqn (|6 . 25[) as follows: 







\-Y{B?\+B^-) 



/3 



V<t>Vl 

' 



E Kl + l£f] + E [w^SJ + wS] , (6-27) 
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and similarly for the analogous terms in eqn. (|6.26p . Some comments are in order with regard to 
eqn (|Q7)1 : 

1. At the order in A we are working in this section Y^p(B^ + Bj L) = 1. Consequently, the second 
term in the first line vanishes. 

2. As concerns the second line of eqn. (|6.27p , the first term is proportional to j{£ a 4> ~^ anything), 
while the second term is the sum of 7(anything — > <p£ a ) and terms that are second order in the 
asymmetry. We know from eqn (|5.23|) that there can be no CP asymmetry in differences of the 
form 7(A <-> anything) — y(X «-> anything). Consequently, these terms do not contribute to the 
CP asymmetry. 

3. The source term (oc e aa ) in the BE arises from the first term of (|6.27p and from its analog in Y^ q . 
We can set the sum of the bracketed branching ratios to 1, because corrections to this approximation 
are second order in the asymmetry. We thus obtain the correct behaviour: no asymmetry can be 
generated in thermal equilibrium (yNi = 1)- 

4. The last term in eqn. (I6.27P does not contribute to the washout of yi a . Nevertheless this term 
should not be dropped since, as it will become clear in section f9.3( it does induce washouts for the 
charge density Y& a in eqn. (|6.9D . 

By subtracting (|6.26|) from (|6.25|) . and including the contribution of the t and u-channel processes 
(in) of eqn. (|6.21[) . we obtain the complete BE at 0(A 4 ): 

s X - ^ / /. <h . 

- 7- 



(*A ia ) = (y^-l)A^-(A Wa + A^)X;(7^+- 

+ E( A ^+ A ^)(^^^) + (^^. (6-28) 



The term 

(W)^' = - E [(! + M( A ^ + A ^ + 2A vMt + {Aye ° - A y^Ht] (6 - 29) 







is straightforwardly obtained by subtracting from eqn (|6.21|) the analogous eqaution for £ a . In eqns 
(|6.28[) and (|6.29[) we approximate the washout rates with their tree level values, and we linearize in the 
asymmetry densities Ay. 

Note that there are no resonant contributions to the washout from the second line of eqn (|6.28[) since 
the on-shell parts contained in 7»' 3 f and 7 f _ cancel in the difference. The washout term in the first 

tot™ ^-a'P 

line of eqn. (|6.28[) can be written as the sum of resonant (0(A 2 )) and non-resonant (C(A 4 )) parts. By 
dropping all the subleading non- resonant terms, we obtain an approximate expression, valid at C(A 2 ): 



fc4 



(y Nl - 1)A7& - (Ay la + A W ) 7 £ 



£ a rp 

1 rn 

7Ni— 2, (6.30) 



(y Nl - l)e aa - \{Ay ia + Ay^)^^ 

2 TO 



where e aa is defined in eqn (|4.2[) . m aa /m = T(N\ — > £ a <j>, £ a 4>)/Tjj, and 7 at 1 ^2 is defined in eqn (|6.17[) . 
Using eqns (|6.10p , (| 13 . 2[1 and (|13.18p , and provisionally neglecting the contribution to the washout of the 
Higgs asymmetry Ay^ (see section [8|), eqn. (|6.30j) can be written in a more explicit form: 



d (Y,r \ - Z ((Y« Y°^ Kl ^ ' 



(YalJj = ir ( (Y Nl - Y^y-^le aa - ^ Kl (z)Y ALa ^ ) T D . (6.31) 



3G 




Figure 6.1: Diagrams for various 2^2 scattering processes: (a) scatterings with the top-quarks, (b), (c) 
scatterings with the gauge bosons (A = B,Wi with i = 1,2,3), (d) AL — 2 scatterings mediated by Aq. 

6.3 The 0{h 2 t X 2 ) and 0{h 2 t X 4 ) terms 

In this section, we include processes involving the top Yukawa coupling h t . Processes involving gauge 
bosons can be included in a similar way and we add them in our final expressions. 

We denote the left-handed third-generation quark doublet by q 3} and the St/(2)-singlet top by t. The 
inclusion of 1 <-» 3 decays and inverse decays such as Ni <-» £ a q3t, and of N\£ a <-> q-$i scatterings mediated 
by Higgs exchange, follows lines analogous to those presented in the previous section. For the 0(h^X 2 ) 
contributions to the evolution of the Ni density, we obtain: 

(Yn 1 ) h = -{y Nl 1) [7A^3 + Ttip 2 ] • (6-32) 

Here, 
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is the contribution from decays into three-body final states, while 



2^2 = V- / ATiifl Njp N iq3 N iq3 N lt , NJ \ , g M) 

13 

is the contribution from Higgs mediated scatterings: the first two terms correspond to s-channel Higgs 
exchange, while the other four (that are all equal at leading order) correspond to t- and u-channel Higgs 
exchange (see fig I6.1f a)). 

Regarding the evolution of the lepton asymmetries, the derivation of the BE is more subtle. Once 
we include the CP violating asymmetries in 1 <-> 3 (inverse) decays, like N\ «-> £ a qzt, and in 2 <-> 2 
scatterings, like N\£ a <-> qgt, we must include also the asymmetries of various off-shell 2^3 scatterings, 
which contribute to the source term at the same order in the couplings. Accordingly, we write the term 
(Y^ a ) n of eqn (jQ2)) as follows: 



where 



Y£ a ] = [Ni «-> £ a q 3 t] + [q 3 t~ Nit a ] + [NJ «-> q 3 £ a ] + [N lQ3 «-> tt a ] ; (6.36) 



2«2 

/ . \ sub — 



13 

+ J2(l+8 a p)[q 3 i^£ a <t>£p}'; (6-37) 

/3 

iV t 



+ { W «-> Va] + «"» V"] + [q 3 4> Vat] + N Va©]} 

+ ^(1 + 5 Q/3 ) { «-» 4?3^] + [<fet? «-> *a*/s] + [930 ~ • (6.38) 

18 

As in the previous section, the asymmetries in the off-shell 2^3 rates in eqn (|6.37[l can be estimated 
by relating them to the asymmetries of the corresponding on-shell parts. However, for 2 <-> 3 scatterings, 
the definition of the on-shell part is more subtle, because after a real Aq is produced in a collision, it has 
a certain probability to scatter before decaying. Consider, for example, the process £/3<f>q3 — > t£ a . The 
contribution to this process from the exchange of an on-shell Aq corresponds to the production process 
lp4> ~ ¥ Ai, followed by the scattering Aq + q 3 — > t£ a that is mediated by a Higgs in the t channel. 
Processes of this kind can generally be written as AX — » Y where A denotes a possible state to which 
Ai can decay. The corresponding on-shell rate is then 

1 o S AX =1 A i pN lX> (6 39) 

where Py lX is the probability that Aq scatters with X to produce Y. Processes in which the on-shell Aq 
can disappear only by decaying (as for example £p4> «-> £ a <l3t or £pt «-> £ a <t lC lz) can generally be written as 
A — > _B or as X — > BY", where both A and B denote possible final states for Aq decays. The corresponding 
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on-shell rates are 



7 b = 

7 BY — 



N 1 



iNi Y "b ■ 



(6.40) 
(6.41) 



Note that because of the fact that in the dense plasma Aq can be scattered inelastically before decaying, 
as described by eqn ()6.39|) . the quantities Pg 1 in eqns (|6.40[) and ()6.41j) differ from the usual notion of 
branching ratios at zero temperature. In particular, scattering rates should also be included in normal- 
izing properly the decay probabilities. The quantities P b a then denote the general probabilities that Ni 
contained in state a ends up producing a state b. In the case under discussion, we have, for example, 



pNll a 
q 3 t 



>y Nl 

ie a <t> 

Tall ' 
_ 7 93« 



Tall 



>e a g 3 t 

Tall 



P 



N iq3 



II', 



^193 
Tall 



pNtt 
r q 3 Z a 



7g 3 e a 

Tall 



(6.42) 



with similar definitions for the probabilities of the CP conjugate processes. The probabilities are normal- 
ized in terms of the sum of all the rates: 



Tall 







, N-l . iVi , Ni , Nil a , iVi^ , iVig 3 

+ ^ + 7^5 3 t + T^f + 7 S3 j + 7 & i + T^f 



•7U +7^+7^)- (6-43) 



To the order in the Yukawa couplings that we are considering, the unitarity condition for the sum of the 
branching ratios of N\ into all possible final states, J2y By 1 



1, is then generalized to Y Py lX = 1- 



In other words, the probabilities for all the possible ways through which Nx can disappear add up to 
unity. 

To include the new sources of CP asymmetries, we now need to subtract from eqns (|6.36I6.37|) the 
analogous equations for Y£ a . For the source term, we obtain: 



Y ALa y ir (Y ALa )y 2+ (Y ALa ) 



s, sub 
2<->3 



(6.44) 



where we neglect the CP asymmetries of the 2 ^> 3 processes with Ni exchanged in the i-channel 
eqn. (|6.38[) that are of higher order in the couplings. For the first term in the r.h.s. of eqn (|6.44[) we have 



(M 



l<-»3 
2«2 



a 7 . 



t a q 3 t 



At 



Nit 

1 93 



'q 3 t 



(6.45) 



Eliminating the subtracted rates by writing their CP asymmetries as minus the CP asymmetries of the 
on-shell rates, and keeping terms up to 0(X 4 hf), we obtain for the second term in eqn (|6.44D 



s, sub 

YAL a )_ =-2A 7 , 



N t 



1 - 







P 



Ni 



,Ni 



At 



Ni 

e a q 3 t 



1 i 



Nii a 

q 3 t 



(6.46) 

Note that, at 0(X 2 hj ), the sum of the branching ratios in eqn (I6.27|) of the previous section is not unity, 
so the first term in that equation does not have the correct thermodynamic behaviour cx (j/jVi — 1); and 
the second term does not vanish. By keeping track carefully of the relevant higher order terms, the source 
term arising from the difference between (|6.27|) and the analog equation for Yj reads 



VNi + 1 



\pNi 
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where J^fsi^^p + -^i^) ^ However, the first term in eqn (|6.46p combines with eqn (|6.47[) to yield for 
the source term involving A the correct behavior oc (yN ± — !)■ Summing up eqns (|6.45[) . (|6.46p and 
(|6.47p . we obtain the final expression for the source term that holds at 0(A 4 /i(): 
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(6.48) 



Regarding the washouts, we neglect the contributions from eqns (|6.37p and (|6 . 38[) since they both 
involve non-resonant 2^2 scatterings that are of higher order in the couplings. We retain only the 
(D(\ 2 ht) contributions of eqn ()6.36|) . Subtracting from eqn (|6.36p the analogous equations for Z Q , we 
obtain the relevant washout term: 

(^)X^Ob (6 - 49) 

where 
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Note that while the contributions from the RIS-subtracted 2 <-* 3 scatterings in eqn (|6.37p must be taken 
into account to obtain the correct form of the source term, neglecting them in the washouts, as we do in 
eqn (|6.50p . does not have a significant effect on the numerical results. 

The inclusion of N\ decays into three body final states is qualitatively required if we want to take into 
account all processes of the same order in the couplings, and to incorporate consistently 2 ^> 3 scatterings 
where the on-shell piece involves a 1 — > 3 decay (for example, tpcj) — > iVi — > £ a q3t). The quantitative 
impact is, however, rather small: while for zero temperature this decay has a large enhancement related 
to small momentum-values for the Higgs propagating off-shell, for the relevant temperatures, the finite 
value of the thermal mass prevents this enhancement, and the decay rate is below 6% of the two-body 
decay rate [185]. We thus neglect the washout term of 3 — ► 1 inverse decays, and the contribution of the 
three-body decay CP asymmetry Ajj^-^ to the source term in eqn (|6.48|) . 

Following the same procedure outlined above, it is possible to include in the BE other relevant pro- 
cesses, such as those involving the gauge bosons [60,85,189]. With all the subdominant terms neglected 
and with the effects of the gauge bosons included, the simplified expression of the BE for the evolution 
of Yjvj reads: 



Y Nl = -(VNx - 1) [lN^2 + 7to7 2 + 7a 



(6.51) 



where A — Wi or B for SU(2) and U(l) bosons respectively. The term involving the gauge bosons is 
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where a sum over the gauge boson degrees of freedom in all the rate densities is understood. In eqn (|6.5ip 
we neglect three-body decays involving the gauge bosons like that are suppressed by phase space 

factors. We also neglect the contributions to the washouts from gauge bosons 2^3 processes. 

We can finally write the simplified evolution equation for the charge-densities Ya q [defined in eqn 
in terms of the source and the washout terms: 
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The source term is given by 
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In the second line we use the approximate equality between the scatterings and the decay asymmetries 
that was discussed in section f5. 41 for example, 



AV 33 ' A>y Nl 

7Ar^ Q lt a <j> 



(6.55) 



The washout term is given by 

(j^T = E[( A ^+ A ^)(i1+<1)+( A ^+ A ^)(i1-<1 



+ E [(! + <WX A ^ + A W/> + 2A V4-hf^ + " A W/»)7^ 

+ (y Wl A Wa - Ay 93 + Ay t )7^/ Q + [2 Ay ia - (y Nl + l)(Ay q3 - Ay t )] ^ 

(y Nl Ay ia + Ay^f" + (y Nl A y<j> + Ay la ) 7 ^f + (A y<t> + Ay la (6.56) 

In the third line we use the equality of the t- and u-channels top-quark scatterings to set J^ 1 / = 7 t ^ 193 - 
Confronting eqn (|6.54[) with eqn (|6.5ip . we learn that the BE for the evolution of the A Q charge 
density can be written as follows: 

Y Aa = Y Nl e aa + (y Ac ,Y ■ (6.57) 

Two comments are in order: 

1 . The washout term of eqn (|6.56p depends on the density asymmetries of various particle types which 
evolve with time. In principle, we need to know their time evolution in order to solve eqn (16.57p . 
However, as discussed in Sections 18.21 and 19.31 and, in more detail, in Appendix [T31 the chemical 
equilibrium constraints from fast SM reactions always allow one to express all the relevant density 
asymmetries Ay a in terms of the Ya - Doing that, one obtains a closed system of differential 
equations involving only the flavoured charge densities. 

2. Eqn (|6.57j) shows that if washouts were neglected (and the value of e aa assumed independent of 
the temperature; see Section [7]), the final value of Ya q would be simply proportional to the initial 
value of Yn 1 . Therefore, in case that the A-interactions are the only source of populating the N\ 
degree of freedom, the final asymmetry would vanish if not for the presence of the washouts. 
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7 Thermal Effects 



At the high temperatures at which lcptogencsis occurs, the light particles involved in the processes relevant 
for the generation of an initial lepton number asymmetry are in equilibrium with the hot plasma. The 
thermal effects give corrections to several ingredients in the analysis: (i) coupling constants, (ii) particle 
propagators (leptons, quarks, gauge bosons and the Higgs) and (hi) CP-violating asymmetries. These 
three effects are discussed in turn in the following three subsections. A dedicated study of thermal 
corrections to leptogenesis processes with a discussion of the leading numerical effects can be found 
in [60]. 

7.1 Coupling constants 

A detailed study of gauge and Yukawa couplings renormalization in a thermal plasma can be found 
in [190]. In practice, it is a very good approximation to use the zero-temperature renormalization group 
equations for the top-quark Yukawa coupling and for the gauge couplings, with a renormalization scale 
A ~ 2irT [60] . The value A > T is related to the fact that the average energy of the colliding particles in 
the plasma is larger than the temperature. 

The renormalization effects for the neutrino couplings are also well known [191,192]. In the non- 
supersymmetric case, to a good approximation these effects can be described by a simple rescaling of the 
low energy neutrino mass matrix m(p) = r ■ to, where 1.2 ^ r 1.3 for 10 s GeV^S /i ^ 10 16 GeV [60]. 
Therefore, RG effects on neutrino couplings can be accounted for by increasing the values of the neutrino 
mass parameters (for example, to) as measured at low energy by « 20% — 30% (depending on the 
leptogenesis scale). In the supersymmetric case one expects a milder enhancement, but uncertainties 
related with the precise value of the top- Yukawa coupling can be rather large (see fig. 3 in ref. [60]). 

7.2 Decays and scatterings 

In the thermal plasma, any particle with sizeable couplings to the background states acquires a ther- 
mal mass that, modulo the renormalization of the relevant couplings, is proportional to the plasma 
temperature. Consequently, decay and scattering rates are modified. The diagrams corresponding to the 
relevant leptogenesis processes for which these corrections should be estimated are given in fig. 15. 11 (1 «-> 2 
processes) and in fig. 16.11 (for the 2 <-> 2 scatterings). 

Thermal corrections to particle masses have been thoroughly studied in both the standard model and 
the supersymmetric standard model [193-198]. The singlet neutrinos have no gauge interactions, their 
Yukawa couplings are generally small and, during the relevant era, their bare mass is of the order of 
the temperature or larger. Consequently, to a good approximation, corrections to their masses can be 
neglected. We thus need to account for the thermal masses of only the lepton doublets, the third genera- 
tion quarks, the Higgs and the gauge bosons (and, in the supersymmetric case, also their superpartners). 
Explicit expressions for the thermal masses that enter the relevant leptogenesis processes are collected in 
appendix B of [60]. For the following qualitative discussion, it is enough to keep in mind that, within the 
leptogenesis temperature range, to^(T) <; m g3j t(T) ^> mi(T). The most important effects relate to four 
classes of leptogenesis processes: 

(i) Decays and inverse decays. Since thermal corrections to the Higgs mass are particularly large 
(m^(T) w 0.4 T), decays and inverse decays become kinematically forbidden in the temperature range in 
which TO0(T) — mi(T) < Mn 1 < m${T) + mi(T). For lower temperatures, the usual processes N\ <-> £cf> 
can occur. For higher temperatures, the Higgs is heavy enough that it can decay: <fr <-> tN\. A rough 
estimate of the kinematically forbidden region yields 2 ;$ T/M\ ^5 5. The important point is that these 
corrections are effective only at T > M\. In the parameter region to > 10 _3 eV, that is favored by 
the measurements of the neutrino mass-squared differences, the N\ number density and its L-violating 
reactions attain thermal equilibrium at T w M% and erase quite efficiently any memory of the specific 
conditions at higher temperatures. Consequently, in the strong washout regime, these thermal corrections 
have practically no effect on the final value of the baryon asymmetry. 
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Figure 7.1: Comparison between scattering rates in the Standard Model with (solid lines) and without 
(dashed lines) thermal corrections. We use m = 0.06 eV and M\ = 10 10 GeV. (a) 7jy t : the rate density of 
the i-channel iVi-exchange scattering ££ <-* (f>4>, normalized to neH. (b) 7^ op and 7^° p : the rate densities 
of the Higgs exchange scaterrings in, respectively, the s-channel {q^i «-> IN) and the t- and it-channels 
((73 N <-> £i and t AT <-» ^93), normalized to un 1 H. (Figures adapted from ref. [199].) 

(ii) AL = 2 scatterings. A comparison of the scattering rates for ££ <-> (/) with and without thermal 
corrections is given in fig. 17. 1( a) (adapted from [199]). The reaction densities are computed for to — 
0.06 eV and Mi = 10 10 GeV and are plotted as a function of decreasing temperature (increasing z — 
M\/T). It is apparent that for this process thermal effects are sizeable only in the high temperature 
range z < 1. For £<j) <-> £(f> scatterings, a new resonant contribution can appear at high temperatures 
due to the fact that A^ can go on-shell when exchanged in the u-channel [60]. As regards the off-shell 
contributions to this process, they are affected by thermal corrections in a way similar to £<p «-> £(f>, that 
is mainly at Mi/T < 1. We conclude that, for the AL = 2 rates, thermal corrections are sizeable only at 
high temperatures. In the theoretically preferred regime, to > to*, the related effects can be neglected. 

(iii) AL = 1 scatterings with top-quarks. Comparisons between the thermally-corrected and -uncorrected 
rates of the AL — 1 top-quark scattering 7^ op = 7(93 t <-> £ Ni) with the Higgs exchanged in the s-channel, 

and of the sum of the t and it-channel scatterings 7^° p = 7(^3 A^ <-* £t) + j(i Ni <-> ^93), are given 
in fig. I7.1f b) (adapted from [199]). In the case of 7^° p , mild corrections are present only at high tem- 
peratures. However, in contrast to the previous cases, the most relevant corrections to 7^° p appear at 
low temperatures, reducing the scattering rates and suppressing the corresponding contributions to the 
AL = 1 washout. This peculiar situation arises from the fact that in the zero temperature limit there is 
a large logarithmic enhancement ~ I^Ma^/to^) from the quasi-massless Higgs exchanged in the t- and 
u-channels. This enhancement disappears when the Higgs thermal mass to^(T) ~ T ~ Mn x is included. 

(iv) AL = 1 scatterings with the gauge bosons (see figs. I6.1f b) and !6.lT c)). The inclusion of thermal 
masses is required to avoid IR divergences that would arise when massless £ (and (f>) states are exchanged 
in the t- and it-channels. A naive use of some cutoff for the phase space integrals to control the IR 
divergences can yield incorrect estimates of the gauge bosons scattering rates and would be particularly 
problematic at low temperatures, where gauge bosons scatterings dominate over top-quark scatterings. 

7.3 CP asymmetries 

As discussed in section HT31 CP asymmetries arise from the interference of tree level and one-loop am- 
plitudes, when the relevant couplings involved have complex phases, and the loop diagrams have an 
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absorptive part. This last condition is satisfied whenever the loop diagram can be cut in such a way that 
the particles in the cut lines can be produced on shell. In the N\ decay asymmetry at zero temperature 
this is guaranteed by the fact that the decay products, the Higgs (j> and the lepton doublet I, coincide 
with the states circulating in the loops. However, at the high temperatures at which the Nx's decay, 
the Higgs and the lepton doublets are in equilibrium with the hot plasma, and their interactions with 
the background particles modify the CP asymmetries and introduce a dependence on the temperature: 
e — ► e(T). Thermal corrections to CP asymmetries arise from various effects: 

i) The possibility of absorption and re-emission of the loop particles by the medium requires the use 
of finite temperature propagators for computing the absorptive parts of the Feynman diagrams. 

ii) The stimulation of decays into bosons and the blocking of decays into fermions due to the dense 
background requires a proper modification of the density distributions of the final states. 

Hi) Thermal motion of the decay particles with respect to the background breaks the Lorentz symmetry 
and affects the evaluation of the CP asymmetries. 

iv) Thermal masses should be included in the finite temperature resummed propagators, and also 
modify the fermion and boson dispersion relations. Their inclusion yields the most significant 
modifications to the zero temperature results for the CP asymmetries. 

The first three effects were investigated in [200] . A rather complete account of thermal corrections to the 
decay CP asymmetries, that includes also the effects of thermal masses, can be found in [60]. In principle, 
at finite temperature, there are additional effects related to new cuts that involve the heavy -^2,3 neutrino 
lines (see fig. 15.11) . These new cuts appear because the heavy particles in the loops may absorb energy 
from the plasma and go on-shell. However, for hierarchical spectrum, M 2j 3 3> Mi, the related effects are 
suppressed to a negligible level by a Boltzmann factor exp(— M2 3/T) that, at the temperatures relevant 
for the Ni decays, is tiny. 



7.3.1 Propagators and statistical distributions 



The real time formalism of thermal field theory [201,202] can be used to compute the particle propagators 
at finite temperature. In this formalism, ghost fields dual to each of the physical fields have to be 
introduced, and consequently the thermal propagators have 2x2 matrix structures. For the one-loop 
computations of the absorbtive parts of the Feynman diagrams, the relevant propagator components are 
just those of the physical fields, that for fermions (£) and bosons (4>) are: 



Se(p,m e ) 
D<t>{p, m ) 



p 2 - mj + «0 H 



P 2 ~ m s + 



2tt ri£ S (p 2 — mj ) 



+ 2irn^S(p 2 ~ m^) 



(7.1) 
(7.2) 



The leading effects in i) are proportional to the factor —rif + n^ — Iniyi^, where n^A — [exp(£ , £^/T± 1] — . 
This factor vanishes when the thermal masses of the leptons and of the Higgs are neglected, because the 
Bose- Einstein and Fermi-Dirac statistical distributions depend on the same argument, Eg — E$ = Mi/2, 
and consequently the thermal correction to the fermion propagator (ne), the thermal piece of the boson 
propagator (n^) and the product of the two thermal corrections (ngn^) cancel each other. This can be 
interpreted as a complete compensation between stimulated emission and Pauli blocking. As regards 
the effects in ii), they lead to overall factors that cancel between numerator and denominator in the 
expression for the CP asymmetry e. 

In the supersymmetric case, the situation is more subtle. First, the singlet neutrino N± decays not 
only to the standard £<j> final state, but also to their superpartners: £(j>. Given that both decay channels 
contribute to the imaginary part of both decay modes, the same cancellation as in the previous case 
occurs (when thermal masses are neglected). Second, a new source of lepton asymmetry comes from the 
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scalar neutrino N\ that decays both into final fermions £(/>, and into final bosons £<f>. Thermal effects 
modify the asymmetry in each channel. One-loop diagrams with the £<f> bosons contribute to the CP- 
asymmetry for decays into fermions, while one-loop diagrams with internal £<p fermions contribute to the 
CP-asymmetry for decays into bosons. As a result, the finite temperature propagator corrections to the 
partial asymmetries of the single decay channels do not vanish [200] . When the statistical functions that 
block and stimulate the final state emission (corrections of type ii)) are taken into account, the branching 
fractions into fermion and boson final states differ. When the CP-asymmetries of the two channels are 
summed up, the effects of the two types of corrections i) and ii) compensate each other, and the zero 
temperature result is again reproduced [200] . 

7.3.2 Particle motion 

We have seen that when a large hierarchy M 2 ^ 3> Mj is assumed, the particle thermal masses are 
neglected, and the decaying particle is considered at rest in the thermal bath, there are no thermal 
corrections to the zero temperature results for the CP-asymmetries. However, since the decaying particle 
is moving with respect to the background with velocity /?, due to their statistics, the fermionic decay 
products are preferentially emitted in the direction anti-parallel to the plasma velocity (for which the 
thermal distribution is less occupied), while the bosonic ones are emitted preferentially in the forward 
direction (for which stimulated emission is more effective). This induces an angular dependence in the 
decay distribution at order 0(0). In the total decay rate the 0((3) anisotropy effects are integrated 
out, and only O{0 2 ) effects remain [200]. Therefore, while the inclusion of the effects of the thermal 
motion of the decaying particle do modify the zero temperature results, these corrections are numerically 
small [60, 200] and generally negligible. 

7.3.3 Thermal masses 

When the finite values of the light particle thermal masses are taken into account, the arguments of 
the Bose-Einstein and Fermi-Dirac statistical distributions are different. It is a good approximation [60] 
to use for the particle energies Ee^ = Mi/2 =p (m| — mf)/2Mi. Since now Eg ^ E^, the prefactor 
n£ — n ( f > + 2ri£n^ that multiplies the thermal corrections does not vanish anymore, and sizeable corrections 
become possible. The most relevant effect is that the CP-asymmetry vanishes when, as the temperature 
increases, the sum of the light particles thermal masses approaches Mi [60]. This is not surprising, since 
the particles in the final state coincide with the particles in the loop, and therefore when the decay 
becomes kinematically forbidden, also the particles in the loop cannot go on the mass shell. The same 
happens in the supersymmetric case for the CP-asymmetry in N± decays. However, this is not the case 
for the decays of the scalar neutrino Ni — > £<j> and Ni — > £<f> for which the particles running in the loop 
arc different from the particles in the final states. Since thermal masses are larger for the scalars than 
for the fermions, the N\ CP-asymmetry vanishes when the decay N\ — > £(f> is still kinematically allowed. 
The relevant analytical expressions for this case and detailed numerical results can be found in [60] . 

When the temperature is large enough that the Higgs can decay (see section [72]) . there is a new source 
of lepton number asymmetry associated with the decay processes <j> — ► £N\ . The CP-asymmetry in Higgs 
decays can be up to one order of magnitude larger than the CP-asymmetry in N\ decays [60] . This is 
mainly due to a kinematical suppression of the tree level decay rate appearing in the denominator of 
that is roughly proportional to the thermal mass difference m \ — mj . While this represents a dramatic 
enhancement of the CP-asymmetry, is non- vanishing only at temperatures T ~ 5M\ , when the 
kinematical condition m^(T) > me(T) + Mi is satisfied. Therefore, in the strong washout regime, no 
trace of this effect survives. On the other hand, rather large A couplings are required in order that Higgs 
decays can occur before the phase space closes: the decay rate can attain thermal equilibrium only when 
m <; {T^/Mifrn* m*, and therefore, in the weak washout regime (to ^ to*), these decays always 
remain strongly out of equilibrium. This means that only a small fraction of the Higgs particles have 
actually time to decay, and the lepton-asymmetry generated in this way is accordingly suppressed. 

In summary, while the corrections to the CP-asymmetries can be significant at T ^5 M\ (and quite 
large at T ^ Mi for Higgs decays), in the low temperature regime, where the precise value of e plays 
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a fundamental role in determining the final value of the baryon asymmetry, there are almost no effects, 
and the zero temperature results still give a reliable approximation. 

Before concluding this section let us mention that effects similar to the ones described above can be 
expected also for the CP-asymmetries of scattering processes, like top-quark or gauge-boson scatterings. 
These asymmetries were considered in [61,64,85] in the approximation in which they are proportional to 
the CP-asymmetry in decays, and were further analyzed in [185] by going beyond this approximation, 
but still in the zero temperature limit. An important difference is that while 2^2 scatterings are always 
kinematically allowed, the absorptive part of the one-loop diagrams vanishes at the same thresholds when 
the decay CP-asymmetries vanish. This can have the peculiar effect of thermalizing the N\ degree of 
freedom without producing an associated lepton asymmetry. To our knowledge, a study of the effects of 
thermal corrections to CP-asymmetries in scattering has not been carried out yet. 
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8 Spectator Processes 



8.1 Introduction 

During leptogenesis, various processes can modify the densities of particle species. Some of these, such as 
the heavy neutrino decays or various interactions that washout the lepton number, occur on a time scale 
comparable to the expansion rate of the Universe, and hence should be accounted for via appropriate 
Boltzmann equations. Other processes can be very fast (depending on the temperature considered) 
and their effect is to impose certain relations among the chemical potentials of various particle species. 
These processes include the gauge interactions, Yukawa interactions involving the heavier fermions, and 
the electroweak and QCD non-perturbative 'sphaleron' processes. They are called 'spectator processes' 
because they do not change lepton number directly. Instead, they affect it indirectly, by changing the 
densities of the lepton doublets and of the Higgs on which the rates of washout processes depend. The 
issue of spectator processes and their effects on leptogenesis was first raised in ref. [203] . A rather complete 
analysis of the numerical relevance of each spectator processes can be found in ref. [204] . The main results 
of this section can be red off from the last column in Table [TJ where the effects of various spectator 
processes are quantified as a percentage variation of the final asymmetry resulting from leptogenesis. 

We work in the scenario in which the singlet neutrino masses are hierarchical, Mi <C A/2,3, and the 
lepton asymmetry is generated mainly via the CP and lepton number violating decays of the lightest 
singlet neutrino Ni . We restrict the discussion to the non-supersymmetric case, since no qualitative new 
features appear in the supersymmetric case. 

Several spectator processes become relevant only in the temperature regime in which lepton flavour 
effects are also important. In particular, reactions mediated by the r Yukawa couplings become faster 
than the expansion rate of the Universe below T ~ 10 12 GeV (10 9 GeV). Then it becomes important to 
know the flavour composition of the lepton asymmetry. For simplicity, in our discussion of spectator 
processes we assume that leptons and antileptons produced in the decays of iVi are aligned with or 
orthogonal to some specific lepton and antilepton flavours £ a and £ a (with a — e, //, r). That is, we 
assume that the only non-vanishing element in the matrix 

PaP = (la\tN 1 ){tfi\tN 1 )*. (8.1) 

that characterizes the flavour composition of the lepton £^ 1 into which N\ decays, is P aa = 1 while all 
the other diagonal and non-diagonal entries are 0. Accordingly, in the following we replace the notation 
£n ± with £ a , that is appropriate to denote pure flavour states. 

We divide the processes that generate and washout the B/3 — L a asymmetry into three classes: 

i) Ni decays and inverse decays, iVi <-> £ a </>, and on-shell iVi-mediated AL — 2 scatterings £ a <fi ^£ a cf> 
(the on-shell part of the first diagram in fig. IG.lf cO): 

ii) AL = 1 Higgs- mediated scattering processes involving the top- Yukawa coupling, £ a Ni <-* q 3 t, 
£ a q 3 *-> ATi t and l a t *-> JVi q 3 (fig. 0»). 

Hi) AL = 1 scatterings with the gauge bosons, such a i a N\ — > Acf>, £ a A — > Ni<j) and £ a (f> — > AN\ with 
A = Wi or B (figs. O^) andEHc)). 

Note that all the rates above depend (at tree level) on a single combination of neutrino Yukawa 
couplings that can be parameterized, as in eqn (|4.7|) . by rhu = rh = [AA^nv^/Mi. Other (subleading) 
washout reactions couple to lepton states that are different from £^ 1 and £n 1 - In particular, we refer 
here to the AL = 2 scatterings which go through off-shell s-channel, £(j> <-> £(f>, and t-channel, ££ <-> 4>4> 
(fig. l6.lf c0'). In the temperature regime T < Mi, and including the contributions from -/V^, the amplitude 
for these processes is proportional to the light neutrino mass matrix, [m] ai g of eqn (|2.5[) . Consequently, 
the fastest rate couples to the lepton doublet containing the heaviest light neutrino state v 3 . However, 
being of higher order in the A couplings, these processes are generally negligible, with a possible exception 
in the high temperature regimes, M\ > 10 13 GeV, where some Yukawa couplings are of order unity. Since 
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flavour effects, as well as the majority of the spectator processes, become relevant only at T < 10 13 GeV, 
in the following we consider l a and l a as the only relevant directions in flavour space. 

By definition, washout processes are lepton number violating reactions that tend to destroy any excess 
of lepton or of antileptons. In general, they depend also on the abundances of other particle species. 

The washout reactions that we consider are contained in the first three lines in eqn (|6.56[) . Note 
that since we are assuming alignment conditions, £ a = £3, some of the terms vanish. Washout reac- 
tions that involve the gauge bosons (figs. I6.1f fc), I6.1f c)) that appear in the last line in eqn (|6.56p are 
equally important. However, since the chemical potential for gauge bosons vanishes, no further density 
asymmetries arc associated with these reactions, and for simplicity we neglect them. As a result of our 
approximations and simplifications, washout rates are controlled by the following density asymmetries: 
Ayi a for the lepton doublets, Ay^ for the Higgs, Ay q3 for the third generation quark doublet and Ay t 
for the top-quark singlet. The latter two quantities always appear in the combination Ayt — Ay qa . It is 
useful to recall that Ayi a , Ay^ and Ay qa represent the sum of the asymmetries in the two components 
of the relevant S'[/(2)-doublet, and that (in the non-supersymmetric case) the densities of the lepton and 
baryon number charges are 

yal = yA La = J2( A y^+ A ve a )Y^, (8.2) 

a— e,fi,r a 

Yab = I (Ay^ + A^+A^jr^, (8.3) 

i=l,2,3 

where u,, di denote the S77(2)-smglet quarks of the i-th generation and e a are the 5't7(2)-singlet leptons. 
Similar relations hold also for the densities of the flavoured charges Ya q = Yab/3 — Yalc that are 
individually conserved by the electroweak sphalerons. 



8.2 Detailed analysis 

Relations among the asymmetry-densities Ayi a , Ay^ and Ay t — Ay q3 are determined by the chemical 
equilibrium conditions enforced by the reactions that are faster than the expansion rate of the Universe. 
The dependence of the expansion rate on temperature is different from that of various particle interaction 
rates, and as the temperature drops down, more and more interactions become faster than the expansion 
rate and 'enter into equilibrium'. Thus, the equilibrium conditions change with temperature. 

Since leptogenesis takes place at temperatures T ~ Mi, the relevant constraints depend on the value 
of Mi. As discussed in Appendix Q3J the density of the charge B/3 — L a (where in this section a is 
the flavour direction into which Ni decays) is conserved by all standard model processes, but not by 
interactions involving All the relevant asymmetry-densities can be expressed as linear functions 
of the charge asymmetry Ya = \Yab — ^al q , and the equilibrium conditions fix the coefficient of 
proportionality for each temperature range: 

Aw„ = -ce Ay 4 = - C4 ,y^. (8.4) 

The two coefficient q and encompass all the effects of the relevant spectator processes (charged lepton 
and quark Yukawa interactions, and the electroweak and QCD sphalerons). When generalized to multiple 
generations, C£ generalizes to a matrix (that is the inverse of the A-matrix [65] that is introduced at the 
end of section [9]). 

We distinguish between six relevant temperature ranges according to the set of interactions that are 
in equilibrium. For each such temperature regime we present the equilibrium conditions. We impose, 
when relevant, various conditions of flavour alignment and calculate the corresponding ce and defined 
in eqn (|8.4p . Note that ci and give a crude understanding of the impact of the respective asymmetries: 
c^/ci gives a rough estimate of the relative contribution of the Higgs to the washout, while q + 
gives a measure of the overall washout strength. The quantitative significance of the different spectator 
processes can be red off from the last column in Table [T] (adapted from ref. [204] ) where for the different 
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temperatures, and assuming a strong washout regime with m = 0.06 eV, we give the percentage variations 
with respect to the case when all spectator processes are neglected (q = 1, = 0). 

1) Only gauge and top-Yukawa interactions in equilibrium (T > 10 13 GeV). 

Since in this regime the electroweak sphalerons are out of equilibrium, no baryon asymmetry is generated 
during leptogenesis. Moreover, since the charged lepton Yukawa interactions are negligible, the lepton 
asymmetry is just in the left-handed degrees of freedom and confined in the i a doublet, yielding Yal = 
Ayt a Y cq — — Y&, a . As concerns Ay^, although initially equal asymmetries are produced by the decay of 
the heavy neutrino in the lepton and in the Higgs doublets, the Higgs asymmetry is partially transferred 
into a chiral asymmetry for the top quarks {Ay t — Ay q3 ^ 0) implying Ayi a ^ Ay^. We see from the last 
column in Table [T] that the inclusion of the Higgs asymmetry yields a sizeable reduction in the surviving 
asymmetry. 

2) Strong sphalerons in equilibrium (T ~ 10 13 GeV). 

QCD sphalerons enter equilibrium at higher temperatures than the corresponding electroweak processes 
because of their larger rate (Tqcd ~ H(a s /ow) 5 r ew [205]). These processes are likely to be in equilib- 
rium already at temperatures T s ~ 10 13 GeV [159,205,206]) and yield the constraint 

(2M 94 - - M*) = ■ (8.5) 

i 

Direct comparison with the previous case allows us to estimate the corresponding effects. The relation 
Yal = Ayi a Y ccl = — Ya q , implying eg = 1, holds also for this case. However, switching on the QCD 
sphalerons reduces the Higgs number asymmetry by a factor of 21/23. This effect yields a suppression of 
the washout that does not exceed the few percent level. 

3) Bottom- and tau-Yukawa interactions in equilibrium (10 12 GeV Si T ^ 10 13 GeV). 

The asymmetries in the SU (2)-singlet b and e T degrees of freedom are populated. The corresponding 
chemical potentials obey the equilibrium constraints 

Mfc = Mg 3 - /V' Mr = IH T ~ M0- (8-6) 

Possibly, /if, and h T Yukawa interactions enter into equilibrium at a similar temperature as the electroweak 
sphalerons [159]. However, in order to quantify separately the impact of these two effects, we first consider 
the possibility of a regime with only gauge, QCD sphalerons and the Yukawa interactions of the whole 
third family in equilibrium. As concerns the flavour composition of the lepton asymmetry, we distinguish 
two alignment cases: first, when the lepton asymmetry is produced in a direction orthogonal to l T 
(P TT = 0) and second, when it is produced in the £ T channel (P TT = 1). When P TT = 0, the lepton 
asymmetry is produced in one of the two directions orthogonal to £ T and therefore it does not 'leak' into 
the SU(2) singlet degrees of freedom, implying that q = 1 still holds. In the case where P TT = 1, the 
washout effects are somewhat suppressed, since the lepton asymmetry is partially shared with e T that 
does not contribute directly to the washout processes. Our results for these two cases suggest that the 
effect on the final value of Y& a associated to the r Yukawa interactions is of the order of 10%. 

4) Electroweak sphalerons in equilibrium (10 11 GeV T ^ 10 12 GeV). 

The electroweak sphaleron processes take place at a rate per unit volume Tew/V °t T^ctyy \og(l/aw) 
[153-155], and are expected to be in equilibrium from temperatures ~ 10 12 GeV, down to the electroweak 
scale or below [159]. Electroweak sphalerons equilibration implies 

E( 3 ^*+WJ=Q- (8-7) 

i 

As concerns lepton number, each electroweak sphaleron transition creates all the doublets of the three 
generations, implying that individual lepton flavour numbers are no longer conserved. As concerns 
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Equilibrium processes, constraints, coefficients and effects on Ya 





T (GeV) 


Equilibrium 


Constraints 


Q 




<KI^aJ) 






no spectators 




1 









> 10 13 


h u gauge 


B = Y^ i (2g I + M l + d l ) = 


1 


2 
3 


-40% 


o 


~ in 13 


■f- ^OJJ-opri 


\ ^ (Or, „, J \ n 

^iK^Qi — u i — di) — U 


1 

1 


14 
23 


— o I /o 


II 


io 12 - 13 


+ h b , h T 


T=lr-<j> 1 Ptt = 1 


1 

3 

4 


3 

8 
1 
2 


-27% 
-17% 




10 U-12 


+ EW-Sph 




98 
115 
78 
115 


41 
115 

56 
115 


-15% 
-10% 


o 














n. 

cq 


10 8-11 


+ h c , h s , hfj. 


C = g2 + t' 


151 
179 
344 
537 


37 
179 
52 
179 


-3% 
+ 12% 




< 10 8 


all Yukawas hi, h 


a Pee = 1 


442 
711 


l(i 
79 


+27% 



Table 1: The relevant quantities in the various temperature regimes. Chemical potentials are labeled here 
with the same notation used for the fields: fi qi = <?i, ni a =i a for the SU(2) doublets, fi Ui — Ui, fj,^ —di, 
/i ei = ei for the singlets and = <p f° r the Higgs. The relevant reactions in equilibrium in each regime 
are given in the second column and the constraints imposed in the third. The conditions adopted for P aa 
are indicated (the appropriate constraints on the conserved quantities = B/3 — Lp with j3 = e, fi, t 
should also be imposed) . The values of the coefficients eg and are given in, respectively, the fourth and 
fifth column. For each regime (and assuming to — 0.06 eV in all the cases), the last column quantifies the 
percentage variation in the absolute value of Ya q with respect to the case when all spectator processes 
are neglected (first line in the table). 



baryon number, electroweak sphalerons are the only source of B violation, implying that baryon number 
is equally distributed among the three families of quarks. In particular, for the third generation, B3 = B/3 
is distributed between the doublets 93 and the singlets t and b. In Table Q] we give the coefficients ct and 
for the two aligned cases: (i) P TT — implying Ya t = 0, and (ii) P TT = 1 implying Ya c = Ya (1 = 0. 
We see that in this case the transfer of part of the lepton asymmetry to a single right-handed lepton (e r ) 
can have a 5% enhancing effect on the final value of Ya q • 

5) Second generation Yukawa interactions in equilibrium (10 8 GeV T 10 11 GeV). 

In this regime, the h c , h s and interactions enter into equilibrium. We consider two cases of alignment: 
(i) P ee = 1 implying Y& — Ya t = 0, and (ii) P ee = 0. To ensure a pure states regime we further assume 
complete alignment with one of the two flavours with Yukawa interactions in equilibrium, for defmitcness 
P TT = 1, and therefore Ya c = Ya„ = 0. The difference in eg between the two aligned cases is larger than 
in the regimes 3 and 4, and accordingly the difference in the corresponding values of 1a q , of the order 
of 15%, is somewhat larger than in the cases in which just the third generation Yukawa couplings are in 
equilibrium. 

6) All Yukawa interactions in equilibrium (T <, 10 8 GeV). 

In this regime, since all quark Yukawa interactions are in equilibrium (actually this only happens for 
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T < 10 6 GeV), the QCD sphaleron condition becomes redundant. Hence ignoring the constraint of eqn 
(|8.5p . as is usually done in the literature, becomes fully justified only within this regime. If, however, 
leptogenesis takes place at T > 10 s GeV, as occurs for hierarchical singlet TV's, the constraint implied by 
the QCD sphalerons is non-trivial, even if the associated numerical effects are not large. 

Due to the symmetric situation of having all Yukawa interactions in equilibrium we have just one 
possible flavour alignment (the other two possibilities being trivially equivalent). We take for dcfinitcncss 
P ee = 1, implying = Y"a t = 0. In this case eg is reduced by a factor of almost two with respect to the 
case in which the spectator processes are neglected (eg = 1) and the final value of Fa q is correspondingly 
enhanced. The reason for the reduction in ci can be traced mainly to the fact that a sizable amount 
of B asymmetry is being built up at the expense of the L asymmetry, and also a large fraction of the 
asymmetry is being transferred to the right-handed degrees of freedom at the same time when inverse 
decays and washout processes are active, reducing the effective value of Aj/£ c that contributes to drive 
these processes. 

To summarize, we considered the possible impact of the spectator processes. Conditions of flavour 
alignment /orthogonality were imposed, to ensure that these effects are cleanly disentangled from lepton 
flavour effects (these are numerically more significant, but of a different nature). A rough quantitative 
understanding of spectator processes can be obtained by relying on the fact that the surviving asymmetry 
is inversely proportional to the washout rate, as discussed in section 21 Hence, the final Ya asymmetries 
obtained in the relevant temperature regimes will be inversely proportional to q + c^. This suggests 
that numerical corrections related to the proper inclusion of spectator processes have at most 0(1) 
effects [204]. Inspecting the table, we learn that when the electroweak sphalerons are not active and all 
Yukawa interactions (except those of the top-quark) are negligible, the Higgs contribution enhances the 
washout processes, leading to a smaller final asymmetry. As more and more spectator processes 
become fast (compared to the expansion rate of the Universe), the general trend is towards reducing the 
value of the washout coefficients and hence increasing the Ya q up to values that can be slightly larger 
then what is obtained when all spectator processes are neglected. 
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9 Flavour Effects 



The aim of this section is to discuss what are flavour effects in leptogenesis, why they arise, and when 
they matter. It makes use of results from Sections HI [5] and Appendices H31 HH and [15l but should be 
independently readable. In Section I^TTl we introduce the puzzle of flavour in leptogenesis calculations, 
which is that the baryon asymmetry depends on the choice of lepton flavour basis. We also give heuristic 
rules for how to choose the correct basis. Sections 19 . 21 and \9 .31 illustrate how flavour effects can modify the 
final baryon asymmetry. Spectator processes (including sphalerons) are neglected in Section l9~2l where we 
focus on how flavour dynamics can strongly enhance the asymmetries in lepton doublets Y A e a ■ The sum 
of the Y A £ a is usually of the order of the final baryon asymmetry; this is discussed in Section 19.31 where 
spectator effects are included, and the relation between the charge densities Y Aa and the asymmetries 
for the various particle species is elucidated. Notice that it is important to distinguish the asymmetries 
in lepton doublets Y A g a from the B/3 — L a charge densities Y Aa , even though they are similar in our 
notation. In Section [9.41 we analyze the main phenomcnological consequences of including the effects of 
the lepton flavours. 

Historically, leptogenesis calculations were performed in the "single flavour approximation" , which 
consists of studying Boltzmann Equations for the B ~ L asymmetry. Reference [65] considered the BE 
for the asymmetries in B/3 — L a , but did not emphasize that the results were significantly different from 
the single flavour approximation. In subsequent years, various authors [173,207,208] noticed that flavour 
effects could be used to enhance the baryon asymmetry in particular models. The flavoured BE are given 
in [66]. The importance of studying lepton asymmetries flavour by flavour, in order to obtain a "reliable" 
estimate of the baryon asymmetry, was recently highlighted in [61-64]. 



9.1 The flavour puzzle 

In the introduction to Section [4j washout interactions were presented as a critical ingredient for thermal 
leptogenesis, and the importance of flavour effects was said to follow from the importance of washout: 
the washout interactions have lepton doublets in the initial state, so to compute the washout rates one 
needs to know which leptons are distinguishable. In Section[B]we became acquainted with the Boltzmann 
Equations, which describe the detailed evolution of asymmetry production and of washout processes. 
Upon solving the BE, one could expect to get a correct description of washout, and discover if flavour 
matters in leptogenesis. Unfortunately, this is not the case. The BE are not covariant in lepton flavour 
space, and solving them in different bases gives different answers. 

A simple example can illustrate this problem. Suppose that Ni decays to £^ and £ T with equal 
branching ratios, so that £m x = (£ /1 + £ T ) / 'V2 (see eqn ([5. lip ). Furthermore, assume that the asymmetries 
in both flavours are the same: e w = e TT . Let us consider the simple BE eqn (|6.30p . and neglect the 
Higgs asymmetry Ay^ and all spectator processes, that are irrelevant for the present discussion. In this 
approximation the total lepton asymmetry coincides with the asymmetry stored in the lepton doublets: 
Yal = Y cq J2 a ^yi a = Y oq Aye. To write down the BE for the total lepton asymmetry we can proceed in 
two ways. We can chose to work in the basis of £n ± and of two other flavours orthogonal to £^ x that, by 
assumption, are not produced in N± decays. In this case the only CP asymmetry is e^ N = J2 a 6aa = e ' 
and furthermore rhf = fh. The BE for the total lepton asymmetry then is 
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(9.1) 



Alternatively, we can chose to work in the "flavour" (charged lepton mass eigenstate) basis, and in this 
case we shall write two BE for the jj, and r asymmetries. By summing them up we obtain 



Y AL = Y A1 
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(9.2) 
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where we have used m MM /m = rh TT /m = 1/2. Eqns (|9.ip and (|9.2p differ by a factor of 1/2 in the 
washout, and therefore they yield different baryon asymmetries. This simple example illustrates what 
was (cryptically) mentioned at the end of section [U the rough estimates for the baryon asymmetry given 
in eqns (|4. 15|) and (|4.16[) depended on the choice of the flavour basis. 

To obtain a reliable result, we either need a formalism that is flavour covariant, or we need to guess 
which is the correct basis. Here we opt to guess, based on physical intuition. (A flavour covariant toy 
model, that motivates our guesses, is discussed in Appendix [151) 

Two basic points provide guidance in making the correct guess: 

1. washout processes are critical in leptogenesis; 

2. interactions whose timescale is very different from that of leptogenesis, drop out of the BE. 

The second point is usual Effective Field Theory: interactions which are strong should be resummed, very 
weak interactions can be neglected. In particular, interactions that are much faster than the timescale of 
leptogenesis and of the Universe expansion rate, are "resummed" into thermal corrections (see Section [7]) 
and into chemical equilibrium conditions (see Section [5]). The second point relates to the first, because 
to perform a correct calculation of the washout rates, one should make the correct choice of basis for 
the lepton doublets which are in the initial state. The fast flavour-dependent interactions, mentioned in 
point 2. are precisely the ones that can resolve the flavour basis ambiguity. 

The lepton doublets of different flavours are distinguished in the Lagrangian by their Yukawa cou- 
plings h a . During leptogenesis, they will also be distinguishable, if the h a mediated interactions are 
fast compared to those of leptogenesis and to the Universe expansion rate. Since the h e ^ T mediated 
interactions are of widely different strengths, they would, for instance, induce differences in the thermal 
masses of the different leptons. The answer to the "flavour puzzle" in the previous example is that when 
(some of) the charged lepton Yukawa interactions are "fast enough", the "flavour" (charged lepton mass 
eigenstate) basis is the correct basis for the BE, and eqn. ()9.2|) should be used. In the opposite situ- 
ation in which charged lepton Yukawa interactions are much slower than the Universe expansion rate, 
leptogenesis has no knowledge of lepton flavours, and the correct BE is eqn (|9.ip . 

To be more quantitative, let us estimate the temperature below which lepton flavor effects cannot be 
neglected. The interaction rate for a charged lepton Yukawa coupling h a [195,209] can be estimated as 

r n ~5x l(T 3 ^T, (9.3) 

(for details see Appendix [T4l around eqn (|14.10|) ). The condition r T ^) H then implies that the rate for 
h r {a) mediated processes becomes faster than the expansion rate of the Universe below T ~ 10 12 GeV 
(10 9 GeV), 20 while for T > 10 12 GeV the charged lepton Yukawa coupling s are irrelevant. Hence, above 
this temperature all flavour effects can be neglected, and the only "basis-choosing" interactions are those 
involving the neutrino Yukawa couplings A. 21 

If leptogenesis occurs below T ~ 10 12 GeV, the /i r -related reactions participate, along with the neu- 
trino Yukawa, in determining a physical flavour basis (via, for instance, the "thermal mass matrix"). 
The conditions under which i T is singled out as a distinct lepton are discussed in Section ["15.41 see e.g. 
eqn (|15.37p . For instance, if at the time the lepton asymmetry starts to survive, the h T interactions are 
also faster than the Ni inverse decays [180,181], then the relevant basis states are t T and the component 
in £ n 1 that is orthogonal to £ T , and the total lepton asymmetry is shared between these two flavours. 
This situation can have significant consequences. In particular, it can happen that the asymmetry in the 
r flavour is larger in size than (and even opposite in sign to) the total lepton asymmetry generated in 
Ni decays (see eqn (15.1010 . If, in addition, the r flavour is weakly washed-out, a sizeable fraction of its 
asymmetry can survive until the end of the leptogenesis era. An example of two cases in which flavour 
effects yield a large enhancement of the baryon asymmetry is given in figure 19.11 

20 In SUSY h T = m T /(v u cos/3), so the tau Yukawa is in equilibrium for T < (1 + tan 2 0) X 10 12 GeV. 
21 AL = 2 interactions mediated by Ni exchange could also become fast at large temperatures. Besides giving additional 
contributions to the washout, this also complicates the issue of the physical basis for the lepton doublets. See e.g. [65]. 
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Figure 9.1: The total baryon asymmetry in the two flavour calculation (upper curves) and within the one-flavour 
approximation (lower curves) as a function of z, for two different sets of washout parameters. Left picture: 
m TT /m, = 10, fhfifjt/m* — 30, m ee /m» = 30. Right picture: m TT /m* = 10, m,j M /m» = 30, m ee /m» = 10~ 2 . In 
both cases e TT = 2.5 x 10~ 6 , e MM = -2 x 10" 6 , e ee = 10" 7 and Mi = 10 10 GeV. These plots are from [64]. 



9.2 Enhancement of the B — L asymmetry 

To simplify the analysis, we focus in this Section on the asymmetries in lepton doublets, and neglect 
spectator and sphaleron processes. We make two assumptions to fix the flavour basis: 

(i) The rate of the /i r -interactions T T is much larger than H and Tjjj, the rates of, respectively, the 
expansion of the Universe, and the Ni inverse decay; 

(ii) During the entire period of leptogenesis, the rate of the /i^-interactions, T^, is either <C H,Tid, or 
^> H,Tj£), so that the flavour basis does not change during leptogenesis. 

When these two assumptions hold, it is a good approximation to work just with the projections of the £n 1 
densities onto (for T M <C H,Tjd) the two-flavour basis (£ , £ T ), where £ is the direction of the component 
in In 1 that is orthogonal to £ T , or (for T^ ^> H,Ti£>) the three flavour basis (£ e ,£^,£ T ). 

The flavour projectors correspond to the diagonal entries of the matrix P in eqn l|8.ip (and P for the 
CP conjugate states) [65], that is 

r(N 1 ^t a <i>) = r(N 1 ^t Nl <t>)p aa . (9.4) 

The CP violating differences, 

^Paa Pact Pact: (^-^) 

are important quantities that account for the misalignment in flavour space of the states and £tf 1 . 
At tree level, AP QQ = and the flavoured washout parameters rh aa can be simply written as 

m aa = mP aa . (9.6) 

At the loop level, however, the AP aa do not vanish in general. We use eqns (|9.4I9.5I9.6[) to rewrite the 
asymmetry in lepton flavour a, ean. (|4.2[ ). as 

e aa = e^ + ]-AP aa . (9.7) 
m I 

The crucial observation is that, while washouts are dominated by the inverse-decay rates, which are 0(X 2 ) 
and do not involve the couplings to N2 3, the CP-violating difference T(Ni — > £ a (f>) — T(A r 1 — > £ a (j)) is 
C(A 4 ) and does involve also the couplings Ap^ with k = 2,3. Therefore, as shown by eqn (|9.7p . in general 
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e aa is not proportional to m aa . This implies that, for fixed values of e and fh, one can have e aa large 
together with fh aa small, and this can yield a strong enhancement of the asymmetry Y A i a . Taking into 
account that ^ Q P aa = ^ a P aa = 1, we obtain the following 'sum rules': 



m aa = fh, ^e aQ =e. (9.8) 

a a 

The first term in eqn (|9.7p corresponds simply to the projection of the total asymmetry e onto the flavour 
a. The second contribution plays a more subtle role. It represents a type of CP-violation that is specific 
of the flavour CP-asymmetries and does not affect the value of e. 
The Boltzmann equation for Y/v is given in eqn (|6.51[) : 

Y N = -(y Nl - 1)[ 7 a^2 + 7t 2 o7 2 + t3T 2 ]- (9-9) 

As concerns the Boltzmann Equations (I6.57|) for the Ya„, we remarked at the end of Section [6] that, to 
write them in closed form, one has to include all the spectator processes. This induces a mixing between 
Ya q and Y Af} that complicates the expressions of the BE. For the sake of clarity, we provisionally neglect 
these effects (we discuss them in Section l9.3p . This approximation allows us to consider, instead of 
coupled equations for the Y Aa 's, independent equations for each of the lepton doublet densities Y A £ a 's: 

Y Aia ~ -Y N e aa + (Y Ae J w . (9.10) 

Note that since we neglect the redistribution of the asymmetry to singlet leptons, we took (Y A i a ) w = 
— (YA a ) w when using eqn (|6.56p . The solutions are of the form of eqn (|4.15[) : 

Y A t = - 4 x 10~ 3 Y t aa Va{m, rh aa ). (9.11) 

a a 

where the summation is over a — e, fj,, t for T < 10 9 GeV, and a = o,t for 10 9 GeV < T < 10 12 GeV. 
Approximate analytic solutions (reviewed in Appendix ITS}) to eqns (|9.9[) and (|9.10p give 



( m aa \ 
\ 2.1m, J 



2m* m. 



fh > m* (g 12) 

m < m*. 



Eqns (|9.1ip and (|9 . 1 2[) allow us to estimate the lepton doublet asymmetry Y A g in various washout 
regimes. In the case that all the flavours are in the strong or mildly strong washout regime, rh aa > to* for 
all a, the effects of the AP term in eqn (|9.7p are most striking. The efficiency for Y A £ a is approximately 
given by rj a ~ m*/(2m aa ). Inserting eqn (|9.7p in eqn (|9.1ip . we obtain 

Ns ^ € + y^^. (9.13) 
2 x 10- 3 1 fh ^ fhaa 2 v ; 

a 

Here Aff denotes the number of lepton flavours that are effectively resolved by the fast charged lepton 
Yukawa interactions. The first term represents an enhancement by a factor A// with respect to the 
doublet asymmetry that would be obtained by neglecting flavour effects. It is clearly independent of 
particular flavour structures and it would be the only effect in the special cases where e aa °c fh aa - Its 
origin is simple to understand. On one hand, the asymmetry in each flavour is reduced by a factor of the 
branching ratio, TO aQ ,/TO. On the other hand, the efficiency factor is increased by the same factor, because 
all the washout rates for Y A £ a are reduced by it. The two factors cancel each other, and an asymmetry 
Y A i a tx £to*/to is generated in each flavour. The sum over flavours yields the Aff enhancement. 

In contrast, the second term depends crucially on the details of the flavour structure. Let us analyze 
its possible effects in the Aff = 2 case, in which AP TT = —AP oa (with AP oa = AP ee + AP W ). If we 
assume to S> to*, then the first term in the r.h.s of eqn (|9.13p is strongly suppressed. However, flavour 
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Figure 9.2: The absolute value of the final |1a(b-l)| (i n units of 10 _5 |e|) as a function of fh TT . We 
use Mi = 10 10 GeV and rh = 0.06 eV. The dashed line corresponds to the case e aa oc fh aa for which 
AP = 0. The solid line gives an example of the results for AP ^ 0. We take AP TT /e oc J fh TT / fh. The 
corresponding values of e TT /e are marked on the upper x-axis. The arrows with labels (a), (b), (c) and 
(d) correspond to the four situations discussed in the text. Note that around (c) Ya(b-l) changes sign. 
(Figure adapted from ref. [63]). 



configurations are possible for which fh TT ~ to* (or fh 00 ~ m*). Then ,the production of an asymmetry 
Ya£ t (or 1a^„) driven by the second term occurs under 'optimal' conditions. This qualitative expectation 
is confirmed by numerical integration of the flavour dependent Boltzmann equations. Fig. 19.21 (adapted 
from ref. [63]) depicts a set of possible cases. The absolute value of the final |Ya(b-zoI ( m units of 10~ 5 |e|) 
computed with Mi = 10 10 GeV and fh = 0.06 eV is shown as a function of fh TT . Two flavour-aligned cases 
at rh TT = (i.e. no r component in £zvi) and at rh TT — rh (i.e. tj^ 1 = £ T ) coincide with the results that 
would have been obtained by neglecting flavour. For rh TT /rh ^ 0, 1 quite different results are obtained. 
Case (a) (dashed line) gives the results obtained by setting AP = 0, that is by taking e aa oc fh aa . As 
predicted in our qualitative analysis, the final B — L asymmetry is practically independent of the flavour 
structure, and gets enhanced by a factor ~ TV/ = 2 with respect to the flavour-aligned results. 

For the more general cases where AP ^ 0, the final doublet lepton asymmetry strongly depends 
on the particular flavour structure. To depict in a simple way the various possibilities, we adopt, as a 
convenient ansatz, the relation AP Tr oc */ rh TT / fh. This is based on the fact that, if we take the tree-level 
Ni — * t r (j) decay amplitude to zero (A T i — > 0) while keeping the total decay rate fixed [(A^A)n = const.], 
r~h TT vanishes as the square of this amplitude (oc |A r i| 2 ) while AP TT vanishes as the amplitude (oc A T i). 
We extrapolate this proportionality to finite values of rh TT because this allows to represent the AP TT 
effects in a simple two-dimensional plot of Ya(b-l) versus rh TT . 

On the upper x-axis of Fig. [9"T2"l we mark for reference the relative value of the r CP-asymmetry e rT /e 
corresponding to the different values of rh TT . The most peculiar features are the two narrow regions 
marked with (b) and (d) where |1a(b-l)| is strongly enhanced. In (b) this happens around the 'optimal' 
value rh TT ~ m* that yields a strong enhancement, even though e TT /e is rather small. The region marked 
with (d) corresponds to rh 00 ~ and to e oa of the opposite sign with respect to e. Here the steep rise 
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of |AYa(b~l)I can reach values up to one order of magnitude larger than the vertical scale of the figure. 
(Note, however, that the ansatz AP TT oc */ rh TT /rh does not yield the required behavior AP — at the 
boundary fh TT /fh = 1 and therefore the continuous line in the plot should not be extrapolated too close 
to this limit.) It is worth noticing that the two peaks correspond to values of Va(b-l)A of opposite sign. 
In fact, the B — L asymmetry changes sign around point (c), and for rh TT /rh <; 0.85 it is of the opposite 
sign with respect to what one would obtain neglecting flavour effects. Similar results are expected in the 
temperature regime below T ~ 10 9 GeV, when both the r and the /j, Yukawa reactions are in equilibrium. 
In this case even larger enhancements are possible because the dynamical production of the asymmetry 
of two flavours can occur in an 'optimal' regime. 

In the weak washout regime for all the flavours, there can be no dynamical enhancement of the 
asymmetry and flavour effects are less important. The efficiencies for this regime are given in eqn (|4.13|) 
(or (|9. 12[) V Using also eqn (|9.7p , we obtain 
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The first sum within square brackets is always < 1 so this contribution by itself is suppressed with 
respect to the unflavoured case. The second sum is always less than or equal to the largest of the AP aa 's. 
Therefore, in this regime a necessary condition to have a large enhancement of the resulting doublet 
lcpton asymmetry is that AP aa 3> e for at least one lepton flavour. 

Finally, in the case where one flavour (3 is weakly coupled to Nx, mpp/m* <g; 1, but other flavours 
a are strongly coupled, rh aa /m* ~ m/m* ^ 1, the singlet neutrinos are brought into equilibrium by 
the reactions involving the strongly coupled flavours, so that hn 1 ~ n 1 . The efficiency factor rip 'loses' 
the factor of m/m*, which corresponds to a suppressed N\ abundance [see eq. (14. 13(1 ] . Consequently, the 
asymmetry for the weakly coupled flavour (3 is given by 
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The asymmetry for the strongly coupled a flavours is again given by eqn (|9 . 1 3[) with A// = 1 or 2. 



9.3 Mixing of lepton flavour dynamics 

In Section [6] we derive BE for the charge densities Y Aa , without taking into account the SM spectator 
processes. For this reason, the dependence of the Ay a 's on the Y& a 's is left implicit. In Section [H] we 
analyze the spectator processes, assuming that TVi decays to pure (flavour) states, that is, neglecting 
flavour effects. In this section we combine these two issues, and we outline how to express the washout 
term in eqn (|6.56p in terms of the Y Aa (see Appendix [T4l for a more complete treatment). 

We assume that asymmetries can be generated in two flavour directions, which we take to be f and 6, 
where the latter is some linear combination of fi and e. Below T ~ 10 12 GeV, baryon and lepton numbers 
violating electroweak sphaleron processes, and /i r -related processes are fast. In this situation, the charge 
densities that are (slowly) evolving because of the leptogenesis processes are 

Y Aa = ]-Y AB - Y ALa = ~ ( A Vv + A Vu* + &yd,)Y c * - {Ay ta + Ay e jY^. (9.16) 

i=l,2,3 

Note that, in the temperature range of interest, all the asymmetry-densities in the right-hand-side of 
this equation are generally non- vanishing. In particular, Ni decays generate a non-vanishing chemical 
potential for the Higgs particles which, in turn (see eqn (|14.5p ). induces a non- vanishing chemical potential 
to the top-quarks. The QCD sphaleron condition of eqn (|8 . 5[) further implies that this happens for all 
the lighter quarks. The only exceptions are Ay e above T ~ 10 9 GeV, and Ay Cc that is not populated 
during standard leptogenesis. 
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Figure 9.3: The evolution of the charges densities Y& T (green dashed line), Y"a (blue dotted line) and 
Ya(B-l) ( re d solid line) in units of e as a function of z — M\/T. We use M\ = 10 10 GeV and 
fhx = 0.06 eV. The b- and r- Yukawa reactions and the electroweak sphalerons are in equilibrium. The 
figure refers to the point labeled (d) in fig. 19.21 that corresponds to AP TT — 0.98 and e TT w 1.2 e. (Figure 
adapted from ref. [63]). 



The asymmetry densities for the quarks and for the SU (2)-singlet leptons can always be re-expressed 
in terms of the asymmetry densities of the S'f7(2)-doublet leptons. This procedure makes use of the 
condition of hypercharge neutrality, eqn (|14.4p . that involves all three generation fcrmions. Consequently, 
once the appropriate substitutions are implemented, one obtains that the densities of the A Q charges in 
cqn (|9.16[) correspond to linear combinations of the asymmetry densities of all the lepton doublets: 

Y Aa = - [A]~l Ay^ Y*\ (9.17) 

The matrix [A] -1 can then be inverted to give Ayi a in terms of the charge densities Y"a q (see Appendix [14"]) . 
The same is true also for the Higgs asymmetry. We can thus write: 

(3 

The matrix A (introduced in ref. [65]) and the vector (introduced in ref. [204]) generalize the coeffi- 
cients ct and c^, defined in eqn (|8.4[) for the case of pure flavour states, to the general flavour case. Their 
numerical values are determined by the relevant set of chemical equilibrium conditions, and therefore 
depend on temperature. 

As a consequence of eqns (|9.18p . the evolution of one density Ya q depends on the other charge 
densities, that is, the BE for the lepton flavours form a system of coupled equations. This phenomenon 
is sometimes referred to as 'electroweak sphaleron flavour mixing', and was studied in refs. [63,210]. 
One example of the possible effects is depicted in fig. 19.31 The plot gives the evolution of Ya t , Ya 
and Ya(b-l) = ^A r + ^Ao f° r the flavour configuration corresponding to the point labeled with (ef) in 
fig- 19.21 Ya t is in the strong washout regime, rh TT w 0.06 eV, and the r CP asymmetry is relatively 
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large, e rr ~ 1.2 e. For Ya oJ the value to do « 10 _3 eV is close to optimal, and the asymmetry is smaller in 
size and opposite in sign to e TT : e ao « — 0.2e. For the relevant temperature, T 10 12 GeV, electroweak 
sphaleron, and /i r reactions are in equilibrium, yielding the matrices of coefficients [in the (Va oJ ^a t ) 
basis] : 

"-as 03 ™)' c * = t>. 56 >- < 9 - w > 

As apparent from fig. 19.31 Ya d gives the dominant contribution to the final B — L asymmetry A 
qualitative explanation goes as follows. For z J> 0.1, a large asymmetry builds up quickly in the r flavour 
and contributes to the early washout of Ya (the sign of Ya t that is opposite to that of Ya combines with 
the negative sign of A OT = —24/230 yielding an increase of the washout rates). Thus, the early (negative) 
asymmetry in Ya„ remains suppressed. It follows that the opposite sign o asymmetry generated from N\ 
decays at z <; 1 is largely unbalanced. As concerns the r asymmetry, typically of strong washout regimes, 
it 'freezes-out' at z > 1 with an absolute value that is sizably smaller than the values reached at earlier 
times. 

Remarkably, as can be clearly seen in fig. 19. 3[ the B — L asymmetry changes sign twice during its 
evolution. This peculiar behavior is qualitatively different from what would be obtained by neglecting 
the sphaleron induced flavour mixing, when there is a single change of sign. Other interesting washout 
effects induced by sphaleron mixing have been analyzed in ref. [211]. 



9.4 Phenomenological consequences 

As discussed in Section 19.21 starting from a specific set of parameters within a given seesaw model, the 
inclusion of flavour effects in the calculation of the baryon asymmetry can change the result by factors 
of a few to orders of magnitude [63,64]. In this section, we take a more "bottom-up" perspective of 
flavoured leptogenesis. As discussed below eqn (|2.5p . there are (currently) 14 unknown parameters in the 
three generation type I seesaw model, of which only a few are accessible to low energy experiments. It 
is interesting, from a phenomenological perspective, to constrain the unknown parameters by requiring 
that leptogenesis works successfully. We will find that the inclusion of flavour effects does not change the 
resulting constraints in a significant way. 



9.4.1 The upper bound on the light neutrino mass scale 

In the "single flavour" calculation, successful thermal leptogenesis requires a light v mass scale ^ 0.2 
eV [135,212]. The argument leading to this bound does not apply in the flavoured calculation. Currently, 
there is no consensus in the literature on the value of the bound in the flavoured case. Analytic arguments 
(which we reproduce below) and the numerical analysis of Refs [62,180,210], suggest that flavoured 
leptogenesis can be tuned to work for m mm up to 0.5 — 1 eV. This is of the order of current bounds 
from Large Scale Structure data and WMAP [119,177-179]. The leptogenesis bound is relaxed because 
there is more CP violation available in flavoured leptogenesis. The upper limit of eqn (|5.17p on the total 
CP asymmetry decreases like Am^ tm /m max , as the light neutrino mass scale m max increases. There is 
therefore an upper bound on m max (~ m min for degenerate neutrinos). However, the limit (|5.17|) does 
not apply to the flavoured CP asymmetries, which can increase with the light neutrino mass scale, as 
shown in eqn (|5.10j) . 

The limit m v 0.2 eV [135], obtained in the single flavour approximation, can be understood as 
follows. For such a high mass scale, the light neutrinos are almost degenerate, with TOi ~ — m^. 
Due to the lower bound rh\ > m m i n , leptogenesis takes place in the strong washout regime. Taking the 
total CP asymmetry to be close to the upper bound of eqn (|5.17|) . the final baryon asymmetry can be 
approximated as 

3 to* to, Aml tm 
Yab ~ 10 e— oc - ^ . (9.20) 

As the light neutrino mass scale is increased, M\ and, correspondingly, the temperature of leptogenesis, 
must increase to compensate the ATO atm /TO^ in suppression. The leptogenesis temperature is, however, 
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bounded from above, by the requirement of having lepton number violating processes out of equilibrium. 
For instance, the requirement that the AL = 2 processes are out of equilibrium (we estimate T ~ 7/n eq , 
with 7 from eqn 113.2!))) gives 

TT^V^ M p l 

which implies M\ ^ 10 10 (eV/m m i n ) 2 GeV. There is therefore an upper bound on the baryon asymmetry 
which scales as l/mj 1 j n . A large enough baryon asymmetry is obtained for m m ; n ^ 0.1 eV. 

Consider now the flavoured case. We define c 2 = P TT and s 2 = P ao (with s 2 + c 2 = 1), so that r rr = 
c 2 r d and r oo = s 2 T]j. The individual flavour asymmetries, e aa , satisfy the bound of eqn (|5.10p . while 
their sum satisfies the stricter bound of eqn (|5.17p . So, if we can arrange — e 00 = e TT ~ 3Afim m i n / (167rw 2 ), 
then for strong flavoured washout, and for sufficiently slow AL = 2 processes, the final baryon asymmetry 
(from e.g. eqn (|4.15p ) is 

, 3M x m* c 4 - s 4 1 1 c 4 - s 4 Mi 

^ ^ 10 ~ 3 X 10 ^lO^eV- 

To show that the m m i n ^ 0.2 eV bound does not apply in the flavoured case, we must demonstrate that 
the three ingredients that entered the derivation of eqn (|9.22j) - (i) an appropriate set of flavoured CP 
asymmetries, (ii) negligible AL — 2 processes, and (iii) strong flavoured washout - can be realized in a 
model. 

(i) To obtain the desired decay rates and asymmetries we use the Casas-Ibarra [101] parametrization, 
in terms of rrn, Mj, U and the complex orthogonal matrix R = vD rr ^^ 2 U T XD 1 J^ 2 (see the discussion 
around eqn (|2.12|) ). This ensures that we obtain the correct low-energy parameters. We take o = /i for 
simplicity, and write R as a rotation through the complex angle 9 — <p + irj. The flavour asymmetries, 
for maximal atmospheric mixing, are 

3Mi Im{U* k RkiU w Rpi} 3Mim min sinh ?7 cosh n 
- ■ cos2<^ , 



16™ 2 EalEi^-^jil 2 16 ™ 2 sinh 2 t] + cosh 2 77 

iM 1 m mia lm{R kl U* k R p iU Tp } ^ 3Mi7n min sinhTycoshTy 

16^« 2 EalEj-Rjl^ajl 2 ~ + 167TW 2 sinh 2 7] + COSh 2 77 



The flavour dependent decay rates and T TT are proportional to 



cos 2ip 



(9.23) 



2 Mim min 2 iWitrw, . ,2 . ,2 . • s 

= — T~2 — -R11 + -R21 = — — (smh 77 + cosh ?7 + sm2^) 
2ir 2ir 

._, _ Mim mill , 2 _-^l m min,. ,2 .2 



| A ]a = o 2 |ii u _ | 2 = L - 2 mm (sinh^ + cosh 2 r? - sin 2y>). (9.24) 
2tr 2t) 2 

We see that we can choose 77 and </? as required to obtain eqn (|9.22[) with c 4 — s 4 <; s 2 c 2 . 

(ii) Eqn (|9.2ip gives the condition for the AL = 2 processes to be out of equilibrium: Mi & 
10 10 (eV/m min ) 2 GeV. Combined with eqn (f9T22|) . this gives 

Y B ~ 3 x 10- 11 C -^4- (—) 2 ■ ( 9 - 25 ) 

(iii) In this strongly washed out area of parameter space, N± decays take place in the flavoured regime 
if T T > Fid [180,181] (see Appendix |15.4[) . Suppose s 2 < c 2 . Then eqn (|9.22|) is a valid estimate, 
provided T T ^> T/u when s 2 Tie> = H: 

T T 10 12 GeV 1 1 

— ,9 - 26) 

We conclude that with careful tuning, thermal leptogenesis can work for m m ; n Si eV. This can be seen 
in Fig. 19.41 For the plot on the left, flavour effects are included, and leptogenesis can work for values 
10~ 4 eV ^ 777! <, 1 eV. At the upper bound we have rh 1 ~mi, so the cosmological bound is saturated. 
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Figure 9.4: Points in the Mi — rhi plane where thermal leptogenesis generates a sufficient baryon asym- 
metry (updated from [210]). These plots clearly show the lower bound on M 1 arising when the Ni are 
hierarchical. The plot on the right corresponds to a single flavour calculation. The plot on the left 
includes flavour effects. The lower bound on Mi is similar in both cases. Leptogenesis works for larger 
values of rhi in the flavoured plot, which is obtained with two-flavour BE, and imposing the constraint 



9.4.2 Sensitivity of Yab to low energy phases 

An important, but disappointing, feature of "single-flavour" leptogenesis is the lack of a model- independent 
relation between CP violation in the leptogenesis processes and the (observable) phases of the lepton mix- 
ing matrix U (see eqn (|2.7[1 ). With three singlet neutrinos N, thermal leptogenesis can work with no 
CP violation in U (see e.g. [213]) and, conversely, leptogenesis can fail in spite of non-vanishing phases 
in U [214]. In some specific models, however, it is possible to establish a relation [214,215]. 

In unflavoured leptogenesis, it is simple to see that the baryon asymmetry is insensitive to phases of 
U [214]. Using the Casas- Ibarra parametrization of eqn (|2.11[) in eqn (|5.13p , e can be written in a form 
where U does not appear: 

6 = E^W Mk 9 \M) ■ (9 ' 27) 

Alternatively, in the top-down parametrization in terms of Vr, Vl and D\ (see Section [2 .1.1 1) , e depends 
on the phases in Vr, while the U phases could arise only from Vl- 

The latter observation can be exploited to relate the U phases to flavoured leptogenesis in models 
where there is no CP violation in Vr 22 [63,64,217-220]. However, for flavoured leptogenesis with three 
iVj's, it is still true that the baryon asymmetry is not sensitive to phases of U [221], in the sense that the 
baryon asymmetry can be accounted for with the phases of U having any value. Conversely, if the phases 
of U are measured, the baryon asymmetry is still not constrained. This can be seen from the scatter 
plots of ref. [221], which show that, for any value of the U phases, the baryon asymmetry can be large 
enough. 

In the two Ni model (see section 12.1.21) . there are three phases. Relations between U phases and 
leptogenesis were obtained in unflavoured leptogenesis [106]. The relation between the flavoured CP 
asymmetries e aa and the phases of U, focusing on various texture models, was discussed in [222]. Possible 

22 Such models may be difficult to construct [216]. 
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relations between low energy observables and the baryon asymmetry, in the general case and for particular 
texture models, were discussed in [64]. 

In minimal left-right models with spontaneous CP violation, there is a single phase which controls all 
CP violation including that for leptogenesis [223,224]. 

9.4.3 The lower bound on T re heat 

In the space of leptogenesis parameters, there is an envelope inside which leptogenesis can work. In the 
single-flavour calculation, the most important parameters are Mi, (or, equivalently to), e and the light 
neutrino mass scale [212]. Including flavour gives the parameter space more dimensions (Mi, e aa , r QQ ...), 
but the envelope can still be projected onto the M\ — fh plane by fixing e to its maximum value of eqn 
(I5.17p . which is a function of M\. 

For thermal leptogenesis to take place, an adequate number density of N\ must be produced by 
thermal scattering in the plasma, suggesting T rc hoat *S Mi/5 [59,60]. There is a lower bound on Mi, from 
requiring a large enough CP asymmetry. The corresponding bound on T re h ea t depends on the fine details 
of reheating, the dynamics of which has been studied in refs. [59,60]. To avoid these complications, the 
lower bound on Mi is usually quoted. 

For non-degenerate light neutrinos, the lowest allowed value of Mi occurs at to ~ m*, where analytic 
approximations are unreliable. It was shown in [129,209,210] that including flavour effects does not relax 
the bound with respect to the unflavoured case [60,212], and we still have 

Mi > 2 x 10 9 GeV, non - degenerate m*. (9.28) 

This is shown in Fig. ULl [210]. 
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10 Variations 



The bulk of this review focuses on thermal lcptogcncsis, with hierarchical singlet neutrinos, and where 
the CP asymmetry is generated in the decays of the lightest singlet neutrino. In this chapter, we give 
brief overviews of various alternative scenarios of leptogenesis. More detailed presentations can be found 
in the original literature, to which we give appropriate references. 

10.1 Supersymmetric thermal leptogenesis 

In the framework of supersymmetric seesaw models, new leptogenesis mechanisms become possible, such 
as Affleck-Dine leptogenesis or soft leptogenesis. In this subsection, we focus on the standard thermal 
leptogenesis scenario with hierarchical singlet iV's, where the presence of supersymmetry (SUSY) makes 
only small qualitative and quantitative differences [54,225]. We have often referred to the supersymmetric 
modifications in previous chapters, so here we collect this information and summarize the differences that 
arise when calculating the baryon asymmetry in the supersymmetric type I seesaw. 

Neglecting small SUSY-breaking terms, a singlet neutrino Ni and its superpartner, the singlet sneu- 
trino Ni, have equal masses, equal decay rates = Tjj. and equal CP asymmetries ejV; = ■ In 
this approximation, estimating the B — L asymmetry in the minimal supersymmetric standard model 
(MSSM) with respect to the SM case (at fixed values of the Yukawa couplings and M Ni masses) amounts 
to a (careful) counting of a few numerical factors: 

1. There are twice the number of states running in the loops, and thus the CP asymmetries and 
e^. are roughly twice the SM value of (see the loop function g(xj) in eqn (|5.15p ). 

2. Including the asymmetry generated in Ni decays gives another enhancement of a factor of two. 

3. In the MSSM, the plasma is populated by twice as many particles as there are in the SM. More 
precisely, we have g^ SSM /gf M — 228.75/106.75 s» 2.14. Thus the lepton asymmetries to entropy 
ratios are reduced by 1/2. 

4. Due to the additional £<fi final states, the N decay rate is twice faster than in the SM. (For the scalar 
neutrinos there are also two channels jVj — > £(fi,lcf).) In the strong washout regime, the associated 
inverse-decay reactions double the washout rates, reducing the asymmetry by a factor of two. In 
the weak washout regime, the production of JV, and Ni is more efficient, and this increases the 
asymmetry by the same factor. 

5. The expansion rate of the Universe is proportional to ^fgl and is therefore roughly a factor of 
y/2 faster in SUSY (mf SSM w V2mf M ). This reduces the time during which the strong washout 
processes can erase the asymmetry, yielding a V2 enhancement. For weak washouts, it reduces the 
time for N and N production, yielding a y/2 suppression of the final asymmetry. 

6. Finally, in both the SM and the MSSM cases, the asymmetry in baryons Y/\b is of order 1/3 of the 
asymmetry in B — L. The exact relation differs only slightly: 28/79 for the SM vs. 8/23 for the 
MSSM (see eqns (|14.17|) and (|14.18|) , respectively). 

Putting all these factors together, we estimate: 



Thus, in the MSSM, in spite of the doubling of the particle spectrum and of the large number of new 
processes involving superpartners [186], one does not expect major numerical changes with respect to the 
non-supersymmetric case. 

Supersymmetric thermal leptogenesis has a potential gravitino problem [76,77,79,80]. After inflation, 
the universe thermalizes to a reheat temperature T rc hoat- Gravitinos are produced by thermal scattering 
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(weak washout). 



(10.1) 
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in the bath, and the rate is higher at higher temperatures. The gravitinos are long-lived; if there are 
lighter SUSY particles (the gravitino is not the LSP), the decay rate can be estimated as [79] 



If too many gravitinos decay during or after Big Bang Nucleosynthesis (at t ~ seconds), the resulting 
energetic showers in the thermal bath destroy the agreement between predicted and observed light element 
abundances [75,78,81]. There are several possible solutions to the gravitino problem: 

1. The reheat temperature is low enough, T re h ea t 10 6 — 10 10 GeV [81,226] that the gravitino density 
is small. The precise bound depends on the gravitino mass and SUSY spectrum. 

2. The gravitinos decay before BBN; sufficiently heavy gravitinos can arise in anomaly-mediated sce- 



3. Late time entropy production can dilute the gravitino abundance, but it also dilutes the B — L 
asymmetry [228]. 

4. The gravitino is the LSP (and the dark matter), in which case the bound on T re heat is less restrictive 
[229,230]. The gravitino can be the LSP in gauge-mediated scenarios. 

As discussed in Section I5T21 there is a lower bound on the reheat temperature T re h ea t <^ 10 9 GeV, to 
produce a big enough baryon asymmetry by thermal leptogenesis (in the three generation type I seesaw 
with hierarchical Ni masses). This is difficult to reconcile with the first solution described above [231], 
which is arguably the most plausible one in gravity mediated scenarios. A large enough baryon asymmetry 
can be produced at lower T re h ea t, if Ma — Mi ^ T_d, the N± decay width. 

10.2 Less hierarchical iV's 

It is crucial for thermal leptogenesis that there are more than one singlet neutrino. In particular, the 
CP asymmetry that is produced in Ni decays comes from interference between the tree diagram and 
the loop diagrams involving the heavier Ni as virtual particles in the loop. In the conventional picture, 
however, one assumes a strong hierarchy, Mj>i 3> M\. Then, an upper bound on the CP asymmetry 
in N\ decays applies (see eqn (|5.10p ). which provides much of the predictive power of the conventional 
leptogenesis scenario. 

However, as mentioned in section f5.2[ when the singlet neutrino masses are not very strongly hierar- 
chical, a term that is higher order in Mi/Mj>i can become important. Specifically, there is a contribution 
to the CP asymmetry that is proportional to Mf / '(M^M^) and does not vanish when the light neutrino 
masses are degenerate [162,163]. 

In addition to this effect of the heavier singlet neutrino masses on the CP asymmetry in the N\ decays, 
there are also interesting effects related to the decays of the heavier singlet neutrinos A^,... themselves. 
These effects are discussed in the following subsection. 

10.2.1 N 2 effects 

When analyzing leptogenesis, the effects of A2 and A3 are often neglected. This is reasonable if T ro h oa t < 
M2,3. However, if T ro hcat > M^, one should not assume that the L-violating interactions of Ni would 
washout any lepton asymmetry generated at temperatures T ^> Mi and, in particular, the asymmetry 
generated in the decays of ^2,3. Under various (rather generic) circumstances, the lepton asymmetry 
generated in Af 2 ,3 decays survives the ATi leptogenesis phase. Thus, it is quite possible that the lepton 
asymmetry relevant for baryogenesis originates mainly (or, at least, in a non-negligible part) from A?2,3 
decays. 

The possibility that A2 leptogenesis can successfully explain the baryon asymmetry of the Universe 
has been shown in two limiting cases: 




3 



(10.2) 



narios [227]. 
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1. The "iVi-decoupling" scenario, in which the Yukawa couplings of N± are simply too weak to washout 
the iV 2 -generated asymmetry [187,208,232]. 

2. The "strong A^-coupling" scenario, where iVi-related decoherence effects project the lepton asym- 
metry from N2 decays onto a flavour direction that is protected against N\ washout [65,71,188,207]. 

It is plausible that the role of N2 leptogenesis cannot be ignored also in the intermediate range for Ni 
couplings, but the analysis in this case is complicated and has not been carried out yet. 

The TVi-decoupling scenario is simple to understand. It applies when Ni is weakly coupled to the 
lepton doublets (see eqn (|4.7|) : 

fhi <C m*, (10.3) 

where = 1.66gl^ 2 8irv 2 /m p i w 10~ 3 eV and 

(10.4) 

In this case, the asymmetry generated in thermal Ni leptogenesis is too small. Furthermore, the Ni 
washout effects are negligible and, consequently, the asymmetry generated in N2 decays survives. 

We next discuss N2 leptogenesis in the strong TVi-coupling limit. We are interested in the case that a 
sizeable asymmetry is generated in N2 decays, while Ni leptogenesis is inefficient. We thus assume that 
the A^-related washout is not too strong, while the iVi-related washout is so strong that it makes Ni 
leptogenesis fail: 

7712^771*, mi »m,. (10.5) 

To further simplify the analysis, we impose two additional conditions [188]: thermal leptogenesis, and 
strong hierarchy, M 2 /Mi. These two conditions together guarantee that 

n Nl (T ~ M a ) « 0, n N2 (T ~ M x ) w 0. (10.6) 

Thus, the dynamics of N 2 and Ni are decoupled: there are neither iVi-related washout effects during N 2 
leptogenesis, nor iV 2 -related washout effects during Ni leptogenesis. 

The N 2 decays into a combination of lepton doublets that we denote by £2- 

|A) = (AtA) i T 1/2 ^A OT .|4). (10.7) 

a 

The second condition in (| 10. 5[) implies that already at T <; Mi the N\ -Yukawa interactions are sufficiently 
fast to quickly destroy the coherence of I2 . Then a statistical mixture of £\ and of the state orthogonal 
to £\ builds up, and it can be described by a suitable diagonal density matrix. On general grounds, one 
expects that decoherence effects proceed faster than washout. In the relevant range, T M 1 , this is also 
ensured by the fact that the dominant C(A 2 ) washout process (the inverse decay £4> — > N\) is blocked 
by thermal effects [60], and only scatterings with top-quarks and gauge bosons, that have additional 
suppression factors of faf and g 2 , contribute to the washout. 

Let us consider the case where both N2 and N± decay at T <; 10 12 GeV, so that flavour effects are 
irrelevant. We also neglect the effects of AL = 2 interactions, which are generically in equilibrium at 
T <; few x 10 12 GeV. (This scenario is interesting to consider, even though it may only occur in a narrow 
temperature range.) A convenient choice for an orthogonal basis for the lepton doublets is (£i,£q,£' ) 
where, without loss of generality, £' satisfies (^0^2) = 0. Then the asymmetry AYf 2 produced in N2 
decays decomposes into two components: 

AY, = c 2 AY, 2 , AY ei = s 2 AY i2 , (10.8) 

where c 2 = | (^0^2) | 2 and s 2 = 1 — c 2 . The crucial point here is that we expect, in general, c 2 7^ and, since 
{£0 l^i) = 0, AYg is protected against Ni washout. Consequently, a finite part of the asymmetry AYg 2 
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from N 2 decays survives through N± leptogenesis. A more detailed analysis [188] finds that AYg 1 is not 
entirely washed out, and that the final lepton asymmetry is given by Yal = (3/2)AY^ = (3/2)c 2 AY^ 2 . 

For 10 9 GeV ^ Mi ^ 10 12 GeV, flavour issues modify the quantitative details, but the qualitative 
picture, and in particular the survival of a finite part of AYg 2 , still hold. On the other hand, for Mi ^ 
10 9 GeV, the full flavour basis {t e ,t^,t. T ) is resolved and, in general, there are no directions in flavour 
space where an asymmetry is protected and Yi 2 can be erased entirely. 

The conclusion is that N 2 and N 3 leptogenesis cannot be ignored, unless at least one of the following 
conditions applies: 

1. The asymmetries and/or the washout factors vanish, €N 2 r l2 ~ and £at 3 ?73 ~ 0. 

2. iVi-related washout is still significant at T 10 9 GeV. 

3. The reheat temperature is below M 2 . 

In all other cases, the A^-rclated parameters play a role in determining the baryon asymmetry of the 
Universe. Consequently, relations between these parameters and neutrino masses are important. Such 
relations were obtained in Ref. [106, 136], and they lead, for the case that the light neutrino masses have 
normal hierarchy, to the following consequences: 

• In the framework with two singlet neutrinos, both N\ and N 2 interactions are in the strong washout 
regime, with both m, > 0.009 eV [106], and at least one > 0.025 eV. 

• In the framework with three singlet neutrinos, at least two iVj's have interactions in the strong 
washout regime, with m, > 0.005 eV and at least one > 0.02 eV. 

The lower bounds are stronger for inverted hierarchy, and even more so in the framework with three 
singlet neutrinos and quasi-degenerate light neutrinos. 



10.2.2 Resonant leptogenesis 

A resonant enhancement of the CP asymmetry in N\ decay occurs when the mass difference between N\ 
and N 2 is of the order of the decay widths. Such a scenario has been termed 'resonant leptogenesis', and 
has benefited from many studies in different formalisms 23 [61,82,85, 131, 132, 184, 189,234-244]. As this 
review focuses on hierarchical Ni, we briefly list the idea and some references. 

The resonant effect is related to the self energy contribution to the CP asymmetry. Consider, for 
simplicity, the case where only N 2 is quasi-degenerate with Ni. Then, the self-energy contribution 
involving the intermediate N 2 , to the total CP asymmetry (we neglect important flavour effects [61]) is 
given by 

/ ,r \ Mi T 2 M 2 AM 2 Tra[(AU)? 2 ] 

e Nl (self -energy) = - — — { ^ )2 - ^ j^^ - (10-9) 

The resonance condition reads 

M 2 -M!=r 2 /2. (10.10) 



In this case 



1 |Jm[(AtA) 



2 



<iail) 

Thus, in the resonant case, the asymmetry is suppressed by neither the smallness of the light neutrino 
masses, nor the smallness of their mass splitting, nor small ratios between the singlet neutrino masses. 
Actually, the CP asymmetry could be of order one. (More accurately, [ejvil < 1/2.) 

With resonant leptogenesis, the Boltzmann Equations are different. The densities of N\ and N 2 are 
followed, since the both contribute to the asymmetry, and the relevant timescales are different. For 
instance, the typical time-scale to build up coherently the CP asymmetry is unusually long, of order 



3 See [233] for a comparison of the different calculations. 
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1/AM. In particular, it can be larger than the time-scale for the change of the abundance of the sterile 
neutrinos. This situation implies that quantum effects in the Boltzmann equations can be significant, in 
the weak or mild washout regime, for resonant leptogenesis [241,242]. 

The fact that the asymmetry could be large, independently of the singlet neutrino masses, opens 
up the possibility of low scale resonant leptogenesis. Models along these lines have been constructed in 
Refs. [61,236,240,245]. It is a theoretical challenge to construct models where a mass splitting as small as 
the decay width is naturally achieved. For attempts that utilize approximate flavour symmetries see, for 
example, [189,237,239,243,246], while studies of this issue in the framework of minimal flavour violation 
can be found in [131,132,244]. 



10.3 Soft leptogenesis 

The modifications to standard leptogenesis due to supersymmetry have been discussed in section 110.11 
The important parameters there are the Yukawa couplings and the singlet neutrino parameters, which 
are all superpotential terms (see eqn l2.21) ): 

W {N) = \LH U N C + ^MN C N C . (10.12) 

Supersymmetry must, however, be broken. In the framework of the supersymmetric standard model 
extended to include singlet neutrinos (SSM+N), there are, in addition to the soft supersymmetry breaking 
terms of the SSM, terms that involve the singlet sneutrinos N, in particular bilinear (B) and trilinear 
(A) scalar couplings. These terms provide additional sources of lepton number violation and of CP 
violation. Scenarios where these terms play a dominant role in leptogenesis have been termed 'soft 
leptogenesis' [86, 87, 247-258] . 

Soft leptogenesis would have taken place even with a single lepton generation. To simplify things we 
work, therefore, in the framework of a single generation SSM+N. The relevant soft supersymmetry terms 
in the Lagrangian are given in eqn (|10.13|) : 

C^l = - (a\L4> u N c + ^BMN C N C + h.c)j . (10.13) 

In addition, the electroweak gaugino masses play a role: 

=-(m 2 \ a 2 \ a 2 +h.c). (10.14) 

Here X 2 (a = 1,2,3) are the SU(2) gauginos. The effects related to Ai, the U(1)y gaugino, are similar 
to (and usually less important than) those of A2, so we do not present them explicitly. 

The Lagrangian derived from eqns (|10.12|) . (|10.13|) and Q10.14|) has two independent physical CP 
violating phases: 

= arg(AB*), (10.15) 
4> w = arg(m 2 5*). (10.16) 

These phases give the CP violation that is necessary to dynamically generate a lepton asymmetry. If we 
set the lepton number of N c and N c to —1, so that A (and AX) are lepton number conserving, then M 
(and BM) violate lepton number by two units. Thus processes that involve A and M give the lepton 
number violation that is necessary for leptogenesis. 

A crucial role in soft leptogenesis is played by the N — mixing amplitude, 

(N\H\Nl) = M 12 - l -T 12 (10.17) 

{M%2 is the dispersive part of the mixing amplitude, while Ti2 is the absorptive part), which induces 
mass and width differences, 

AM M H -M L AT T H -T L 

'= — = — f — ' V =n=^r— ' (iai8) 
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(r is the average width) between the two mass eigenstates, the heavy \Nh), and the light \Nl) 

\N L<H ) = p\N)±q\N^). 
The ratio q/p depends on the mixing amplitude ratio: 
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2M 12 - iT 12 

Also of importance are the decay amplitudes of N or into various final states X: 

A x = (X\H\N), A x = (X\H\N*). 
The decay width of the singlet sneutrino is given by (for |M| ^> \A\) 

| MA 2 1 
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while the mixing parameters are given by 
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The quantity of interest is the CP asymmetry in singlet sneutrino decay, 

_ r(Z) + r(L)-r(Zt)-r(£) 
r(Z) + r(L) + r(Lt) + r(L)' 



1/4 



(10.24) 



where T(X) is the time-integrated decay rate into a final state with leptonic content X. Here L(L) is the 
(anti)lepton doublet and L(L^) is the (anti)slepton doublet. 

Qualitatively, the most special feature of soft lcptogcncsis is that it gets contributions that are related 
to CP violation in mixing. This is a phenomenon that is analogous to the CP violation in K — K mixing, 
where it leads to the Kl — * tttt decays, the process where CP violation was first discovered. A relative 
CP violating phase between Mi 2 and Ti2 of cqn (|10.17p gives [see eqn (|10.19|) ] \q/p\ ^ 1 which, in turn, 
leads to a situation where the mass eigenstates Nh,l are not CP eigenstates. One such contribution is 
given by [86,87] 
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Here m SUS y is the scale of the soft supersymmetry breaking terms (mi ~ m S usy and A, B ^ m SU sy), and 
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(10.26) 



where Af s (A//) are phase space factors for final states involving scalars (fcrmions). At zero temperature, 
A s f = C(m 2 usy /M 2 ), but for temperature at the time of decay that is comparable to the singlet neutrino 
mass we have A s / » (Af s — A/"/)/(A/" s + A//) = C(l). The difference between A// and Af s at finite 
temperature arises from the Pauli blocking of final state fermions and Bose-Einstein stimulation of decays 
into scalars. 

The contribution (I10.25|) stands out among the soft leptogenesis contributions as the only one that 
is linear in the ratio m susy /M. (All other contributions are quadratic in this ratio.) Thus, for M ^> 
10 2 m susy , this is the only contribution that is potentially significant. Indeed, it could account for the 
observed baryon asymmetry if the following conditions are all fulfilled: 
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1. The lightest singlet sneutrino is light enough, M ^ 10 9 GeV. 

2. The Yukawa couplings are small enough, A 10~ 4 . The lighter is M, the smaller the Yukawa 
coupling must be. 

3. The B parameter is well below its naive value, B/m susy & 10~ 3 . The lighter is M, the more 
suppressed B must be. 

For M ^ 10 2 m susy , there are several other contributions that can be significant [86]. All of these 
contributions involve A2 and are therefore proportional to the weak coupling cti. In addition, as mentioned 
above, they are proportional to (m susy / M) 2 . In particular, there are contributions related to CP violation 
in decay (a phenomenon that is analogous to the one giving 1Ze(e') in K — > irir decays and to the vertex 
contribution in standard leptogenesis) , and to CP violation in the interference of decays with and without 
mixing (analogous to S^k m B decays). These two contributions are suppressed by OL2(jn susy / M) 2 but 
no other small factors. In particular, unlike the contribution e mix of eqn (|10.25[) . these contributions do 
not vanish when i> 1 and therefore they allow B ~ m susy . 

Soft leptogenesis is an interesting scenario in the framework of the supersymmetric seesaw for several 
reasons. First, the relevant new sources of CP violation and lepton number violation appear generically 
in this framework. In this sense, soft leptogenesis is qualitatively unavoidable in the SSM+N framework, 
and the question of its relevance is a quantitative one. Second, if M 10 9 GeV (in the supersymmetric 
framework, this range is preferred by the gravitino problem), then standard leptogenesis encounters 
problems, while soft leptogenesis can be significant. Third, if M ^ 10 2 m susy , then it is almost unavoidable 
that soft leptogenesis plays an important role. 

10.4 Dirac leptogenesis 

The extension of the Standard Model with singlet neutrinos allows two different ways of giving the active 
neutrinos their very light masses. First, one can invoke the seesaw mechanism which gives Majorana 
masses to the light neutrinos. This extension has at least three attractive features: 

• No extra symmetries (and, in particular, no global symmetries) have to be imposed. 

• The extreme lightness of neutrino masses is linked to the existence of a high scale of new physics, 
which is well motivated for various other reasons {e.g. gauge unification). 

• Lepton number is violated, which opens the way for leptogenesis. 

The second way is to impose lepton number and give Dirac masses to the neutrinos. A-priori, one might 
think that all three attractive features of the seesaw mechanism are lost. Indeed, one must usually impose 
additional symmetries. But one can still construct natural models where the tiny Yukawa couplings that 
are necessary for small Dirac masses are related to a small breaking of a symmetry. What is perhaps most 
surprising is the fact that leptogenesis could proceed successfully even if the neutrinos are Dirac particles 
and lepton number is not broken [259,260] (except, of course, by the sphaleron interactions). Such 
scenarios have been termed 'Dirac leptogenesis' [260-266]. Actually, the success of Dirac leptogenesis is 
closely related to the extreme smallness of the neutrino Yukawa couplings in this scenario. 

An implementation of the idea is the following. A CP violating decay of a heavy particle can result 
in a non-zero lepton number for left-handed particles, and an equal and opposite non-zero lepton number 
for right-handed particles, so that the total lepton number is zero. For the charged fermions of the 
Standard Model, the Yukawa interactions are fast enough that they quickly equilibrate the left-handed 
and the right-handed particles, and the lepton number stored in each chirality goes to zero. This is not 
true, however, for Dirac neutrinos. The size of their Yukawa couplings is A 10~ n , which means that 
equilibrium between the lepton numbers stored in left-handed and right-handed neutrinos will not be 
reached until temperatures fall well below the electroweak breaking scale. To see this, note that the rate 
of the Yukawa interactions is given by T\ ~ A 2 T. It becomes significant when it equals the expansion 
rate of the Universe, H ~ T 2 /m p \. Thus, the temperature of equilibration between left-handed and 
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right-handed neutrinos can be estimated as T ~ A 2 m p i ~ (A/10 _8 )Tewpt- By this time, the left-handed 
lepton number has been partially transformed into a net baryon number by the sphaleron interactions. 

More specifically, consider a situation where the CP violating decays of some new, heavy particles 
have produced a negative lepton number in left-handed neutrinos, and a positive lepton number, of 
equal magnitude, in right-handed neutrinos. The sphalerons interact with the left-handed neutrinos, 
violating B + L and conserving B — L. Consequently, part of the negative lepton number stored in 
the left-handed neutrinos is converted to a positive baryon number. At much lower temperatures, when 
sphaleron interactions (and B and L violation) are frozen, the remaining negative lepton number in left- 
handed neutrinos equilibrates with the positive lepton number stored in the right handed-neutrinos. Since, 
however, the negative lepton number in left-handed neutrinos has become smaller in magnitude than the 
positive lepton number in right-handed neutrinos, a net lepton number remains. The final situation is 
then a Universe with total positive baryon number, total positive lepton number, and B — L = 0. 

A specific example of a natural, supersymmetric model with Dirac neutrinos is presented in Ref. [261]. 
The Majorana masses of the ./V-superfields are forbidden by U(1)l (lepton number symmetry). The 
neutrino Yukawa couplings are forbidden by a U(1)n symmetry where, among all the SSM+N fields, 
only the N superfields are charged. The symmetry is spontaneously broken by the vacuum expectation 
value of a scalar field x that can naturally be at the weak scale, (x) ~ v u . This breaking is communicated 
to the SSM+N via extra, vector-like lepton doublet fields, <f> + <fi, that have masses much larger than 
v u . Consequently, the neutrino Yukawa couplings are suppressed by the small ratio (%)/M^. The CP 
violation arises in the decays of the vector-like leptons, whereby T(<f» — ► NH^) ^ T((f> — ► N C H U ) and 
T((f> — > Lx) r(0 — * L c x c )- The resulting asymmetries in N and in L are equal in magnitude and 
opposite in sign. 

The main phenomenological implications of Dirac leptogenesis is the absence of any signal in neutri- 
noless double beta decays. 



10.5 Triplet scalar leptogenesis 

As explained in Subsection 12. 2[ one can generate see-saw masses for the light neutrinos by tree-level 
exchange of SU (2)-triplet scalars. Since this mechanism necessarily involves lepton number violation, 
and allows for new CP violating phases, it is interesting to examine it as a possible source of leptogenesis 
[111,249,254,267-281]. One obvious problem in this scenario is that, unlike singlet fermions, the triplet 
scalars have gauge interactions that keep them close to thermal equilibrium at temperatures T ^ 10 15 
GeV. It turns out, however, that successful leptogenesis is possible even at a much lower temperature. 
This subsection is based in large part on Ref. [274], where further details and, in particular, an explicit 
presentation of the relevant Boltzmann equations can be found. 

The CP asymmetry that is induced by the triplet scalar decays is defined as follows: 

e T . 2 r(f^)-r(r^ 
Tt + r T 

The overall factor of 2 comes because the triplet scalar decay produces two (anti)leptons. 

To calculate ex, one should use the Lagrangian terms given in eqn (|2.15|) . While a single triplet is 
enough to produce three light massive neutrinos, there is a problem in leptogenesis if indeed this is the 
only source of neutrinos masses: The asymmetry is generated only at higher loops and in unacceptably 
small. 

It is still possible to produce the required lepton asymmetry from a single triplet scalar decays if there 
are additional sources for the neutrino masses, such as type I, type III, or type II contributions from 
additional triplet scalars. Define mu (mi) as the part of the light neutrino mass matrix that comes (does 
not come) from the contributions of the triplet scalar that gives ex- 

m = mu + mi. (10.28) 
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Then, assuming that the particles that are exchanged to produce mi are all heavier than T, the CP 
asymmetry is given by 

e T = -. fVB L B H [ \ 11 (10.29) 

where Bl (Bh) is the tree-level branching ratio to leptons (Higgs doublets). If these are the only decay 
modes, i.e. B L + B H = 1, then B L /B H = Tt(X L X\)/(X H A^). There is an upper bound on this 
asymmetry: 

h ^ ^^y^E<- ( 10 - 3 °) 

Note that, unlike the singlet fermion case, \ex\ increases with larger m Vi . 

As concerns the efficiency factor, it can be close to maximal, rj ~ 1, in spite of the fact that the gauge 
interactions tend to maintain the triplet abundance very close to thermal equilibrium. There are two 
necessary conditions that have to be fulfilled by the two decay rates, T — » ££ and T — > 00, in order that 
this will happen [274]: 

1. One of the two decay rates is faster than the TT annihilation rate. 

2. The other decay mode is slower than the expansion rate of the Universe. 

The first condition guarantees that gauge scatterings are ineffective: the triplets decay before annihilating. 
The second condition guarantees that the fast decays do not lead to a strong washout of the lepton 
asymmetry: lepton number is violated only by the simultaneous presence of T — > 11 and T — > 00. 

Combining a calculation of r\ and the upper bound on the CP asymmetry (|10.30|) , successful leptoge- 
nesis implies a lower bound on the triplet mass Mt varying between 10 9 GeV and 10 12 GeV, depending 
on the relative weight of m\i and m\ in the light neutrino mass. 

Interestingly, in the supersymmetric framework, "soft leptogenesis" , namely one that is driven by the 
soft supersymmetry breaking terms, can be successful even with the minimal set of extra fields - a single 
T + T - generating both neutrino masses and the lepton asymmetry [249, 254]. 



10.6 Triplet fermion leptogenesis 

As explained in Subsection 12.31 one can generate see-saw masses for the light neutrinos by tree-level 
exchange of S l f7(2)-triplct fcrmions. This mechanism necessarily involves lepton number violation, and 
allows for new CP violating phases so we should examine it as a possible source of leptogenesis [162,270, 
282]. This subsection is based in large part on Ref. [270], where further details and, in particular, an 
explicit presentation of the relevant Boltzmann equations can be found. 

As concerns neutrino masses, the formalism and qualitative features are very similar to the singlet 
fermion case. As concerns leptogenesis, there are, however, qualitative and quantitative differences. 

With regard to the CP asymmetry from the lightest triplet fermion decay, the relative sign between 
the vertex loop contribution and the self-energy loop contribution is opposite to that of the singlet fermion 
case. Consequently, in the limit of strong hierarchy in the heavy fermion masses, the asymmetry in triplet 
decay is 3 times smaller than in the singlet decay (see the discussion in Section 15.21 for the singlet case 
and for definitions). On the other hand, since the triplet has three components, the ratio between the 
final baryon asymmetry and £77 is 3 times bigger. The decay rate of the heavy fermion is the same in 
both cases. This, however, means that the thermally averaged decay rate is 3 times bigger for the triplet, 
as is the on-shell part of the AL — 2 scattering rate. 

A significant qualitative difference arises from the fact that the triplet has gauge interactions. The 
effect on the washout factor r\ is particularly significant for m <C 10 -3 eV, the so-called "weak washout 
regime" (note that this name is inappropriate for triplet fermions). The gauge interactions still drive 
the triplet abundance close to thermal equilibrium. A relic fraction of the triplet fermions survives. The 
decays of these relic triplets produce a baryon asymmetry, with 

77 w Mi/10 13 GeV (for fh < 10~ 3 eV). (10.31) 
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The strong dependence on Mi is a result of the fact that the expansion rate of the Universe is slower at 
lower temperatures. On the other hand, for rh 3> 1CP 3 eV, the Yukawa interactions are responsible for 
keeping the heavy fermion abundance close to thermal equilibrium, so the difference in 77 between the 
singlet and triplet case is only 0(1). 

Ignoring flavour effects, and assuming strong hierarchy between the heavy fermions, Ref. [270] obtained 
a lower bound on the mass of the lightest heavy triplet fermion: 

Mi £ 1.5 x 10 10 GeV (10.32) 

When the triplet fermion scenario is incorporated in a supersymmetric framework, and the soft break- 
ing terms do not play a significant role, the modifications to the above analysis is by factors of 0(1). 
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11 Conclusions 



During the last few decades, a large set of experiments involving solar, atmospheric, reactor and ac- 
celerator neutrinos have converged to establish that the standard model neutrinos are massive. The 
seesaw mechanism extends the standard model in a way that allows neutrino masses. It provides a nice 
explanation of the suppression of the neutrino masses with respect to the electroweak breaking scale, 
which is the mass scale of all standard model charged fermions. Furthermore, without any addition or 
modification, it can also account for the observed baryon asymmetry of the Universe. The possibility of 
giving an explanation of two apparently unrelated experimental facts - neutrino masses and the baryon 
asymmetry - within a single framework that is a natural extension of the standard model, together with 
the remarkable 'coincidence' that the same neutrino mass scale suggested by neutrino oscillation data 
is also optimal for leptogenesis, makes the idea that baryogenesis occurs through leptogenesis a very 
attractive one. 

Leptogenesis can be quantitatively successful without any fine-tuning of the seesaw parameters. Yet, 
in the non-supersymmetric seesaw framework, a fine-tuning problem arises due to the large corrections 
to the mass-squared parameter of the Higgs potential that are proportional to the heavy Majorana 
neutrino masses (see Section 12. 4p . Supersymmetry can cure this problem, avoiding the necessity of 
fine tuning. However, the gravitino problem that arises in many supersymmetric models requires a low 
reheat temperature after inflation, in conflict with generic leptogenesis models (see Section ["l0.1[) . Thus, 
constructing a fully satisfactory theoretical framework that implements leptogenesis within the seesaw 
framework is not a straightforward task. 

From the experimental side, the obvious question to ask is if it is possible to test whether the baryon 
asymmetry has been really produced through leptogenesis. Unfortunately it seems impossible that any 
direct test can be performed. To establish leptogenesis experimentally, we need to produce the heavy 
Majorana neutrinos and measure the CP asymmetry in their decays. However, in the most natural seesaw 
scenarios, these states are simply too heavy to be produced, while if they are light, then their Yukawa 
couplings must be very tiny, again preventing any chance of direct measurements. 

The possibility of indirect tests from measurements of the asymmetries produced by leptogenesis is also 
ruled out. This is because there are too many high energy parameters that are relevant to leptogenesis 
and, even if we adopt the most optimistic point of view, there are at most four observables. These 
are the value of the baryon asymmetry that is known with a good accuracy (see Section and the 
three cosmic neutrino flavour asymmetries that we do not know how to measure. A measurement of 
the total lepton asymmetry would be very valuable. If a value of the order of the baryon asymmetry 
were found, it would provide strong evidence that electroweak sphalerons have been at work in the early 
Universe. 24 Even if we have no reason to doubt that sphaleron processes occurred in the thermal bath 
with in-equilibrium rates, they arc the fundamental ingredient of the whole leptogenesis idea, and thus 
the importance of an experimental test should not be underappreciated. Measurements of the single 
lepton flavour asymmetries would provide some information on leptogenesis parameters (on the flavour 
CP asymmetries if leptogenesis occurred in the unflavoured regime, or on a combination of the flavour 
CP asymmetries and of the flavour-dependent washouts otherwise). However, at present, we have not 
even observed the relic neutrino background, and the possibility of revealing 0(1O -10 ) asymmetries in 
this background seems completely hopeless. 

Lacking the possibility of a direct proof, experiments can still provide circumstantial evidence in 
support of leptogenesis by establishing that (some of) the Sakharov conditions for leptogenesis (see 
Section fl.2|) are realized in nature. 

Planned neutrinoless double beta decay (Ovj30) experiments (GERDA [285], MAJORANA [286], 
CUORE [287]) aim at a sensitivity to the effective Qvj3(3 neutrino mass in the few x 10 meV range. 
If they succeed in establishing the Majorana nature of the light neutrinos, 25 this will strengthen our 
confidence that the seesaw mechanism is at the origin of the neutrino masses and, most importantly, will 

24 The opposite result would not disprove sphaleron physics: see [283, 284] for a model in which production of a large 
lepton asymmetry after sphalerons freeze-out can yield Yal orders of magnitude larger than Yab- 
25 This is likely to happen if neutrinos are quasi-degenerate or if the mass hierarchy is inverted 
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establish that the first Sakharov condition for the dynamical generation of a lepton asymmetry, that is 
that lepton number is violated in nature, is satisfied. 

Proposed SuperBeam facilities [288, 289] and second generation off- axis SuperBeams experiments 
(T2HK [290], NO^A [291]) can discover CP violation in the leptonic sector. These experiments can only 
probe the Dirac phase of the neutrino mixing matrix. They cannot probe the Majorana low energy or 
the high energy phases, but the important point is that they can establish that the second Sakharov 
condition for the dynamical generation of a lepton asymmetry is satisfied. 26 

In contrast to the previous two conditions, verifying that the decays of the heavy neutrinos occurred 
out of thermal equilibrium (the third condition) remains out of experimental reach, since it would require 
measuring the heavy neutrino masses and the size of their couplings. 

Given that we do not know how to prove that leptogenesis is the correct theory, we might ask if there 
is any chance to falsify it. Indeed, future neutrino experiments could weaken the case for leptogenesis, 
or even falsify it, in two main ways: Establishing that the seesaw mechanism is not responsible for the 
observed neutrino masses, or finding evidence that the washout is too strong. 

By itself, failure in revealing signals of Quf3[3 decays will not disprove leptogenesis. Indeed, with 
normal neutrino mass hierarchy one expects that the rates of lepton number violating processes are below 
experimental sensitivity. However, if neutrinos masses are quasi-degenerate or inversely hierarchical, and 
future measurements of the oscillation parameters will not fluctuate too much away from the present best 
fit values, the most sensitive 0v(3(3 decay experiments scheduled for the near future should be able to 
detect a signal [72]. 27 If instead the limit on \mp^\ is pushed below ~ lOmeV (a quite challenging task), 
this would suggest that either the mass hierarchy is normal, or neutrinos are not Majorana particles. The 
latter possibility would disprove the seesaw model and standard leptogenesis (see however section fL0.4|) . 
Thus, determining the order of the neutrino mass spectrum is extremely important to shed light on the 
connection between Qvf3(3 decay experiments and leptogenesis. 

Long baseline neutrino experiments can achieve this result [293, 294] by exploiting the fact that os- 
cillations in matter can determine the sign of the atmospheric mass difference Ato 2 3 . I n fact, if the 
oscillations of atmospheric neutrinos involve the heaviest (lightest) state, corresponding to normal (in- 
verted) hierarchy, then <-> v e oscillations are enhanced (suppressed) while «-> z> e , oscillations are 
suppressed (enhanced). 

Cosmology could also provide important information in establishing if the neutrino mass hierarchy is 
inverted. Cosmological observations arc sensitive to the sum of the three neutrino masses m cosmo = m\ + 
?Ti2 + 7713. In the near future, improved CMB data and more precise large scale structure measurements 
can push the sensitivity on m cosmo down to the 0.1 eV level. For an inverted hierarchy a signal should be 
detected around or above this value. 

In summary, if it is established that the neutrino mass hierarchy is inverted and at the same time no 
signal of Oz^/3/3 decays is detected at a level \m ee \ & lOmeV, one could conclude that the seesaw is not 
at the origin of the neutrino masses, and that leptogenesis is not the correct explanation of the baryon 
asymmetry. 

The neutrino mass scale has also more direct implications for leptogenesis. A quantitatively successful 
leptogenesis prefers a mass scale for the light neutrinos not larger than a few times the atmospheric neu- 
trino mass scale (see Section l9.4.ip . In particular, in the standard type I seesaw model (see section |2~T]) . 
if leptogenesis occurs in the unflavoured regime and the heavy Majorana neutrinos are sufficiently hi- 
erarchical, leptogenesis would fail to produce enough baryon asymmetry when the neutrino mass scale 
exceeds values as low as ~ 0.1 — 0.2 eV. Indeed, if past experiments had detected neutrino masses of 
order a few eV, leptogenesis, if not completely ruled out, would have certainly lost most of its theoretical 
appeal. In the future, a discovery of a neutrino mass at the /3 decay experiment KATRIN will also put 

26 As discussed in [292], it is conceivable (but not natural) that the Dirac phase is solely responsible for the matter- 
antimatter asymmetry. 

27 0v(3f3 decay experiments are sensitive to the effective mass parameter nupp = ^2i^ei m i' ^ n * ne mver ted hierarchy 
scenario the s 2 3 rri3 contribution to mpp can be neglected, and approximating for simplicity c 2 3 ~ 1 one gets mpp n 
mie^^Cjj + m2e 2m Sj 2 . Since 612 deviates from maximal mixing, we can then predict that mpp should be of the order of 
the atmospheric mass scale, independently of the values of the Majorana phases a and f). 
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the type I seesaw scenario for leptogenesis in serious doubt, since it will yield a mass value a bit too 
large (the claimed discovery potential is for an effective /?-decay mass parameter = \ U e i\ 2 rnf of 
about 0.35 eV [295]). However, cosmological observations have already reached a sensitivity to m cosmo 
somewhat below 1 eV [296-298] , implying that it is unlikely that leptogenesis could be put to doubt by 
a direct detection of neutrino masses in laboratory experiments. 

As concerns CP violation, a failure in detecting leptonic CP violation will not weaken the case for 
leptogenesis in a significant way. Instead, it would mean that the Dirac CP phase (or the #13 mixing 
angle, or both) is small enough to render CP violating effects unobservable. 

Finally, the CERN LHC has the capability of providing information that is relevant to leptogenesis. 
First, electroweak baryogenesis (see Section ll.2[) can be tested at the LHC. It will become strongly 
disfavored (if not completely ruled out) if supersymmetry is not found, or if supersymmetry is discovered 
but the stop and/or the Higgs are too heavy. Eliminating various scenarios that are able to explain the 
baryon asymmetry will strengthen the case for the remaining viable possibilities, including leptogenesis. 
Conversely, if electroweak baryogenesis is established by the LHC (and EDM) experiments, the case for 
leptogenesis will become weaker. 

Second, new physics discoveries at the LHC can play a fundamental role in establishing that the 
origin of the neutrino masses is not due to the seesaw mechanism, leaving no strong motivation for 
leptogenesis. This may happen in several different ways. For example (assuming that the related new 
physics is discovered), the LHC will be able to test if the detailed phenomenology of any of the following 
models is compatible with an explanation of the observed pattern of neutrino masses and mixing angles: 
supersymmetric R-parity violating couplings and/or L- violating bilinear terms [299,300]; leptoquarks [301, 
302]; triplet Higgses [303,304]; new scalar particles of the type predicted in the Zee-Babu [305,306] types 
of models [307-309]. It is conceivable that such discoveries can eventually exclude the seesaw mechanism 
and rule out leptogenesis. 

To conclude, the seesaw framework provides the most natural and straightforward explanation of 
the light neutrino masses and has, in principle, all the ingredients that are necessary for successful 
leptogenesis. This makes leptogenesis arguably the most attractive explanation for the observed baryon 
asymmetry. This scenario has limited predictive power for low energy observables, so it is unlikely to 
be directly tested. Yet, future experiments have the potential of strengthening, or weakening, or even 
falsifying the case for leptogenesis. 
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12 Appendix: Notation 
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(.ilillUllolUIllL-cih Lllllc Val IdUlc 1V1\ I ± 


IIO. 1UI 


\ ' 

I 


\sti\l z)\ai j az) 


Ib.lUI 


1 I 1 I 


Hubble rate, Hubble rate at 1 = Mi 


11 O TI 11 O 111 

lid. 71 113.111 


s 


entropy density 


\l6.b\ 


Y x 


ratio ol A number density to entropy density 


113.61 


cq 

n x 


number density lor particle X m kinetic equilibrium 


113.21113.51 


Y x 


Yx in kinetic equilibrium 


10.31 113. Ol 


vcq 


ratio of v a and e a density to entropy density ( ge = 2) 


alt ei 10. ol 


yaa 


density or doublets and singlets 01 flavour a 


f6?fl 


VX 


Y x/Y x 


lo.3l 




yx - yx 


10.31 




ratio ol B / 6 — L a asymmetry to entropy density 


\A Al rTTTTTl 

16.91114.111 


Y A b 


ratio ol baryon asymmetry to entropy density 


11. 21 14. 11116. 101 


9*i 9*s 


degrees ol Ireedom m the thermal bath 


110 oi 1 1 ni 

113.81 Ho.91 


9i 


internal degrees of freedom of particle i 


1 1 „ 1 1 ol 

belowllo.dl 


Va 


efficiency parameter for flavour a 


n n 11 >c ti 11 c 01 
I4.ll lib. n llb.ol 




CP asymmetry m flavour a 


I r 11 1 c: 101 

15.11 15.131 


r i 

[to] 


liglit neutrino mass matrix 


12.21 




largest, smallest liglit neutrino mass 


1 C 1 Al IC 1 71 

I5.1UI 15.171 


rh, m aa , m* 


(rescaled) N decay rates, Hubble expansion rate 


14.71 


A 


matrix relating the Y/\i a to the Ya q 


I14.14II14.15I 
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neutrino Yukawa coupling 


lo ol 
I2.3I 


rl x 


1 Lillet W a LU LlJJllllfc Ul Lidl LlClC X I11UL a LLll 111U 1 




V u 


XllggS VeV [V U — 1/4 KjG V J 


lO ql 


1 L> 


LOtai cecay rate 2^a. \ — ^P^oli y^ol) 


I/i /ii |i q oni 
14.41 llo.ZUI 




total inverse decay rate 


i/i 1 ni 

I4.1UI 


p 

1 act 


pai Liai ciecay laic lo iiavoui aiiQ aiiLiiiavoui ex. i \i\ — > cpt a: cpz a j 


I4.1UI 


B x 

D xz 


ui aliening lanu 1 ^yv ? y/, ) I X 


■ r= 1 01 
I0.1UI 


1 A B 


rate density for A — > S 


|13.14| 




rate density for two body decays of ATj 


113.211 


r T , 7r 


rate, rate density of r Yukawa interaction 


|14.10| 


A 7 £ 


CP difference: 7^ - 7 J 


[Ql 




difference between process and its time reversed 


l6~4l 




7^ with exchange of on-shell particle removed 


1631 


#i(z), X 2 (z) 


Bessel functions 


113.191113.31 


X, .4 


matrix element, amplitude 


I5~2l 


1 


4-momcntum conservation: (27r) 4 <5 4 (p ai + ... — Pb t — ■■) 


l5~6l 




phase space d 3 p/(16i?7r 3 ) 


EH 



Table 2: Parameters and variables introduced in the text, and equations where they are defined. Simple 
BE for leptogenesis that neglect scattering processes are given in eqns (]6 . 16[) . and (|6 . 28[) . More complete 
BE that include top-quark and gauge bosons scatterings are given in eqns (|6.51|) , (]6.53[) , (|6.54p and (|6.56jl . 
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13 Appendix: Kinetic Equilibrium 



The interested reader can find a more complete set of formulae for example in the Appendices of refs. 
[31,310]. 

13.1 Number densities 

Of the particles relevant to leptogenesis, all but the N^s carry charges under the standard model gauge 
group and are, therefore, subject to gauge interactions. We take the distributions fi(p) of the (non- 
singlet) particle species i to have their kinetic equilibrium form, because gauge interactions are fast in the 
relevant temperature range, fi(p) = (rii/n° q )f° q (p). (The N^s may also be produced in an equilibrium 
distribution because the initial state particles are in kinetic equilibrium, or kinetic equilibrium may be 
obtained due to fast Yukawa interactions, if A is in the strong washout regime.) With some additional 
assumptions (see eqn ll3.12")) . this allows to solve the BE for the total number density rather than mode by 
mode (as was done in [137]). The number densities of particles and anti-particles are considered separately. 
For some other applications, it may be the convention to count the particles and anti-particles together 
(for instance for g*,g*s> see eqns (|13.8I) . (|13.9|1 ) so some care with factors of 2 is required in matching. 
Using Maxwell-Boltzmann statistics for 

f^m(p)=e-^-^ T , (13.1) 
the equilibrium number density of particles of species i is given by 

9 — f «bW = { ^B 2 ^ * = ^ °« " = 0, (13.2) 



i,MB " (27r)3 J » "Jt,MB wv - | a£ m _ = Q 

Here Ki{z) is a Bessel function (see rcf. [311] eqns 8.432.3 and 8.472.2, and ref. [59] for useful analytic 
approximations) : 



z 2 K 2 (z)=J xe x y/x 2 - z 2 dx -> j /15 + z ?/wr e -z rl:L]i 



and gi is the number of internal degrees of freedom of the particle. For example, g^ — 2 for the iVj's, 
because they are Majorana fermions, g eR = 1 for the SC/(2)-singlet leptons, and gi — 2 for the SU(2)- 
doublet leptons and g<f, = 2 for the Higgs. Antiparticle number densities have gj = gi. 
Using Fermi-Dirac (+) or Bosc-Einstcin (— ) statistics, 

gives the following equilibrium number densities: 

ffiT 3 f C(3) + f C(2) + ... (bosons) 



H,± 



(2tt) 3 



d?pft%{p)^{ "I |C(3) + f^+ (fermions) " (13 . 5) 

\MB < T ' m i > T ' 



where ((2) = tt 2 /6, C(3) = 1.202, and C(4) = tt 4 /90. Using Maxwell-Boltzmann instead of Fermi-Dirac 
statistics makes a difference in n oq of order 10% at T = m. 

13.2 Rates 

Number densities decrease due to the expansion of the Universe. This can be factored out of the BE by 
taking Yx to be the comoving number density: 

= s= 9 -^T\ (13.6) 

s 45 
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The entropy density s is conserved per comoving volume: 



dS ITT 

— = —SHs 
dt 



The expansion rate of the Universe is given by 

H = 



3m 2 pl 



m p i 



(13.7) 



The density and the pressure of the gas of relativistic particles are given by 



P : 



30 



P 



90 



The (temperature-dependent) parameters g*(g*s) are the effective number of degrees of freedom con- 
tributing to the energy (entropy) density: 



<J* 



g*s 



o E 9f 



fen 



fermions 



gj 



£ 5fc (i;) 



E 5fc 

bosons 



(13.8) 
(13.9) 



where the particle and anti-particle should both appear in the sums over species. (Recall that we defined 
<7i as the internal degrees of freedom for the particle only.) The SM Higgs and anti-Higgs contribute 
50 + 9^ = 4. Each flavour of neutrino and its anti- neutrino contribute 7/4. At leptogenesis temperatures, 
g*s = .<?*• The contribution to the present entropy density sq comes from the photons at T~, and from 
the neutrinos and anti-neutrinos 28 at ~ (4/ll)T 3 . Thus sq oc g*sT^, with g*s = 2 - 



A convenient time variable for leptogenesis in Ni decay is 

z = Mi/T, 

for which 

dz Mi dT H x 



8 11 



3.9. 



(13.10) 
(13.11) 



where H± = H(T = Mi) is a constant. 

Suppose the number density of X-particles can be changed by various interactions. The time evolution 
of Y x is 



dYx 
dt 



E 

int 



dHxfxdn a f a ...\M(X + a 



l sdii i (i±fhdii j (i±f^). 



• • ■ \M(i+j + . . . -> X + a + . . .)\ 2 ~S (1 ± fx)dU x (l ± f a )dU a 



, (13.12) 



where dUx and 5 are defined in eqn (|5.6p . and \A4\ 2 is summed over the internal degrees of freedom of the 
initial and final state. This equation simplifies if the particles are assumed to be in kinetic equilibrium, 
and if the final state Bose enhancement/Pauli blocking factors are ignored. Then one obtains: 



dY x 
dt 



E 

int 



Y x Y a . 



YY 

1 3 



-7(A + a + 



i + ] + ■■■) 



{sYP){ S Y^) . . l{l + J + 



X + a + ...) 



(13.13) 



28 The momentum distribution of the neutrinos freezes out when they are relativistic, so even though their mass may be 
greater than their temperature, they contribute to s like a relativistic species. 
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Here 7 is an interaction density: 

= j dU x f c x q dU a f^ . . . \M(X + a + ...^i+j + .. .)\ 2 S dU.dUj , (13.14) 

where we drop final state phase space factors. Recall that \A4\ 2 is summed (not averaged) over the 
internal degrees of freedom of the initial states. 

The rates in this Appendix are obtained in T = field theory. Calculating rates at finite temperature 
[60] gives O(l) effects on the number densities as a function of time, although it can change individual 
rates in a more significant way, as discussed in Section [7] It is useful, for the purpose of calculating decay 
and scattering rates, to have the expression for two body phase space: 

SdUpdUg = [ J Pp } r dn p 

J 167rVs 

\p P -p q \ \/^p-Pi) 2 - m >\ (1;n . ) 

8TTy/s 47TS 



where the first line has the familiar form of the center of mass frame, while the second provides the 
Lorentz-invariant form. 

Consider the two-body decay rate density, 

A + ifi, = 7(JVi - htM)- ( 13 - 16 ) 

There are four possible final states, vptfi , ep4> + ^pcfp , ep(j)~ . The four rates are all equal: 

\M{N X - v (tP)\ 2 = 2\X px \ 2 p N -p l = \X pi \ 2 (M 2 + m 2 - m 2 H ). 
We thus obtain for the decay rate density: 

d 3 p N 



Hip + 



-b-En/t 4 |A /J1 | 2 Mf J 6 dTl^dlle 



2tt 2 

where K\ is a Bessel function (see ref. [311], eqn 8.432.3): 



2E n (2tt) s 

1 T 2 f°° _ x r- \X pi \ 2 M 2 

= 7T / e x y'x 2 -z 2 dx LJ ^ - 

2 2tt 2 J z V 2tt 

= sY^^-T{Nt ^ Ip^l^) (13.17) 
' /A r \ 2 K 1 {z)T{N 1 tphlrf), (13.18) 



3C 



zK x {z)= I e-*y/x*-z*dx^\ 2_ z (13.19) 



and (<?at = 2) 

r(JVl -> ^0, 1^) = |A f |2Ml . (13.20) 
The inclusive two-body decay rate density appears frequently in our calculations: 

= Y,(l(N - 0^) + 7 (JV -> = ^-^X x ( 2 )r D . (13.21) 
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In the literature it is frequently denoted by 7^. It was evaluated with quantum statistics in [60], in which 
case it does not have a simple closed form, and differs by ~ 10%. 
Consider the scattering rate density, 

ltin = 7{i+3^rn + n). (13.22) 

It is given by 

lt n = J </l I, </l I, I; 7 .// / \M(i + j —> m + n)\ 2 5 dU m dU n 

= 4. 9i9j j dn 4 dn j e-& t+B iV T ,Jfa ■ Pj ) 2 - m 2 m 2 a(( Pi + Pj f), (13.23) 

where the g+gj -factor appears because cross-sections are usually averaged over initial state internal degrees 
of freedom. The square root cancels the initial particle flux factor of the cross-section. It is convenient to 
separate the center of mass from initial state phase space integrals, by multiplying by 1 = j d i Q5' l (Q — 
Pi —Pj)- Using eqn (|13.15p gives: 

, x f dQ d 3 Q e~Q°/ T u 2 2 2 2 
7(«+j->m + n) = g l g j J 4 ^— {{Pi-Pj) ~m imj \a(Q) 



9*9 j T 
32 



^^^(^f 1 '^)^ (13 - 24) 



where A(a, b, c) = (a—b — c) 2 — Abe, zK\{z) is given in eqn (|13.19p . and we used d 3 Q = y/Q 2 , — s ds dSl/2 
to obtain the second equality. 

In the massless limit, \(l,x,y) — > 1 and s can be integrated from to 00, so the definite integral 

x n Kiips) = 2 n ~ 1 r(l + n/2)r(n/2) (13.25) 
is useful for obtaining 7. For interaction terms in the Lagrangian of dimension 4 and 5, one obtains 

{ g, gj T 4 _ 1 

^ ;rr (13.26) 
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14 Appendix: Chemical Equilibrium 



In this section, the conditions of chemical equilibrium are used to derive the A-matrix and the factor of 
12/37 relating (within the Standard Model) AY B to AY B -l- 

Various Standard Model interactions can change the number of particles of different species. For 
instance, the Lagrangian term h a £ a 4> c efi a changes a charged SU (2)-singlet lepton into a Higgs and an 
S'C/(2)-doublet lepton. If such interactions are fast, compared to the expansion rate H, they lead to an 
equilibrium state where the comoving number densities of the participating particles remain constant. 
This is described by conditions of chemical equilibrium (see chapter 10 of ref. [312]): The sum of the 
chemical potentials, over all particles entering the interaction, should be zero. For example, if the charged 
lepton Yukawa interaction is fast, we have 

Me Q - H a + M</> = 0. (14-1) 

Thus, the set of reactions that are in chemical equilibrium enforce algebraic relations between various 
chemical potentials [41,313]. 

The chemical potentials can be related to the asymmetries in the particle number densities, by ex- 
panding the distribution functions of eqn (|13.5p for small n/T: 

fermions (14.2) 

bosons (14-3) 

where gi is the number of degrees of freedom of the particle. 

At the temperatures where the lepton asymmetry is generated, the interactions mediated by the top- 
quark Yukawa coupling h t , and by the SU(3) x SU (2) x U(l) gauge interactions, are always in equilibrium. 
This situation has the following consequences: 

• All components of a gauge multiplet share the same chemical potential, and the chemical potential 
of gauge bosons is zero. 

• Hypercharge neutrality implies 

(m?« + 2 M« 4 -MdJ - (W* ~ Me Q ) + 2/i0 = . (14.4) 

i=l,2,3 a=e,fi,T 

• The equilibrium condition for the Yukawa interactions of the top-quark ji t = fi q3 + /j,^ yields: 

Ay t -^ = ^±, (14.5) 

where the relative factor 1/2 between Ay 93 and Ay^ can be traced back to eqns (|14.2|) - (|14.3|) : the 
relation of the number asymmetry to the chemical potential is different by a factor of 2 between 
bosons and fermions. 

Using these relations, the asymmetries in Bi and L a can also be expressed in terms of the chemical 
potentials, for example, 

T 3 

Yab 3 = gj (2 figs +Mt R . 

Yal t = ^-(2M£ T +MtJ, (14.6) 

and similarly for the lighter generations. We use the notation Yax for the asymmetry Yx — Y Xl because 
we think it is clearer, though not standard: for instance Yb is usually the baryon asymmetry. However, 
this should not cause confusion. 



n f - nj 



n b - ni = 



6 

9iT 3 
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The baryon flavour asymmetries can always be taken to be equal: 



Y 



AB 3 



— Yar„ — Yi 



ABi 



(14.7) 



because they are conserved before flavour-changing quark interactions come into equilibrium, and we 
assume them to be equal as an initial condition. (Once flavour-changing quark interactions enter equilib- 
rium, these asymmetries are driven to be equal.) The conservation rules (|14.7|) impose two conditions on 
the sixteen chemical potentials of the SM fields, and hypercharge neutrality, eqn (|14.4[) . imposes a third. 

We work in the approximation that each of the SM interactions is either negligible, in which case there 
are additional conserved quantum numbers, or in chemical equilibrium. Either case yields a restriction on 
the chemical potentials. The interaction rate densities for processes involving SM Yukawa couplings can 
be obtained from eqn Q13.24p , assuming the Yukawa coupling at one vertex, and a gauge or top Yukawa 
coupling at the other. Then the interaction rate for a Yukawa coupling that brings into equilibrium an 
5E/(2)-singlet r can be estimated as T r ~ j n l{ij — * rn)/Yj q . 

The rate of the charged lepton Yukawa interactions is important for leptogenesis. A careful study [195] 
of the relevant processes, which takes into account thermal masses, finds that the Higgs decay gives the 
fastest rate ~ 10~ 2 /i 2 T. Using the results obtained in [195,314] for a = e, we find that the h T ,hfj,,h e 
Yukawa interactions are in equilibrium for T < 10 12 , 3 x 10 9 , 10 5 GeV, respectively. Here, as an example, 
we estimate the (subdominant) ft T -related scattering rate at zero temperature. It can be calculated from 
interactions such as + £ T — > t r + t R and W + £ T — > t r + H . In the massless limit, the cross-sections 
for tr production satisfy 



v t tr) 



9 Z K 



I6irs ' 



a(y T t a R — > T R b a L ) = a[y T bh — * T R t R a ) = <r(b L a tR a ■ 
where a is a colour index, and 

o{v T W° -> T R <j>+) = o(v T W~ - T R ctP d ) = 

Including all these processes, summing over color and SU(2) indices, and using eqn (|13.26[1 . we obtain 

(9h 2 t + 2g 2 )h 2 T Tyi28ir i 



(14.; 



(14.9) 



r =^- = 
T i 



T 3 /tt 2 



5 x 10 - /if T. 



(14.10) 



The temperature ranges in which the various SM interactions are faster than H are given in table [TJ 
Consider the temperature range 10 12 GeV > T > 10 9 GeV, with SM particle content. The strong 
and electroweak sphalerons, and the top, bottom, tau and charm Yukawa interactions are taken to be in 
chemical equilibrium. This imposes six conditions, in addition to the three conditions of eqns (|14.7|) and 
l|14.4p . Moreover, the asymmetries in fi R , e R , u R — d R , and d R — s R , are conserved and presumably zero. 
So the system of equations can be solved for the three asymmetries 



Y Aa - 3 Y AB - Y AL 
which are conserved in the SM, as a function of the /j,i a : 




-22/9 


-4/9 


-4/9 


/ W e 


-4/9 


-22/9 


-4/9 


Mi* 


-2/9 


-2/9 


-139/45 


\ Mr 



(14.11) 



(14.12) 



However, in the range of temperatures that we consider, 
potentials should be summed, /i = + jit , giving 



and £ e are indistinguishable, so their chemical 



Y Aa + Ya„ 

Y At J 6s 



-26/9 -8/9 
-2/9 -139/45 



too 

to-r 



(14.13) 
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The A- Yukawa interactions can generate an asymmetry in the lepton doublets. However, SM in- 
teractions rapidly change the doublet densities, so it is better to write Boltzmann Equations for the 
asymmetries Ya„ , which are conserved by all the SM interactions. The N interactions can add a lepton 
to the asymmetry, or remove a doublet lepton i a . Since the washout interactions can only destroy the 
part of the B/3 — L a asymmetry stored in t a , we need to express the Y^ a ~ 2 r ig a T 3 /(6s) asymmetry 
as a function of the charges Y"a« eqn (|14.11[) . Inverting eqn (|14.13[) gives 



Yai 
Yau 



-417/589 120/589 
30/589 -390/589 



Ya c + Y A „ 
Ya t 



(14.14) 



and the matrix is called "the ^-matrix" [65]. 

At temperatures below ~ 10 9 GeV, the strange and muon Yukawa interactions are in equilibrium, so 
the asymmetries in d,R — sr, and in /i^ are no longer conserved. There are three distinguishable flavours, 
and the A-matrix is: 




-151/179 20/179 20/179 
25/358 -344/537 14/537 
25/358 14/537 -344/537 




(14.15) 



These expression for the A-matrix agree with refs. [63,64]. Ref. [63] defines Ya ( per gauge degree of 
freedom, that is 1/2 of what we use here. Therefore the ^-matrices they quote differ by a factor of 2. 
Our expressions differ from ref. [62], where the u and d Yukawa interactions are taken in equilibrium, 
and from ref. [65], where the strong sphalerons are not included. 

The coefficient relating Yab to the Ya„ should be calculated at the temperature when the sphalerons 
go out of equilibrium. The reason is that, after the period of leptogenesis, the A a 's are conserved. 
Baryon number B is, however, not conserved in the presence of the EW sphalerons (see eqn l6.13j) . and 
the relation of Yab to the Ya q depends on which other interactions are in equilibrium. Using the chemical 
equilibrium and charge conservation equations, one obtains the following relation, depending on whether 
the sphalerons go out of equilibrium above or below the electroweak phase transition [313]: 



Yab = r A Q x J 

a { 



24+4m rp rp 
66+13m 1 > J EWPT, 

32+4m rp rp y / 

98+13m 1 < ^EWPT, 

where m is the number of Higgs doublets. For the Standard Model (where m — 1), one obtains [41,313] 

[ 37 J EWPT- 

In ref. [315], the relation is given for any value of the Higgs VEV, valid through the EWPT. In the 
minimal supersymmetric Standard Model (where m = 2), one obtains 

= J* *;>w. (1418) 

I 31 1 < -'EWPT- 
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15 Appendix: Evolution of Flavoured Number Operators 



The aim of this appendix is twofold. First, to understand which flavour basis is appropriate for the 
Boltzmann Equations. Second, we investigate whether the physics that defines the preferred flavour basis 
introduces additional significant flavour-related effects on the baryon asymmetry. 

To address the first question, one would like to have a formalism which is covariant under flavour 
transformations, where it is possible to "derive" the Boltzmann Equations. One possible approach is 
to study Schwinger-Dyson equations for the non-equilibrium Green's functions of the particles involved 
in leptogenesis. The resulting Kadanoff-Baym [316] equations for Ni two-point functions were first 
discussed in [317]. Since the N% population must be out of equilibrium for leptogenesis, the aim of 
this paper was to study quantum non-equilibrium effects. The Boltzmann Equations were found to be 
a good approximation for non-relativistic A^. 29 Non-equilibrium Schwinger-Dyson equations were also 
used in [241,242] to show the importance of time-dependent CP violation in resonant leptogenesis. An 
elegant approach to flavour effects in leptogenesis would use the finite temperature Schwinger-Dyson 
equations for flavoured fermionic propagators, with a flavour off-diagonal e a p as a source for flavour off- 
diagonal propagators. An alternative approach, which should in principle be equivalent, is to study the 
equations of motion, in field theory, for a flavoured number operator (following, for instance, chapter 17 
of [319] . For a pedagogical introduction, see e.g. [320].) In this framework, however, neither the notation 
nor the formalism are transparent (similar to perturbative calculations without Feynman diagrams and 
rules). Here, we simplify this formalism and study a toy model of Simple Harmonic Oscillators (SHOs), 
then guess an extrapolation to field theory. We show that, in this toy model, fast interactions choose a 
basis where the equations of motion have the form of Boltzmann Equations. Thus, the flavour basis is 
determined by which interactions are fast at the time of leptogenesis. 

The Boltzmann equations for number densities treat interactions as a series of quantum processes of 
a classical particle. They neglect quantum effects, such as oscillations. In field theory, effects of quantum 
mechanical evolution appear in the equations of motion for the (flavour-dependent) number operator 

ffip) = <HP)4(P) - ^J(p)a a _(p), (15.1) 

which counts the asymmetry in lepton doublets. (We denote the operator / with a hat, to distinguish 
it from its expectation value.) The operator (o_ ) creates lepton doublets (anti-lepton doublets) of 
flavour a. Notice the inverted flavour order between the particle and antiparticle number operators [321]. 
The operator is a matrix in flavour space, analogous to the density matrix of quantum mechanics, and is 
sometimes called a density matrix. The diagonal elements are the flavour asymmetries stored in the lepton 
doublets. The trace, which is flavour-basis-independent, is the total lepton asymmetry. The off-diagonals 
elements encode (quantum) correlations between the different flavour asymmetries. 

Variants on such an operator have been studied in the context of neutrino oscillations in the early 
Universe [321-324], and, in particular, the generation of a lepton asymmetry by active-sterile oscillations 
[325-329] . The first discussion of such an operator in the context of thermal leptogenesis was made in 
ref. [65]. 

For simplicity, we replace f^{p) by the number operator of "flavoured" harmonic oscillators. Our 
aim is to separate the (possibly quantum) flavour effects from the (assumed classical) particle dynamics. 
That is, we extract the flavour structure from a toy model of harmonic oscillators, and input the particle 
dynamics and the Universe expansion, by analogy with the Boltzmann equations. The system contains 
harmonic oscillators of two "flavours", coupled to match the interactions of the Lagrangian in eqn. (|2.3p . 
By taking expectation values of the number operator in a "thermal" state, this model can be reduced to 
a two-state quantum system. We then extrapolate the flavour structure of these equations of motion to 
obtain "flavoured" equations of motion for the lepton asymmetry number operator in the early Universe. 

29 See however, [318] for numerical comparisons of the Kadanoff-Baym and Boltzmann Equations in a related model. 
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15.1 A toy model 



Consider a system of two simple harmonic oscillators with "flavour" labels 3 and 2, which have different 
frequencies in the free Hamiltonian: Hq = (w2o4 a 2+ Li; 3 a 3 a 3)^- Usually, the number operator is introduced 
as a\a,2 + 0303. However, since we are interested in changing the basis of the SHOs, we introduce the 
number operator as a matrix /: 



/ = 



a 3 a 2 



4 a 3 
4°3 



(15.2) 



whose trace is the usual number operator. The commutation relations 30 are [02, 4] = 1) [ fl 3>4] = !• 
The number operator evolves according to the Hamiltonian equations of motion: 4£ = +i[H, /]. Since 
our Hamiltonian conserves particle number, we can take expectation values - for instance, in a thermal 
bath - and reduce our toy model to a two-state system described by the amplitude of the particles to be 
3's, or 2's. The density matrix of this two-state system, is 



f = 



(a\a 2 ) 

(4°2> 



(4a 3 ) 
(4 a s) 



satisfying 



df 
dt 



(wa - 0J 3 )f 23 

(L0 3 -U 2 )f 32 



(15.3) 



If at t — the system is created in some state, then the probability to be found in this same state at a 
later time t is Tr{/(0)/(£)}, which can oscillate. 

It is useful to make an analogy to the use of this formalism in the context of neutrino oscillations [330] . 
We take the expectation value of the neutrino number operator in a beam of neutrinos of momentum 
p. This gives a density matrix / representing the flavour of the beam. For two neutrino species, this 
is a two-state quantum system. For an initial state \v(t — 0)) = \v^) = 023^2) + 523^3), (so that 
f(t = 0) = |^/j)(i / /i|), the density matrix at t = is given by 



/(* = 0) = 



c '23 
C23S23 



C23S23 

C 2 
L 23 



The solution to eqn. (|15.3[) with the initial condition (|15.4p is 

iAt 



fit) 



°'23 
C23S23C 4 



iAt 



C23S23e 

r 2 
c 23 



A 



2E 



We obtain the survival probability: 

P^{t) = Tk{/(0)/(t)} 



1 - sin 2 26*23 sin 2 At 



(15.4) 



(15.5) 



(15.6) 



Returning to lepton flavour in the early Universe, consider doublets £ a of momentum p 7 in a temper- 
ature range where the rate associated to the r Yukawa coupling is fast (IV ^ iJ), while the muon and 
electron Yukawa interactions can be neglected. The £ a 's have thermal masses [197,331]: 



m 2 (T) 



32 



n 

32 



16 



AT) 



32 



32 



(15.7) 



Notice the difference with respect to neutrino oscillations in the lab: In the early Universe, the mass 
basis is the flavour basis. The contribution of the neutrino Yukawa coupling to the thermal mass has 
been neglected in eqn. (|15.7[) . because we assume that h T |A a i|. For h T < |A Q i|, an additional term 
cx AqiA^/jvj should be included. This case is discussed in section [15.31 



30 In studying flavour, we consider "bosonic" leptons. We anticipate no difficulties in generalizing to fermionic leptons. 
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The two harmonic oscillators in the early Universe can be labeled r and o, where o is the projection 
onto e and /j, of the direction into which Ni decays: o oc X e ie + X^ip,. The eigenvalues of the free 
Hamiltonian are the particle energies u>(p), and the time evolution of the number operator is determined 
by the energy differences: 

h 2 T 2 

c^-^ = M_. (15.8) 

By analogy with the neutrino oscillation example, we anticipate that "flavour oscillations" might affect 
the lepton asymmetry. Imagine the excitations of the oscillators to be the lepton asymmetry, carried by 
particles of energies oj t ~ uj ~ E. Then uj t — uj ~ h 2 T 2 / (32E). We treat the production and washout 
of the lepton asymmetry as, respectively, initial condition and subsequent measurement on the system. 
In other words, we consider the production of an asymmetry in some linear combination of o and t, we 
allow it to evolve, and later turn on the inverse decays. Then the inverse decay rate can be suppressed 
or enhanced because the asymmetry changed flavour during time evolution (analogous to a survival 
probability in neutrino oscillations) [62]. Notice that the oscillation timescale is of order (/i^T) -1 , which 
is parametrically the same as the timescale of the /i r -mediated scattering rates that cause decoherence. 
This is different from, e.g. the MSW effect in matter, where the rate for decohering scattering can be 
neglected on the oscillation timescale. 

The flavour-blind gauge contribution is neglected in cqn. p5.8[) . This is subtle, because the early 
Universe version of A (see eqn. (U5.5[) ) depends on the energy of the lepton, and within the oscillation 
timescale, a lepton participates in many energy-changing gauge interactions. Ref. [62] used the thermally 
averaged energy (E) ~ 3T to estimate A [325]. In path integral language, this means approximating 
the integral (i J Edr) along the path from one lepton- number violating interaction to the next, to be 
i(E) J dr. This is in agreement with the analytic and numerical analysis of [332], which indicates that 
fast gauge interactions do not affect the coherence of flavour oscillations: if the timescale for transitions 
between different energy levels is much shorter than the oscillation timescale, then a particle spends 
time at many different energies during an oscillation timescale. The probability to have an energy E is 
proportional to the distribution function / of eqn. (|13.5|) . Therefore, on the oscillation timescale, the 
particles all have the thermal average energy, and oscillate coherently. In any case, this does not affect 
the principal claim of this appendix, that charged lepton Yukawa interactions choose the flavour basis 
for the Boltzmann Equations. 

We now include interactions among the SHOs. This allows one to take into account the production 
and annihilation of particles, and can lead to decoherence. Recall that in ordinary neutrino oscillations, 
decoherence happens anyway after a few oscillation lengths, due to "spreading of the wave packet" (see 
e.g. [333]). Here we are looking for decohering interactions, which, in the quantum mechanical analogy, 
collapse the wavefunctions onto an eigenbasis. We will see that in this basis, the equations of motion 
look like Boltzmann Equations. 

The decays and inverse decays due to the Yukawa interactions can be included by perturbing in an 
interaction Hamiltonian, iJj. We first consider the production and washout of leptons, due to A, in a 
model of four harmonic oscillators: N, <j>, and two leptons flavours {£ a }. Although we use the index a, 
this discussion is covariant in flavour space, so any basis choice for doublet flavour space is valid. We are 
looking for behaviors analogous to the lepton number production and washout in the decay/inverse decay 
N «-► (j>t a . We take </> and £ a to be massless and N massive, and ignore the free Hamiltonian (considered 
previously) which does not change the particle numbers. The interaction Hamiltonian is 

Hi = Xaia^a^a^ + X* al a N a^al a . (15.9) 
We perturbatively expand the Heisenberg equations of motion and obtain 



d_ 

= -Aai \* pl (Uff - f N US p0 - f N f! p - f N s Pp ) + (f N s ap + f N us ap + ]n!7 - Uft p )K^h 



frf" = -[Hi,[ H ii fe 13 ]] = -^pi^aiiaNa^al^alfCijatp] + A p i A^ [a^a^a^, ajyaj,aj 



~A Q i A* 



pi 



Uff {I + In) - Mi + UW P + ff) 
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+ Mi + UW + ft") - UffiX + In) x pl x* 01 



(15.10) 



This has the form of inverse decays and decays of N, with Bose enhancement factors (recall that the toy 
model leptons are bosons) for the final states. For simplicity, we now drop the Bose enhancement factors. 
A similar calculation can be performed for the tau Yukawa coupling h T (\h T \ = h T ). The interaction 
Hamiltonian is 

Hi.h = h T a i tj} a T oa tp + h* T p a rj> a\ c a\ p . (15.11) 

We assume that the processes 4> <-» I t t c are allowed, because the thermal mass of Higgs gets contributions 
from the top Yukawa interaction. (To avoid appealing to thermal masses, one could go to higher orders 
in Hi, and consider scattering such as W£ — ► T c <f>.) We neglect Bose enhancement factors and take 
expectation values. We obtain, in the charged lepton mass eigenstate basis, 



d_ 

Of 



f OO f 

Ji Ji 
Ji Ji 



or 

t 



-/r« 
+ 2/0 





\h T \ 2 



\h T \ 2 



£00 for 

Je Je 

£TO fTT 

Je u 



fro 



foo for 

Ji J I 

fro frr 

Je Je 





\hJ 2 



It is reassuring that the i-independent equilibrium solution of this equation imposes f T "fJ T 
Combined with eqn. (|15.3p . this gives damped oscillations for the off-diagonal elements: 

dff/dt = -i(w - w T - i\h T \ 2 f T c)f t T , 



(15.12) 

U = 0. 



(15.13) 



We conclude that exchanging doublet with singlet leptons destroys the coherence between different doublet 
flavours. In the "flavour" basis, the off-diagonal elements of the l a number operator vanish when the 
charged lepton Yukawa interactions are fast. The charged lepton Yukawa interactions are therefore 
the desired "basis-choosing" physics, that collapses the equations of motion of the density matrix to 
Boltzmann Equations for the number densities on the diagonal. 

It is convenient to write these equations in the notation of projectors, introduced in eqn. (|8.ip . In 
an arbitrary basis, we define, at tree level, the projectors onto the decay direction of N%, and onto the r 
flavour: 



-—a b 



P — = P ab = 



61 



h a h* 

pab ^r' l r / 15 w \ 

When loop contributions are included in the decays of iVi, a CP asymmetry can arise. It can be rep- 
resented by a difference between P and P, where P is the projector onto the anti-lepton direction into 
which N\ decays. 

Equations for the anti-lepton density can be derived in a similar way. For simplicity, we take (j) and 
N to be their own anti-particles, so there are now eight coupled oscillators: N,(f>,£ a ,£ a ,T c and f c . The 
expectation value of the difference between the two equations, again in the charged lepton mass eigenstate 
basis, is 

§- t fl = -\^\ 2 P ap Uf p A ~ \Xx?Ufl P P pP + 2(f N - f«)e<*P 



(15.15) 



where p,a, (3 = o, r. The purpose of this equation is to motivate the flavour index structure we use in 
the Boltzmann equations. We now briefly discuss the various terms. 

• The first two terms describe washout by inverse decays. By construction, there is no washout by 
AL = 2 scattering (we will introduce the resonant part by hand). Notice that the usual partial decay 
rates T(N — » <j>£ a ) become the diagonal elements of a matrix. The total decay rate, which is flavour-basis 
invariant, is the trace of this matrix. 
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• The third term, which is oc e, has been obtained artificially. It is well-known in field theory that 
no asymmetry is created at second order in A (see eqn. 15.191) . A CP asymmetry is proportional to the 
imaginary part of the loop amplitude times the imaginary part of the product of tree and loop coupling 
constants (see eqn. l5.5|) . We include the loop as an effective (non- unitary) interaction in the Hamiltonian: 
(<$A) Q f2 ajya^a^ + (<5A)* f2 0,^0^0^ . Here ft encodes the imaginary part of the loop amplitude; it is flavour 
independent, of little direct interest here, and can be obtained by matching to the calculation of e aa in 
section [5TT1 Just as the decay rate T becomes a matrix in flavour space, the CP asymmetry e also becomes 
a matrix in flavour space, proportional to [(5X) a Xp — X a (5X)%\ oc F Q/3 — P Q/3 . This gives 

e Q(3 oc Im{\ a (S\)* p - (SX) a X*p} = ^lm{X a m^ p X p - m ap X* p X*p}. (15.16) 
Using the equilibrium condition, 

2f N q [e} = ^([e][h + m + [h + fm) 

(where square brackets are matrices in flavour space), we obtain a production term 2(/jv + /^)e Q/3 - This 
cannot be complete, because it gives an asymmetry in thermal equilibrium. The resonant part of AL = 2 
scattering should be subtracted [31], which ensures that Tr /a vanishes in thermal equilibrium. Working 
in the charged lepton mass basis, it is straightforward to obtain eqn. (|15.15j) by following section \6?2\ or 
ref. [31]. 

• The remaining terms of eqn. (|15.15[) describe the effects of the tau Yukawa coupling. The asymmetry 
oscillates in flavour space, as described at the beginning of the section, because the off-diagonal flavour 
matrix elements have time-dependent phases. We drop this oscillatory term from the equation of motion, 
when extrapolating to the early Universe. (See [62] for a discussion including these oscillations). 

The tau Yukawa coupling also causes the off-diagonal elements to decay away, via the last two terms, 
This happens separately for the lepton and anti- lepton densities, as shown in eqn. (|15.13j) . However, 
the equation is not so simple to solve for the asymmetry, because the last term depends also on the 
asymmetry in the S'[/(2)-singlet r c s. The coupled equations for the singlet and doublet asymmetries 
should be solved. Instead, we attempt to include the singlet asymmetry via the A-matrix, as is done in 
the Boltzmann Equations, and argue in the following that the singlet asymmetry has little effect on the 
decay of the flavour-off-diagonal terms. 

The last two terms of eqn. (|15.15l) drive the diagonal asymmetry Y/\t T to the correct relation with 
Ya t , when the tau Yukawa coupling is in equilibrium. (In our toy model, where </> is a real scalar, this 
means (f T c — f fc ) = f^. In the SM, there is a Higgs asymmetry [41,65,313].) This can be clarified by 
studying the evolution of 

F°0 = f¥-P?0f AT ., (15.17) 

where P T is the projector onto the r direction in some arbitrary basis for lepton doublets, and /at c — 
f T o — /-<=. The trace of F is the total lepton number stored in ol, tl and tr, so studying F in the toy 
model is analogous to replacing the asymmetry in lepton doublets with the asymmetry in B/3 — L a [65] 
in usual leptogenesis calculations. For the remainder of this subsection, we focus on the decohering effects 
of the r Yukawa. This means we neglect the oscillation term and the neutrino Yukawa A. The equation 
for the doublet asymmetry is 

|[/A] = -\hr\ 2 {[Pr},[fA}} ifTC+ 2 hc) -M^PrWt + jf]} ^"'^ , (15-18) 

where square brackets denote matrices in doublet space. The equation for the singlet asymmetry is 

^/Ar= = -\h T \ 2 (f^+Mf T A T -\hr\ 2 fAAfi + fI) T \ (15-19) 



8<) 



where f T £ 



ft 1 • [/a] 
8 



dt 



[F] 



h. The equation for the lepton asymmetry F a/3 is 

(/r« + /r») 



|^| 2 
I ft 



({[P r ],[/A]}-2[P r ]Tr[P T /* 



r\ 2 ({[Pr], \h +//]}- mmpAh + m) (/rc 2 hc) 



(15.20) 



where the second line is of the same order as the first (see eqn. I15.15[) . The trace of the right hand side 
is zero as expected, since Tr(F) is conserved in the absence of A. In the flavour basis, eqn. (|15.20[) can 
be rewritten as 
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foo for 
J A J A 







J A 


dt 


£TO fTT f 

J A J A ~ /Ar c 
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fro 
J A 






(/t«+/pO 





fr + f? 



for j_ fro 

h + Jj 




(15.21) 

The first term causes to decay. The second term is of the same order, and can be of either sign. To see 
that, nevertheless, and /™ are driven to zero as h T comes into equilibrium, notice that eqn. (|15. 13[) 
implies that fj° + fj T vanishes exponentially fast. So the problematic second term vanishes as h T comes 
into equilibrium, and the equilibrium (time-independent) solution is = 0. 

In summary, the Yukawa coupling h causes the flavour off-diagonal asymmetries to decay away, as 
it did for the off-diagonal number densities in eqn. (|15.12p . It has no effect on the asymmetry in either 
flavour. To obtain eqn. (|15.22p . we simply drop the second line (so we neglect the asymmetry in SU{2)- 
singlets). 



15.2 Extrapolating to the early Universe 

A matrix equation describing the asymmetry AYe in the early Universe can be guessed by matching the 
flavour structure of eqn. (|15.15|) to the Boltzmann Equations: 



d_ 

dz 



1 A 



1 A 



1 A 1 A 



z 

7hI 



Y Nl 

vcq 



1 



2Y e eq 
1 



OO OT 

IN^2 IN^2 



1 A 



1 A 



1 A 1 A 



\roo 
1 Ai 



VOT 

1 Ai 



1 Ai 1 Ai 

, Tr[ 7r y A ] 

y, oq 



5 TT 



(15.22) 



where {•, •} stands for anti-commutator. We now discuss each of the three terms on the right hand side 
of eqn. p5.22[) in turn: 

1. The first term describes the creation of the asymmetry. The CP asymmetry in the a-flavour is e aa , 
and e = X)q £aQ - Notice that e a/3 is normalized by the total decay rate (so that [e] transforms as a 
tensor under i basis rotations). The matrix [e] is now defined as 



^Im{(A la )[AAt]yA^- (A^)[A*A T ] 1J A JQ } ff ^ 



(16tt) [AAt 
3 M± 

(167T-U 2 ) [\^t] 11 Im ^ A " Ag -[ m ^' T ' 3 ~~ K[ m A<y a Xf3\ (for Mi < M 2 . 3 ) 



where g(x) is the loop function of eqn. (|5.14p . 



(15.23) 



2. The second term is the washout of the asymmetry by decays and inverse decays. The matrix [77^2] 
is defined as 

AlaA*^ 



a/3 j-,r: > 

1 N ^2 = 7iV->2-P P = 1N-.2 



a/3 



E P iAi P r 



(15.24) 
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where 7jv^2 = ^2 p 1 P n^2 1S ^ e total (thermally averaged) decay rate. For a = f3, 7^^ 2 ^ s the decay 
rate of Ni to the a-flavour. This term is written as a function of the asymmetry in doublets, AYf , 
rather than the asymmetry in B — L: Y A — hfvY AT c . One way to proceed is to solve simultaneously 
this equation and the equation for 1at c - Alternatively, one can include the singlet asymmetry in an 
approximate way via an A-matrix, as one does in the Boltzmann equations. Since the /i T -interaction 
is entering equilibrium, the A-matrix is time-dependent: 

A al3 {t) = Af >e - r ^ + Af^l - e" 1 ^), (15.25) 

where A^Py (A^.) are the matrix elements before (after) h T comes into equilibrium, and T T = 
Yf q j(£ <-> tr) is the interaction rate associated to the tau Yukawa coupling, see eqn. (|14.10|) . We 
use this to obtain eqn. (|15.29|) . 

3. The last term describes the decay of the quantum correlations between o and r, that happens as 
the h T interaction comes into equilibrium and makes the r's distinguishable. In the flavour basis, 
we have 

" 



[7r] 



(15.26) 



lr 

where 7 r is the interaction density for the r Yukawa coupling (see discussion around eqn. I14.10jl : 

7r ~ 5 x lO- 3 h 2 T Tn £ q . (15.27) 
These interactions collapse the density matrix onto its diagonal elements in the flavour basis. 

Following Stodolsky [322,330], we write eqn. (|15.22[) in a more transparent way by expanding the 
various matrices on the basis provided by the u-matrices, = (/, a). A hermitian 2x2 matrix Y can 
be written as 

Y = Y»a^ Y^ = ^Tr{Ya^}. (15.28) 

The asymmetry matrix in the flavour basis is then replaced by a four-vector with components = 
(l/2)Tr [F A ], which represents half of the total lepton asymmetry, Y£ = (1/2) (Y^° - Y^ T ), which is half 
of the difference of the asymmetries in the o and r flavour densities, and Y A ' V , which account for the 
quantum correlations between the lepton asymmetries. 

In this notation, the Boltzmann equation p5.22[) becomes a system of equations of the form 



dY A _ z (Yj^ ^ [AY A ]° 1 _ 

r iVi 



s//| . I ^ - 1 1 e- 1D L Y^T- - yeqTjV-^ I^a] +7r x 7t x Y A ) . (15.29) 



Had we kept the thermal masses in our equations, we would find that the components Y£ and Y A 
precess around the z-direction with an angular velocity set by the thermal mass [62]. At the same time, 
such a precession is damped by the tau- Yukawa interactions at a rate ~ 7r , as described by the term 
7r x 7r x Y A . 

For 7t jn—>2, the asymmetry Ya is projected onto 7V before the washout interactions have time 
to act. So the last cross product term involving 7V can be dropped from the equations, provided that 
Ya II "TV is imposed. To summarize, we find that for "f T > 7 jv— >2j the baryon asymmetry can be obtained 
by solving the Boltzmann Equations in the flavour basis. 



15.3 Decoherence due to A 

In this section we consider the case where the N± decay rate is faster than, or of order of, the rate of r 
Yukawa interactions [180,181]. The evolution of the asymmetry density operator in this transition region 
was studied in [180], including the oscillations due to the thermal mass (115. 7|) . 



91 



z 


dz 


2sH 1 


dYl e 


z 


dz 


sHi 


dYl 


z 


dz 


sHi 


dYl 


z 


dz 



For 7at^2 > It, the distinguished direction in the doublet-lepton flavour space is £nh the direction 
into which Ni decays. Therefore the "single flavour" leptogenesis calculation applies, even if T T > H . To 
see this effect analytically, we use the orthonormal basis {Inii^p} °f the plane spanned by £n x and f (p 
stands for perpendicular: -£ p = 0). In this basis for the two-dimensional flavour space, the four- vector 
components are e° = e z — e/2, "f%^ 2 = 7jv^2 = 7jv^2/2, and r y% s ^ 2 = 7tv^2 = 0- The lepton asymmetry 
matrix Y^ e can be written as the four- vector Y^ eq (Ay°, Ayf,/S.y v l: Ay|). The equations (|15.29|) become 

" A - 7iV-*a(i/jv 1 -l)c-7^ 2 (Ay? + A«f) (15.30) 

: 77V^2(yw 1 -l)e-7w-2(Ay £ + Ay|) + 2( 7r xf r x Ayi) ■ i (15.31) 

: 7w^2 (y Wl - 1) e y - : ^Ay» + (> x 7V x Aft) • y (15.32) 

: 7^2 (y Nl -l)e x - ^y^Ay| + ( 7r xf r x Ay,) • £. (15.33) 

where the time variable z = M/T has exceptionally been represented as z, to avoid confusion with the 
four-vector index. In the above equations, the A-matrix has been neglected for two reasons. First, it 
includes the effects of other interactions, in the approximation that they are fast compared to N\ decays, 
while we are interested in the case where the interactions rates of Ni , tr and the sphalerons are of the 
same order. Second, it is simpler to consider the equations for the doublet asymmetry Y&£ rather than 
the B/3 — L a asymmetries Ya q . The ^4-matrix is included in [180]. 

The final baryon asymmetry is proportional to the total asymmetry in the doublets, which is the 
trace Y^ e = Y^ cq Ay°. The second washout term in eqn. Q15.30P depletes the asymmetry less effectively 
for Ay°+Ay| ~ 0. When Ay°+Ay| < 0, it contributes to generating the total asymmetry by reprocessing 
a flavour asymmetry. This change of sign must be accomplished by 7t , before the washout interactions 
have time to act. 

Consider an asymmetry produced in N\ decay. In a short timescale, before the £ have time to interact 
via A or h T , Ay" oc e a (where a = 0, x, y, z). Flavour effects can modify the evolution of the asymmetry 
if Ay| is determined by the h T interactions (recall that z — jn^2)- 

(7t x > x e) ■ 7at^ 2 = (TV • e It ■ 1n^2 - e°)l7r| > 7JV^2e°- (15.34) 

Notice that |e] ^> e° is possible, corresponding to e 00 ~ — e rr ^> e. Eqn. (|15.34p is satisfied for |e| 7T ^> 
7a , ^2£, provided that certain misalignment conditions hold. For the purpose of rough estimates, this 
condition can be taken to be 7 T 3> 7_/v^2, because eis washed out by 7 at -+2, so the approximation that 
Ay^ cx e 3 is only valid in a time-step shorter than the N% inverse decay time. The misalignment conditions 
are expected: if for instance N± decays only to rs (f • £n 1 = 1), the single-flavour analysis is correct. The 
misalignment is probably more transparent in flavour space than in 4-vector space. 



15.4 Summary 

The charged lepton Yukawa couplings are relevant for leptogenesis, because they can determine the initial 
state of washout interactions, and washout is crucial for thermal leptogenesis. 

In this appendix, we explicitly saw that the equations of motion of the "asymmetry number operator" , 
which are covariant in flavour space, are projected onto the flavour basis by the interactions of the 
charged lepton Yukawa couplings. In this flavour basis, the equations for the diagonal elements of the 
asymmetry operator become the familiar Boltzmann Equations, in which the charged lepton Yukawa 
interactions do not appear. We expect this projection to occur if the charged lepton Yukawa interactions 
are fast compared to leptogenesis timescales. Indeed, the flavour off-diagonal elements are negligible in 
the equations when 

It > 7jv^2 ■ (15.35) 
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Intuitively, this makes sense: if a lepton £ interacts many times via the h T Yukawa coupling before 
interacting with N±, it will participate in the washout process as a flavour eigenstate (£ r or £ a ). 

The temperature ranges where flavour effects should be included can be estimated by comparing 
rates, rather than rate densities. The timescale for leptogenesis is H , because the "non-equilibrium" 
is provided by the Universe expansion, so the Yukawa couplings can be neglected when they are out-of- 
cquilibrium, that is, slower than H. Comparing H to the rates for h T - or /i^-mcdiated interactions, such 
as q^R — * £tTRj one finds 

T Q ~ HT 2 h 2 a T > H for T < | ™" ^ 22^. ( 15 - 36 ) 

Below T ~ 10 12 GeV, h T is in equilibrium, and there can be two distinguishable flavours down to T ~ 10 9 
GeV. Below T ~ I0 9 GeV, is also in equilibrium, and there can be three. 

To impose the flavour basis in the washout processes, the charged lepton Yukawa interactions should 
also be faster than the inverse decays of N\. (Notice that the Boltzmann factor e~ Ml / T , heuristically 
included in Tjd, appears in the rate density 7jv-+2 of eqn. I|15.35|) .) We thus define a temperature Tfl, at 
which r r = Fid. Below Tjj, the lepton mass eigenstates are the flavour states. In flavoured leptogenesis, 
one can estimate a temperature Tb, after which the baryon asymmetry is generated. If Tb < Tk, then it 
is consistent to calculate the baryon asymmetry with flavoured Boltzmann Equations. 

With strong washout for all flavours, it is more convenient to keep track of the fraction of the N\ 
population remaining, oc e~ Ml ' T . At Tfl, when the thermal mass eigenstates become the flavour states, 
this fraction is ~ T T /Tn- The fraction remaining when a flavoured lepton asymmetry can survive is 
~ H/T aa , where T aa = min{r oo , T TT }. (T T is the /i T -mediated scattering rate, and T TT — T(Ni — > T<j)).) 
So leptogenesis can be calculated using BE written the flavour basis, when 

(15.37) 

J- D I- act 

(Recall that we are discussing the case of T aa > H .) For instance, for Mi ~ 10 10 GeV, and m ~ m atm , 
one finds that T T (Mi) > r£>(Mi), so the baryon asymmetry is generated after Tfl (for strong washout in 
all flavours). 

There are, however, regions of parameter space where Tb > Tfl [181], such as the strong washout 
regime at large M\, In this case, the baryon asymmetry could be divided into three parts: 

i) An unflavoured contribution from the decays taking place before Tm ~ IV . In strong washout and 
assuming T T > H, these are negligible. 

ii) A contribution from the transition region where Tijj ~ r T , to be obtained by solving the equations 
of the covariant asymmetry density operator. 

Hi) A flavoured contribution from the decays taking place after Tj£, ~ T T . 

In the limit where ii) is neglected, the final baryon asymmetry can be estimated analytically: It is 
the flavoured asymmetry that is produced in the decays of the Aq's that remained at Tfl. See [180] for a 
detailed discussion based on solving the equations for the asymmetry density operator. 
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16 Appendix: Analytic Approximations to Yab 



In this Appendix, wc reproduce analytic estimates of the baryon asymmetry, derived in ref. [64] . These are 
approximate solutions of the BE which take into account decays, inverse decays, and AL = 1 scatterings 
with Ni on an external leg. CP violation in all these processes is included. The off-diagonal elements of 
the A-matrix are neglected (they are included in Section [973]) . as is AL = 2 scattering, so the equations 
for the B/3 — L a asymmetries are decoupled. This simplifies the solutions. Ref. [64] estimates the overall 
uncertainty related to the various approximations to be of order 30%. 

A detailed discussion of analytic approximations (as well as numerical evaluations) can be found 
in ref. [59], for single-flavour leptogenesis and without CP violation in AL = 1 scatterings. We here 
modify the formulae of [59] to include flavour and CP violation in scattering, providing simple and useful 
approximations . 

The Boltzmann equations can be written as follows: 



dz 
d 

dz 



S{z), 



W aa (z)Y Aa 



(16.1) 
(16.2) 



where the source and washout functions can be read from eqns (|6.54|) and (|6.56|) . Their functional form 
can be found in ref. [59], or, with finite temperature effects, in ref. [60]. 

As discussed earlier, the values of the light neutrino mass-squared differences suggest that fh > m*. 
In this strong washout regime, there is a useful approximation to Yn ± (z) — Y^(z) at large z (see chapter 
6.4 of ref. [310]). Indeed, the asymmetry is produced when the inverse decays go out of equilibrium, 
which happens at z 3> 1. The relevant equation is then (j6TT6j) : 



dz 



Y, 



Z 

7th 



Yi 



V cc l 



1N- 



V c 1 



dz Wl 



(16.3) 



When the interaction rate jn^2/(sY^) is fast compared to H, the differential equation is "stiff" and 
can be solved approximately by setting the right-hand-side to zero, and taking dY^/dz 



Yj 



1 N l 



sYfi H lycq _ zK 2 {z)H 1 
1N^2 z Nl 4g*r_o 



(16.4) 



The last equality uses eqn (|13.2[) and, from eqn (|13.17p . / (sY^ 1 ) — Fn- 

Using the trick of rewriting eqn (|16.2[) as (fY)' — fSe aa , with /' = W aa f, it has the formal solution 



Y Aa {z) 



dxS(x)e-£ dyW "« (y) . 



(16.5) 



We are interested in Y Aa {z — > oo). This can be approximated via the "steepest descent", or "saddle- 
point" evaluation of an integral along the contour C [334] : 



J(x) 



„xF(t: 



dt 



2tt 



-xF"(t ) 



cF(t ) 



(16.6) 



where t is an extremum of the integrand. In the case of eqn (|16.2[) . the extremum is at W aa — 
d/dz(ln \S\). Using the saddle point approximation, and simple forms for the functions S and W aa , 
ref. [64] obtained the following approximations for the efficiency parameter : 



| -A aa | Tfl aa 

2.1m* 

\A aa \ih aa fh 
2.5m* m, 



2 A 



aa " L aa 



-1.16 



(fh > m») 
(m < m*) 



(16.7) 
(16.8) 
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where to* and 

Tfiaa ftxc defined in eqn (|4.7p , and the -A QQ &re defined in ecm (|9. 1 8|) . The B/3 — L a 
asymmetries are given by 

Y Aa = 4 x 10- 3 e aQ 77 Q (16.9) 
and finally the baryon asymmetry (in the SM) is 

Yab^ ^E 1 ^- (16.10) 

a 

The temperature range where tau Yukawa interactions are faster than the expansion rate H is T < 
(1 + tan 2 (3) x 10 12 GeV (tan/? = 1 in the SM). If the tau Yukawa interactions are also faster than Tjrj 
when the asymmetries are generated (see Appendix 115.41) , then there are two relevant CP asymmetries - 
£oo = £ee + and e TT - and two relevant partial decay rates - rh 00 — rh ee + m MAt and m TT . So the sum 
in eqn (|16.10[) is over a = o,t and one should use eqn (|14.14[) for the A aa . For T < (1 + tan 2 (3) x 10 9 
GeV, the muon Yukawa interactions are also faster than H and there are three distinguishable flavours. 
The sum in eqn Q16.10P is over a — e, /i, r, and one should use eqn (|14.15[) for the A aa . 
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